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The application of random-matrix theory (RMT) to compound-nucleus (CN) reactions is re- 
viewed. An introduction into the basic concepts of nuclear scattering theory is followed by a sur- 
vey of phenomenological approaches to CN scattering. The implementation of a random-matrix 
approach into scattering theory leads to a statistical theory of CN reactions. Since RMT applies 
generically to chaotic quantum systems, that theory is, at the same time, a generic theory of 
quantum chaotic scattering. It uses a minimum of input parameters (average S-matrix and mean 
level spacing of the CN). Predictions of the theory are derived with the help of field-theoretical 
methods adapted from condensed-matter physics and compared with those of phenomenological 
approaches. Thorough tests of the theory are reviewed, as are applications in nuclear physics, with 
special attention given to violation of symmetries (isospin, parity) and time— reversal invariance. 
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I. INTRODUCTION 



Random-matrix theory (RMT) as developed in the 
1950s by Wigner an d Dyson (see ( Dvsonl . Il962allbl 3: 
Wignerl . [l955al [l957t ) and the reprint collection ( Porteil 
19651)) plays an important role not only in the anal- 
ysis of nuclear spectra. Random matrices and chaos 
play perhaps an even bigger role in the theory of nu- 
clear reactions. The resulting "Statistical Theory of 
Nuclear Reactions" is the topic of the present review. 
It completes the review by two of the present au- 
thors o£^aiidflm_jnAtrices_a^ nuclear struc- 
ture (jWeidenmuller and Mitchelil . l2009t) . 

In proposing RMT, Wigne r was probably inspired 
by Bohr's idea ( Bohrl . Il936| ) of the compound nu- 
cleus (CN). After the early experimental confirmation 
of Bohr's idea as implemented in th e Hauser-Feshbach 
formula (jllauser and Feshbachl . Il952( ) , th e field received 
a boost in 1960 by Ericson's prediction ( Ericsonl Il960l ) 
of statistical fiuctuations in nuclear cross sections. Sub- 
sequent intense experimental work on a number of topics 
in CN reactions (Ericson fluctuations, isobaric analogue 
resonances, isospin mixing in nuclear reactions, tests of 
time-reversal symmetry) reached saturation at the end 
of the 1970s, to be followed later only by studies of parity 
violation in nuclear reactions. Theoretical work extended 
the H auser-Feshba c h form ula to the case of direct reac- 
tions ( Kawai et all [l973l ) . At the same time, theorists 
set out to connect the statistical models of CN scatter- 
ing with RMT. That turned out to be a very challenging 
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problem. Motivated by the fundamental interest in CN 
scattering and by the need in other flelds of physics (neu- 
tron physics, shielding problems, nuclear astrophysics, 
etc.) to have a viable theory of CN reactions with pre- 
dictive power, that work was carried on for a number 
of years and led to partial insights into the connection 
between RMT and CN scattering. 

Theoretical efforts at constructing a comprehensive 
theory of CN reactions received a strong stimulus in the 
beginning of the 1980s by developments in the theory of 
chaotic motion, i.e., the theory of non-integrable classi- 
cal systems and their quantum counterparts. It is not 
a coinci dence that th e Boh igas-Giannoni-Schmit con- 
jecture ( Bohigas et aLl . [l984l ) (which connects properties 
of quantum spectra of class ically chaotic systems with 
RMT) and the first pape rs (|Verbaarschot et al I ll984aL 
Il985l : IWeidenmiilleil . [l98l establishing a firm connection 
between RMT, CN scattering theory, and chaotic quan- 
tum scattering, appeared almost simultaneously. Techni- 
cally, progress became possible by combining scattering 
theory based on the shell model with novel techniques 
using a supersymmetric generati ng func t ional borrowed 
from condensed-matter physics ( Efetovl . Il983() . Physi- 
cally, CN scattering was recognized as a paradigmatic 
case of chaotic quantum scattering, and the theory of 
CN scattering was seen to apply to chaotic scattering 
processes in general. Actually the theory of CN scatter- 
ing is richer than that for most other cases of chaotic 
scattering, for two reasons. First there exist conserved 
quantum numbers (spin and parity). The nuclear cross 
section is the square of a sum of resonance contributions 
each carrying these quantum numbers. Chaotic motion 
only affects resonances carrying the same quantum num- 
bers. Second, the nucleus has internal structure. This 
leads to a strong increase of the number of open chan- 
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nels with excitation energy. Inelastic processes (where 
the masses, charges and/or excitation energies of the re- 
action products differ from those of target and projec- 
tile) add complexity and richness to the theory: Nuclear 
reactions can be studied versus scattering angle, versus 
bombarding energy, and for different final fragment con- 
figurations. 

The statistical theory of CN reactions is generic and 
applies likewise to many other cases of chaotic scattering. 
That fact is borne out by applications to elect ron trans- 
port t hrough dis o rdered meso s copic samples ( Alhassidl . 
I2OOOI: iBeenakkeii Il997l: llmrvl . |2002| ). and to the pas- 
sage of electromagnet i c wav es through microwave cav- 
ities ( Fvodorov et all |2005() . Some of these cases are 
treated below. However, the theoretical developments 
reviewed in this paper do not cover all aspects of the 
theory of chaotic scattering. In systems with fe w degrees 
of freedom, semiclassical periodic-orbit th eory ( Gaspardl . 
I1991I: iGutzwilleil . I1990I: ISmilanskvl . Il99ll ) plays a promi- 
nent role. That branch of scattering theory has not been 
much used in nuclear many-body physics. (It has found 
applications, for instance, in the scattering of two heavy 
ions where it applies to relative motion). That is why 
it is not dealt with here. The connection between the 
RMT approach and periodic-orbit theory is discus sed, 
for instance, in ( Lewenkopf and Weidenmiillerl . [l99ll) . 

In view of the very general applicability of the statis- 
tical theory to chaotic scattering processes governed by 
RMT, we aim at a presentation which is accessible to 
readers not familiar with the topic. We have in mind, 
for instance, physicists working in other areas of nuclear 
physics, or in chaotic scattering. Therefore, we begin 
this review in Section |ll] with a summary of some basic 
facts and concepts of nuclear reaction theory. We also 
present the central ideas and models that were developed 
with the help of plausible albeit intuitive arguments and 
that were used to treat CN scattering before the con- 
nection to RMT was established. These are the Hauser- 
Feshbach formula, the Weisskopf estimate for the average 
total width of CN resonances, Ericson fluctuations, and 
modifications of the Hauser-Feshbach formula due to di- 
rect reactions. In order to avoid repeating our arguments 
later in slightly different form, we present the arguments 
in modern terminology. 

The modern access to the statistical theory of CN reac- 
tions is based on RMT and presented in Sections IIIIl IIVI 
IVl While the aforc-mcntioncd models do not all survive 
close scrutiny, at least their results are vindicated and 
their ranges of validity are established. Needless to say, 
additional results are also obtained. Tests of the theory 
and applications to a number of topics are reviewed in 
later sections of the paper. 

The field has not been reviewed comprehensively 
for man y years. A rev iew of RMT in nuclear 
physics ( Brodv et all Il98l[) contains sections on the 
statistical theory of nuclear reactions. Shorter re- 
views may be found i n dBohigas and Weidenmiiller ^ 198S ; 
iFvodorov and Savinl I2OIOI: iMahaux and Weidenmiilleii. 



Il979l : IWeidenmulleil.l2002l). We refer to Part I of this re- 
view (IWeideumfiller and Mitchell . 120091) with the letter 
"I" so that equations, figures or Sections in that paper 
are referred to, for instance, as Eq. (1.34), as Fig. 1.16, 
or as Section I. II. A. As in part I, we have preferred cit- 
ing a review over giving a large number of references: 
Readability of the article was our primary concern. 



II. BASIC FACTS AND CONCEPTS. STATISTICAL 
MODELS 

A stable nucleus with mass number A possesses a dis- 
crete spectrum of levels that extends from the ground 
state up to the lowest energy where decay by particle 
emission is possible (the first "particle threshold" ) . (Here 
we disregard the small widths of levels due to beta or 
gamma decay). In most cases the first particle thresh- 
old corresponds to nucleon emission (A — > {A — 1) + n) 
and typically has an excitation energy of 6 or 8 McV 
in nucleus A. The levels above that threshold have fi- 
nite widths for particle decay and appear as resonances 
in the scattering of a nucleon by the nucleus with mass 
number {A — 1). The density p{E) of nuclear levels in- 
creases roughly exponentially with excitation energy E 
(more precisely, p{E) ex exp V aE where a is a mass- 
dependent constant) and the average spacing d = 
of resonances decreases accordingly. The number of de- 
cay channels also grows with E since the density of states 
available for decay in neighboring nuclei likewise grows 
nearly exponentially. As a result, the average total decay 
width r of the resonances grows strongly with excitation 
energy. The Weisskopf estimate given in Eq. ([TT|) below 
shows that T/d is roughly given by the number of decay 
channels over 27r. Thus, we deal with isolated resonances 
(r ^ d) at the first particle threshold and with strongly 
overlapping resonances (F ^ d) several MeV higher. A 
comprehensive theory of nuclear reactions should cover 
the entire range from F ^ c? to F ^ d. 



A. Resonances 

Resonances in the cross section play a central role in 
the theory. In Part I, the empirical evidence was dis- 
cussed showing that isolated resonances measured near 
neutron threshold or near the Coulomb barrier for pro- 
tons display stochastic behavior: Spacings and widths of 
resonances with identical quantum numbers (spin, par- 
ity) are in agreement with predictions of RMT, more pre- 
cisely, with the predictions of the Gaussian Orthogonal 
Ensemble of Random Matrices (GOE). That evidence is 
here taken for granted and not reviewed again. The con- 
nection between RMT and chaotic motion was also re- 
viewed in Part I. We will use the term "chaos" as syn- 
onymous with spectral fluctuation properties of the GOE 
type. The theory of nuclear reactions makes use of the 
GOE properties of nuclear resonances. It is postulated 
that the stochastic features found for isolated resonances 
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also prevail at higher bombarding energies where reso- 
nances overlap. The stochastic description of resonances 
then applies for all bombarding energies where resonance 
scattering is relevant. Actually, statistical concepts are 
used in one form or another to describe all collision pro- 
cesses between atomic nuclei with center~of-mass ener- 
gies between and about 100 McV, except for reactions 
between pairs of very light nuclei (where the density of 
resonances is too small for a statistical approach). 

A stochastic description of CN resonances is not only 
physically motivated, but a practical necessity. While 
typical spacings of neighboring levels near the ground 
state are of the order of 100 keV, the resonances seen in 
the scattering of slow neutrons have typical spacings of 
10 eV, see Fig. I.l. This is a consequence of the nearly 
exponential growth of the average level density with E 
mentioned above. Put differently, there are about 10^ 
to 10^ levels between the ground state of the CN and 
the isolated resonances seen in the scattering of a slow 
neutron. There is no viable theoretical approach that 
would allow the prediction of spectroscopic properties of 
such highly excited states. Needless to say, the situation 
becomes worse as the excitation energy increases further. 

Within the framework of a statistical approach, one 
does not predict positions and widths of individual res- 
onances. Rather, the GOE predicts the distribution of 
spacings between resonances, and the distribution of par- 
tial and total widths for decay into the available channels. 
By the same token, the statistical theory of nuclear re- 
actions does not aim at predicting the precise form of 
some reaction cross section versus energy or scattering 
angle. Rather, it aims at predicting the average values, 
higher moments, and correlation functions of cross sec- 
tions obtained by averaging over some energy interval. 
That interval must encompass a large number N of res- 
onances. For isolated resonances, the resulting finite- 
range-of-data error is expected to be inversely propor- 
tional to N. For overlapping resonances, N is replaced 
by the length of the averaging interval divided by the 
average width of the resonances. With the exception of 
the lightest nuclei, such averages are the only theoreti- 
cal predictions presently available for nuclear reactions 
in the regime of resonance scattering. 



Dynamic Origin of Resonances 

It is useful to have an idea how the numerous res- 
onances dominating CN reactions come about dynami- 
cally. Subsequent theoretical developments then do not 
appear as purely formal exercises. We use the nuclear 
shell model reviewed in Section I.IV.A as the funda- 
mental dynamical nuclear model. In that model, two 
mechanisms lead to resonances. The dominant mecha- 
nism is that of formation of bound single-particle shell- 
model configurations with energies above the first par- 
ticle threshold. These states become particle-unstable 
when the residual two-body interaction is taken into ac- 



count. We refer to such states as to quasibound states. 
iMahaux and Weidenmiillerl ( 19691 ) use the term "bound 
states embedded in the continuum" (BSEC). Details fol- 
low in the next paragraph. A second and less important 
mechanism is due to barrier effects of the shell-model 
potential. The angular-momentum and Coulomb barri- 
ers cause the occurrence of more or less narrow single- 
particle resonances. These occur typically within the 
first MeV or so of the continuous spectrum of the single- 
particle Hamiltonian. There is at most one such narrow 
single-particle resonance for each angular-momentum 
value; resonances at higher energies are too wide to mat- 
ter in our context. It is clear that with only one single- 
particle resonance for each value of angular momentum 
we cannot account for the numerous CN resonances with 
equal spins and average spacings of about 10 eV ob- 
served at neutron threshold. Therefore, the quasibound 
states are the main contributor to the large number of 
resonances in CN scattering. Single-particle resonances 
can be incorporated i n the description of CN resonances 
as quasibound states ( Mahaux and Weidenmiille"il . Il969l ) 
and are not mentioned explicitly again. 

To describe the quasibound states in more detail, we 
consider nuclei in the middle of the sd-shell, i.e., nuclei 
with mass number 28 and with 12 valence nucleons. That 
same example was extensively discussed in Part I. Our 
considerations apply likewise to other nuclei. The spac- 
ings of the energies of the single-particle states (1^/2, Si/2j 
and (^3/2 (in spectroscopic notation) in the sd-shell are 
of the order of MeV. In ^^O, for instance, the spacing 
between the lowest (^5/2) and the highest ( d^/?) single- 
particle state is 5.08 MeV (jZelevinskv aUllEIQd) . Thus 
the spectrum of the bound sd-shell-model configurations 
of 12 nucleons (which typically comprises 10^ states for 
low values of total spin J) extends over several 10 MeV 
while the first particle threshold in these nuclei has a typ- 
ical energy of only several MeV. Inclusion of the residual 
interaction within this set of bound states spreads the 
spectrum further and produces chaos. But the two-body 
interaction also connects the resulting bound many-body 
states with other shell-model states defining the open 
channels. Such states are obtained, for instance, by lift- 
ing one nucleon into the s-wave continuum of the shell 
model and diagonalizing the residual interaction among 
the remaining 11 nucleons in the sd-shell. The antisym- 
metrized product of the ground state or of the nth ex- 
cited state containing 11 nucleons with the single-particle 
s-wave continuum state would define channels with dif- 
ferent threshold energies. The coupling of the diagonal- 
ized quasibound 12-nucleon states to the channels causes 
the states above the first particle threshold to turn into 
CN resonances. Lifting one or several nucleons out of 
the sd-shell into bound states of the p/-shell increases 
the number of bound single-particle configurations and, 
thus, of CN resonances. This shows how the (nearly ex- 
ponential) increase of the nuclear level density leads to a 
corresponding increase of the density of CN resonances. 
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B. S-Matrix 

We need to introduce some elements of nuclear reaction 
theory. We assume throughout that two-body fragmen- 
tation dominates the reaction. (Three-body fragmenta- 
tion is a rare event in the energy range under consider- 
ation) . Mass and charge and the internal states of both 
fragments are jointly referred to by Greek letters a, /3, . . .. 
The energy of relative motion of the fragmentation a at 
asymptotic distance is denoted by Ea while E denotes 
the excitation energy of the CN. By energy conservation, 
Ea is trivially related to E and to the energy Ep of rela- 
tive motion of any other fragmentation j3 via the binding 
energies of the fragments. Channels a,b, . . . are speci- 
fied by the fragmentation a, by the angular momentum 
£ of relative motion, and by total spin J and parity H. 
A channel a "opens" when Ea = 0, the corresponding 
value of E is called the threshold energy for fragmen- 
tation a. A channel is said to be open (closed) when 
Ea > {Ea < 0), respectively. We speak of elastic scat- 
tering when incident and outgoing fragments are equal 
{a = (3) because then Ea = Ep. But elastic scattering 
does not imply a = b (while J is conserved, £ is not when 
the fragments carry spin). We use the term "strictly elas- 
tic" for a = b. Some of these concepts are illustrated in 
Fig.ffl 
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FIG. 1 Level scheme of O. The center part of the figure 
shows the low-lying levels, beginning with the ground state. 
On the left and on the right, various thresholds are indicated, 
plus t he energy depende nce of some reaction cross sections. 
From (jTillev et aU 1 19951 ). 

The central object in the theory of nuclear reactions is 
the scattering matrix S{E). It is a matrix in the space 



of open channels and carries fixed quantum numbers J, 
n and, if applicable, isospin T. It is defined in terms of 
the asymptotic behavior of the scattering wave functions 
1 4'+). These are solutions of the Schrodinger equation 
subject to the boundary condition that there is an incom- 
ing wave only in channel a. We simplify the notation by 
considering only channels with neutral fragments. With 
£c the angular momentum of relative motion in channel 
c, Tc the radial coordinate of relative motion, kc the wave 
number, rric the reduced mass, and the spherical Han- 
kel functions, the radial part of in channel b has the 
form 5ab{ma/k'J'^)h~^ {kaTa)- {mb/kl^'^)Sba{E)hj^ {kbn). 
Thus, the element Sba{E) of S gives the amplitude of the 
asymptotic flux in channel b for unit incident flux in chan- 
nel a at energy E. The straightforward generalization to 
charged fragments is obtained by replacing the spherical 
Hankel functions by the Coulomb wave functions. The 
deflnition of a channel includes the quantum numbers J 
and n. In writing SabiE) we omit J and 11 (which must 
be identical for a and b). In general, channels a and b will 
belong to different fragmentations a and /3, respectively. 

The nuclear Hamiltonian is invariant under time re- 
versal. Therefore, the amplitude of the asymptotic flux 
in channel b for unit incident flux in channel a is equal 
to the amplitude of the asymptotic flux in channel a 
for unit incident flux in channel 6, and S is symmetric, 
Sab{E) = Sba{E). The conservation of total flux implies 
the unitarity relation 

Y,SaciE)S;,iE)^Sab. (1) 

c 

The dimension of S is equal to the number A of open 
channels. With increasing E that number grows nearly 
exponentially because so does the number of states in 
the residual nuclei into which the CN may decay. We 
can work with a matrix S of fixed dimension A only for 
E in an energy window with end points given by two 
nearest threshold energies. The size of that window de- 
creases nearly exponentially with E. That constraint is 
usually neglected in applications of the theory because 
the fiux into channels that have just opened is typically 
small due to angular-momentum and Coulomb-barrier 
effects. Moreover, because of the nearly exponential in- 
crease of the level density the number A/" of resonances 
within each window changes only algebraically (and not 
exponentially) with energy E, in spite of the exponen- 
tially decreasing size of the window. At neutron thresh- 
old, where d is of the order of 10 eV and the spacing 
of thresholds is of the order of several 100 keV, we have 
TV w 10^. At 15 MeV above neutron threshold (the Er- 
icson regime, see Section riI.D|) M is estimated to be one 
or two orders of magnitude smaller. 

Given S, the differential cross section of any nuclear 
reaction is obtained bilinear form in the elements 
SabiE) of the S-matrix, and similarly for other observ- 
ables such as fragment polarization. In addition to the 
elements of the 5'-matrix, the formulas involve kine- 
matical factors and angular-momentum coupling coef- 
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ficicnts ( Blatt and Biedenharnl . Il952[) and are not repro- 
duced here. The statistical theory of nuclear reactions 
focuses on the calculation of the statistical properties of 
the iS-matrix. 

For later use wc introduce the eigenvalues of the uni- 
tary and symmetric matrix S. These have unit mag- 
nitude and are written as exp{2iSc{E)). The real phase 
shifts Sc{E) are called the eigenphase shifts of S. A theo- 
rem analogous to the Wigner-von Neumann non-crossing 
theorem for the eigenvalue s of a Hermitean ma trix also 
applies in the present case ( Weidenmiilleil . Il967t) . There- 
fore, generically no two eigenphase shifts coincide, and 
the matrix 5* can be written as 



Sab{E) = ;^ac(i^)cxp(2z<5,(i;))a,(i;) 



(2) 



Here the real orthogonal matrix 0{E) induces a transfor- 
mation from the physical channels to the "eigenchannels" 
of S. The scattering wave functions in the eigenchannel 
representation are given by \{lc{E)) = ^a^ca\'^t{E)). 
The radial part of \ilc{E)) in the physical channel a has 
the form Oca[hg^{kara) - exp{2iSc)h^^{kara)], in keeping 
with the fact that S is diagonal in the eigenchannel basis. 
The eigenchannels are unphysical because there is an in- 
cident wave in every physical channel. Nevertheless the 
eigenchannels are helpful theoretical constructs. 



C. Bohr Assumption and Weisskopf Estimate 

Statistical concepts have governed the th eory of CN 
scattering from its inception. iBohrl ( 19361 ) introduced 
the idea of the CN as an equilibrated system of strongly 
interacting nucleons. The incident nuclcon shares its en- 
ergy with the nucleons in the target. The system equi- 
librates and "forgets" its mode of formation. It takes a 
long time (long in comparison with the time it takes a 
nucleon with Fermi velocity to traverse the nucleus) for 
the CN to accidentally concentrate the available energy 
back onto a single nucleon which can then be re-emitted. 
Therefore, formation and decay of the CN are indepen- 
dent processes ("Bohr assumption"). The decay of the 
CN is assumed to be governed by statistical laws (with 
the proviso that energy, spin, and parity are conserved). 

That intuitive picture found its first quanti- 
tative formulation in the "Hauser-Feshbac h for - 
mula" (|Hauser and Feshbachl . Il952t IWolfensteinl . Il95l[ ) 
for the average differential cross section. The average 
is taken over an energy interval containing a large 
number of resonances. We list the assumptions that 
were used and defer a discussion of their validity to 
Section |Vl Because of the presence of resonances (which 
behave stochastically) the scattering matrix S fluctuates 
randomly in energy and is accordingly decomposed into 
an av erage part and a fluctuating part (jFeshbach et all 
[1951, 



The average over energy is indicated by angular brack- 
ets. By defi nition, we have (S^h(E )) = 0. The standard 
assumption ( Feshbach et all Il954[ ) on the average part 
is that it vanishes for a b, 



(SabiE)) = 5ab{Saa{E)) 



(4) 



It is also assumed that in the CN, S-matrix elements 
pertaining to different conserved quantum numbers are 
uncorrelated, 

{S^,iE)iS^,iE)r) = for J^.r and/or U ^ W (5) 

and that even when the quantum numbers are equal we 
have 



{S^,iE)iS!,{E)r)=0 

unless a — c, b ~ d or a = d, b = c 



(6) 



where we have used the symmetry of S. Intuitively, the 
assumptions ^ and ^ are related to random phases for 
the contributing resonances (only absolute squares sur- 
vive the averaging process). The decomposition (|3]) of 
the scattering matrix implies a corresponding decompo- 
sition of the average cross section. The average consists 
of a sum over terms of the form {SabiE)S*^{E)) . We use 
Eq. (U) and focus attention on the part which is bilinear 
in the elements S^. It is that part which contributes to 
the average CN cross section (a™). With the assump- 
tion (0), (crCN^ contains only terms (^a&(^)('5'cd(^))*) 
where all channel indices refer to the same quantum 
numbers J and 11. That statement implies that av- 
erage CN cross sections are symmetric about 90 de- 
grees in the center-of-mass (cm.) system. The terms 
{S^bi^)i^cd(.E))*) are subject to assumption ^ which 
reduces the average CN cross section to a sum over av- 
erages of squares of iS-matrix elements. Each such term 
(|S'^^(£')P) describes the formation of the CN from chan- 
nel a and its decay into channel b, or vice versa. The Bohr 
assumption (independence of formation and decay of the 
CN) is used to write {\S^f^{E)\'^) in factorized form as 

{S^,{E){S^,{E)r) - (SacSbd + SadSbc) Tafb ■ (7) 

Here Ta denotes the probability of formation of the CN 
from channel a and is defined by 



Ta = l-\{Saa{E))\'' . 



(8) 



SabiE) = (SabiE)) + S^,iE) 



(3) 



The factor fb gives the relative probability of CN decay 
into channel b and is normalized, J2b fb ~ 1- Using the 
symmetry and unitarity of S and neglecting a term that 
is inversely proportional to the number of channels, we 
obtain from Eqs. ^ and ^ that /fa = TblY.c ^c- Thus, 

{S^,iE)iS^^iE)r) ^ (SacSbd + SadSbc) TaTb/Y^T, (9) 



while the average of the CN part of the cross section takes 
the form (we recall that a, /? denote the fragmentation 
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while a,b, . . . denote the channels) 
(da^^/duj) = ^ coefficients (^(l+,5a6)T,Tb/^r, 



xP^(cos0) 



(10) 



Here w is the solid angle, Q the scattering angle, and Pi 
the Legendre polynomial. The sum extends over all val- 
ues of J,£,a^b and contains geo metric and kinematical 
coeffic ients not specified here, see lBlatt and Biedenharr] 
(|1952[ ). Except for the fac tor (1 + Sg b) . Eg (fTO]) 
w as originally proposed by IWolfensteinI ( 195l[ ) and 
bv lHauser and FeshbachI ( 19521 ) and is commonly referred 
to as "Hauser-Fe shbach formu la" . The factor (1 + 6ab) 
was later shown ( Vagerl . Il97ll) to be a necessary conse- 
quence of the symmetry of S and for obvious reasons is 
referred to as the "elastic enhancement factor" . In the se- 
quel we apply the expression "Hauser-Feshbach formula" 
to both, Eq. © and Eq. ([TU| . 

Because of the unitarity of S, the average S'-matrix 
is subunitary. The "transmission coefficients" Ta defined 
in Eq. ([5]) measure the unitarity deficit of (5) (we recall 
Eq. dH)). The Ts obey < < 1. It is natural to 
interpret the unitarity deficit as the probability of CN 
formation (or, by detailed balance, of CN decay) from 
(into) channel a. We speak of weak (strong) absorp- 
tion in channel a when Ta is close to zero (to unity), 
respectively. The Ts arc central el ements of the theory . 
In th e early years of CN theory ( Blatt and Weisskopj . 

the CN was assumed to be a black box, so that 
Ta = 1 (Ta = 0) for all channels with angular momenta 
( < 4iax {£ > ^max, respectively). Here ^max is the an- 
gular momentum corresponding to a grazing collision be- 
tween both fragments. According to Eq. (jH)), the assump- 
tion Tc = 1 implies (5) = 0, and the decomposition 
was in fact only introduced when Feshbach, Porter and 
Weisskopf ( Feshbach et all Il954l ) proposed the optical 
model of elastic scattering. The model was originally 
formulated for neutrons but soon extended to other pro- 
jectiles. The model changed the view of CN reactions: 
The CN was not a black box but was partly transpar- 
ent. The transmission coefficients Ta in Eq. (|8]) were not 
put equal to unity but could be calculated from the op- 
tical model which provided the first dynamical input for 
CN theory (aside from the average level density that is 
needed to calculate the energy dependence of CN emis- 
sion products). The optical model for nucleons and the 
shell model are closely connected concepts: Both involve 
a single-particle central potential, see Section llV.DI 

The average time r for decay of the CN or, cquiv- 
alcntly, the average width F ~ h/r oi the CN reso- 
nances, can be estim ated using an argument due to Weis- 
skopf ((Blatt and Weisskopf . .19521 . For bound levels with 
constant spacing d, the time-dependent wave function (a 
linear superposition of the eigenfunctions) is, aside from 
an overall phase factor, periodic with period d/{2-Kh). 
The wave function reappears regularly at time intervals 
27:%/ d at the opening of any channel a where it escapes 



with probability Ta- For the time Ta for escape into chan- 
nel a this gives = 2T:h/{dTa), the partial width for 
decay into channel a is Fq = Ti/ra = {d/2n)Ta, and the 
total width F is 



(11) 



Although the derivation of Eq. ([TT|) is based on equally 
spaced levels, it is also used for CN resonances, with d 
th e average resonanc e spac ing. A precursor to Eq. (fTTj) 
in lBohr and Wheeleil (|l939t) . F = M/{2n), was based on 
the assumption of strong absorption Ta = 1 in all A open 
channels. 

Eq. ([5]) predicts symmetry of the CN cross section 
about 90 degrees cm. while the Hauser-Feshbach for- 
mula with Tc = 1 in all channels predicts that the energy 
distribution of CN decay products is proportional to the 
density of states in the final nucleus. Early tests of the 
Bohr assumption were focused on these predictions. Ex- 
amples are shown in Figs. [2] and [3] 



Rh""(p.n)Pcl' 
f,-8 MeV 5 = 80" 
r = 780 keV 




4 5 
(McV) 



FIG. 2 Spectrum of neutrons ("evaporation spectrum") emit- 
ted from the CN ""Pd in the (p,n) reaction on ^"^Rh at 
an angle of 80 degrees versus neutron energy En (semilog- 
arithmic plot). With all transmission coefficients put equal 
to unity ("black box" model for the CN), the cross section 
oc exp[— i5„/r] mirrors the level density in the residual nu- 
cleus and per mits the determination of the nuclear tempera- 
ture T. From (|Holbrow and Barschal] . llQeaD . 
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150° 180° 

FIG. 3 Angular distribution of the compound-elastic cross 
section for the ''''Si(p,p) reaction at a bombarding en- 
ergy Ep = 9.8 MeV. The data (dots), taken in the 
regime F >> d, show symmetry of the CN cross sec- 
tion about 90 degrees cm. The solid lines (a) and (b) 
are predictions of the Hauser-Feshbach formula ^ with- 
out an d with an elastic enhancement f actor two, respectively. 
From (|Kretschmer and Wangleij . Il978l 'l . 



D. Ericson Fluctuations 

In the early days of CN theory it was widely held that 
in the "continuum region" of strongly overlapping res- 
onances, F ^ d, t he cross section would be a smooth 
function of energy (|Blatt and Weisskop j . [l952l) . The nu- 
merous overlapping resonances contributing randomly at 
each energy were thought to yield a scatte ring ampli- 
tude that varies slowly with energy. lEricson (1960) re- 
alized that this is not the case. He predicted strong and 
random fluctuations of CN cross sections that would be 
correlated over an energy interval of length F, the aver- 
age width of the CN resonances as defined in Eq. ([TT|). 
Eric son refined his earlier conjecture in a seminal arti- 
cle (lEricsonl 1963 ). Toge ther with theoret i cal wo rk by 
iBrink and StephenI ^9Mj and iBrink et al] (|l964f) . this 
paper became the basis for a large number of experi- 
mental investigations. The almost simultaneous advent 
of electrostatic Tandem van-de-Graaff accelerators that 
produced ion beams with sufhciently small energy spread 
and with the beam energy required to produce compound 
nuclei in the Ericson regime, led to intense experimental 
and theoretical activity and confirmed Ericson's conjec- 
ture. A recent example is shown in Fig. 1141 below. "Eric- 
son fluctuations" , as the phenomenon came to be known. 



have since been found and investigated in many areas 
of physics and constitute one of the most characteristic 
features of chaotic scattering in the regime of strongly 
overlapping resonances. 

To estimate the magnitude of cross-section fluctua- 
tions, Ericson argued that for F ^ all resonances have 
approximately the same width F. To justify that state- 
ment he used the fact that for each resonance, the total 
width can be written as the sum over fluc t uatin g con- 
tributions from the open channels ( Ericsonl . Il963() . The 
fluctuations of the sum are, thus, inversely proportional 
to A, the number of open channels. From the Weisskopf 
estimate (fTTjl we see that F ^ d is possible only for 
A ^ 1 . Thus the fluctuations of resonance widths should 
become negligible for F ^ c? (see, however, Scction lV.Fp . 

We simplify the presentation by excluding the case of 
strictly elastic scattering, and by assuming that Eq. (|4]) 
applies. Taking all resonance widths to be equal to F, we 
write the S'-matrix in the form 



SabiE) = -i 



^ E- E,, 



ia^b). (12) 



The parameters are the partial width amplitudes for 
the decay of resonance fj, into channel a. The fluctua- 
tions of the resonances are due to the stochastic nature 
of the resonance parameters and jafj.- The resonance 
energies E^ and the complex partial width amplitudes 
Ifia = laij. are assumed to be uncorrclated random vari- 
ables; the 7S are complex Gaussian random variables with 
mean values zero; pairs of 7s carrying different indices are 
assumed to be uncorrclated; we write (|7a^P) = 27rt;^ 
where the angular bracket now stands for the ensemble 
average. We refer to Eq. (|12p and to these statistical 
assumptions jointly as to the Ericson model. 

According to Eq. Sab{E) is, at fixed energy E 

and for F 3> rf, a sum over many terms, each term 
containing the product of two Gaussian-distributed ran- 
dom variables, all terms being statistically independent 
of eac h other. We use the central limit theorem to con- 
clude ( Brink and StephenI . Il963l ) that Sab{E) is a Gaus- 
sian random process (i.e., the generalization of a ran- 
dom variable to a random function of some parameter, 
here the energy E). A Gaussian distribution is com- 
pletely defined by its first and second moments. Hence we 
need to determine only {Sab{E)) and {Sab{E)S*j^{E + e)) 
where e denotes the difference between the energy argu- 
ments of S*j^ and of Sab- (We use that quite generally 
(SabScd) = {Sab){Scd), see Section EVl). Then all higher 
moments and correlation functions of S are known. 

Rather than energy averages we actually calculate the 
ensemble averages of (SabiE)) and (SabiE) S*^{E + e)) . 
This is done by averaging over the distribution of the 7^^ 
and of the i?^. Energy and ensemble averages give identi- 
cal results but ensemble averaging seems physically more 
transparent. It is obvious that (SabiE)) = 0, in keeping 
with Eq. (HI). In calculating (SabiE)S*j^iE + e)), we first 
carry out the ensemble average over the 7s. In the re- 
maining summation over E^ we use that for an arbitrary 
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function / we have /(^^m) = / ^E' f{E') 5{E' - 
Ef,). We also use the definition p{E) = (X)^ S{E - E^,)) 
of the average level density p{E) of the resonances, and 
we assume that p{E) is constant over an interval of length 
r. This yields 



{Sab{E)S*,l{E + e)) = {5ac5M+5ad5bc 

x^^vlvlp{E) 



1 



T + ie 



(13) 



We compare Eq. ([T5|) for e = with the Hauser-Feshbach 
formula ([7]). Complete agreement is obtained when 
wc use the Weisskopf estimate Eq. (fTTj) . the identity 
d = l/p{E), and write for the transmission coefficients 
Ta = 47r^v^/c»(_E). (In anticipation we mention that that 
result agrees with Eqs. and (gS]) if Ta < 1 for all 
channels). This shows that the Ericson model yields the 
Bohr assumption. It docs so, however, only upon aver- 
aging. It is obvious that \Sab{E)\'^ as given by Eq. 
docs not factorizc as it stands. Calculating the energy av- 
erage (rather than the ensemble average) with the help 
of a Lorentzian weight function of width /, wc find that 
the averaging interval / must be large compared to F for 
Eq. ([13]) to hold. Thus, independence of formation and 
decay of the CN hold only for cross sections averaged 
over an interval of length / ^ F. CN cross sections mea- 
sured with particle beams of sufficient energy resolution 
are expected to deviate from the Hauser-Feshbach for- 
mula. Since Eq. ([5]) involves an average, we expect the 
same statement to apply to the symmetry of CN cross 
sections about 90 degrees cm. In contrast to the pre- 
sentation chosen in Section III. C[ it had not always been 
clear prior to Ericson's work that the Bohr assumption 
holds only upon averaging. 

The magnitude of cross-section fluctuations is es- 
timated by calculating the normalized autocorrelation 
function of \Sab{E)\'^. We use the Gaussian distribution 
of the S'-matrix elements and Eq. (fT3|) and obtain 



1 



{\Sa,iE)\'\Sa,iE + e)\')^{{\Sa,{E)\')r 

i{\Sab{E)\^)y 



1 



(£/F)2 ■ 
(14) 

Eq. ([H)) shows that the fluctuations have the same size as 
the average cross section, and a correlation width given 
by F, the average width of the CN resonances. These 
results suggested that F could be measured directly us- 
ing Eq. (HH). The Lorentzian on the right-hand side of 
Eq. (|T4| signals that the CN decays exponentially in time. 
The hfctime is h/T. 

The Ericson model can be e xtended to include the case 
of strictly ela stic scattering ( Brink and Stephenl . Il963l : 
lEricsonl (l963l ). Moreover, the arguments can straight- 
forwardly be extended to the fluctuations both in en- 
ergy and angle of the ac tual CN cross section as given by 
Eq. ^ (|Ericsonl[l963l) . 

In summary, in the regime of strongly overlapping reso- 
nances (F ^ d) the Ericson model leads to the following 
conclusions: (i) The elements Sab{E) of the scattering 
matrix arc Gaussian random processes, (ii) The Bohr 



assumption (i.e., the Hauser-Feshbach formula) and the 
symmetry of the CN cross section about 90 degrees cm. 
hold only for the energy-averaged cross section but not 
for cross sections measured with high energy resolution. 
The averaging interval has to be larger than F. (iii) CN 
cross sections (including the angular distributions) fluc- 
tuate. The fluctuations have the same magnitude as the 
average cross section. The correlation length of the fluc- 
tuations is F. This fact can be used to measure F. 



E. Direct Reactions 

Experiments in the 1950s performed with poor energy 
resolution showed that CN cross sections are not always 
symmetric about 90 degrees cm. These results were con- 
firmed with better resolution. For the case of elastic scat- 
tering, an example is shown in Fig. 21 Although the data 
show the cross section versus lab angle, it is clear that 
the asymmetry persists also in the cm. system. Similar 
results were also oblained for inelastic processes. Such 
deviations were attributed to a failure of Eq. i.e., to 
inelastic scattering processes without intermediate for- 
mation of the CN. Theoretical efforts to develop a the- 
ory of such "direct reactions" dominated nuclear reaction 
theory for some years and led to explicit expressions for 
the non-diagonal parts of (S) . That raised the question 
how the Hauser-Feshbach formula has to be modified in 
the presence of direct reactions, i.e., when (S) is not di- 
agonal. 




FIG. 4 Elastic cross section for protons at 5 MeV scattered 
on ^^Mg versus scattering angle. The data points with error 
bars are compared with predictions of the Hauser-Feshbach 
formula (NC) without (solid line) and with elastic enhance- 
ment factor (dashed l i ne) a nd of the optical model (ID). 
From (|Gallmann et d.l . ll966l ). 
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After contributions due to lMoldaueil (|l96llll963lll964l ) 
and lSatchleil (I1963D . a definitiv e answer was propo sed by 
Kawai, Kerman, and McVoy ( Kawai et al\ . Il973l ). We 
briefly describe the approach for F ^ d. The starting 
point is a decomposition of the scattering matrix of the 
form 



Sab = {Sab) - 



9^a9^b 
E-£a 



(15) 



The first term on the right-hand side is the average 
S-matrix, and the last term (which is identical to 5*^) 
represents the resonance contributions, with £^ denot- 
ing the complex resonance energies. By definition, that 
term averages to zero. The decomposition (jlSp im- 
plies ( Kawai et al\ . Il973[ ) correlations between g^a and 
gfj^f, for a 7^ 6 if {Sab) 7^ for a ^ h. In that respect the 
approach differs fundamentally from the Ericson model. 
Subsequent developments are similar, however, to those 
sketched in Section flLDI In particular, the total widths 
of the resonances arc effectively assumed to be constant 
(independent of /x). The average CN cross section is ex- 
panded in powers of d/T. The term of leading order is 



\Slb?) 



1 



TraccP 



PaaPbb + PabPal 



Here 



Pab = 5ab - '^{Sac){{Shc)* 



(16) 



(17) 



is Satchler's transmission matrix (ISatchleillTgel . It gen- 
eralizes the transmission coefficients defined in Eq. ([8]) to 
the case where {S) is not diagonal. It meas ures the uni - 
tarity deficit of {S). According to (Kaw ai et al\ . 1197,'^ . 
Eq. ([T5|) replaces the Hauser-Feshbach formula (jH]) when 
direct reactions are present. We observe that if {S) is 
diagonal, Eq. ([TB| reduces to the Hauser-Feshbach for- 
mula, including the elastic enhancement factor. 



F. Limitations of the Compound-Nucleus Picture 

Nuclear reaction cross sections begin to show devia- 
tions from the CN picture some 20 MeV or so above the 
first particle threshold. This failure of the CN model is 
attributed to "preequilibrium" or "precompound" pro- 
cesses. The deviations occur because basic assumptions 
on characteristic time scales of the reaction made in CN 
theory do not apply any longer. In Eq. ([3]) such an as- 
sumption is implicitly made. The average taken to cal- 
culate {S) must obviously extend over an energy interval 
containing many resonances. By the uncertainty rela- 
tion, a large energy interval relates to a short time in- 
terval for the duration of the reaction, see Section IV. CI 
Thus, {S) describes the fast part of the CN reaction. We 
estimate the minimum duration time of the reaction by 
the time of passage vp/R of a nucleon with Fermi veloc- 
ity vp through a nucleus with radius R. In contrast, the 



fluctuating part of S with its rapid energy dependence 
describes the slow part of the CN reaction: Decay of the 
CN with an average lifetime fi/T, with F given by the 
Weisskopf estimate pT|) . For the decomposition ^ to 
be meaningful, these two time scales must be well sepa- 
rated, R/vp < h/T. 

An even more stringent constraint on time scales is 
hidden in the assumption that the CN resonances obey 
GOE statistics. As shown in part I, that assumption is 
experimentally validated in the regime of isolated reso- 
nances. But does it also hold for higher excitation en- 
ergies where F 3> d? The assumption implies that the 
partial width amplitudes of the CN resonances (defined 
here as the eigcnfunctions of the GOE Hamiltonian in- 
troduced in Eq. ([36)) below) have a Gaussian distribu- 
tion in each channel. In other words, the couplings of 
all resonances to a given channel are, within statistics, 
equally strong, and there does not exist a preferred state 
or configuration, or a group of such states. (The coupling 
strengths may differ, of course, for different channels). 
This is the formal expression of Bohr's assumption that 
the CN "equilibrates" . To display the time scale hidden 
in that assumption, we consider a nucleon-induced reac- 
tion. The nucleon shares its energy with the nucleons in 
the target in a series of two-body collisions. Configura- 
tions of ever greater complexity are created. It takes sev- 
eral or perhaps even many such collisions until the energy 
of the incident particle is shared among many nucleons 
and the situation described by assuming GOE statistics 
for the resonances, is attained. The time elapsed between 
the first collision and the attainment of equilibrium is 
the "equilibration time" Tcq. The GOE description of 
CN scattering holds if decay of the CN sets in after the 
nucleus is equilibrated, i.e., if Tcq is smaller than the av- 
erage decay time h/T. Obviously, wc have vp/R < r^q 
and the condition r^q < h/T for applicability of RMT is 
usually more stringent than the condition vp/R <C h/T 
deduced from the decomposition ([3]). We observe that 
h/Teq can be interpreted as the spreading width of the 
simple configuration created in the first collision between 
the incident particle and the target nucleus. The condi- 
tion Tcq < h/T then requires that spreading width to be 
large compared to the CN decay width F. 

Precompound decay sets in when these conditions are 
violated. Simple estimates show that Teq <C h/T for 
isolated resonances where T <^ d. But the complexity 
of CN resonances and, therefore, Teq increases with E. 
More importantly, F increases strongly with excitation 
energy, and h/T decreases correspondingly. Thus, while 
Tcq <ti h/T for isolated resonances, h/T becomes compa- 
rable with and eventually smaller than Tcq as E increases: 
Particle decay is possible before the equilibrated CN is 
reached in its full complexity. The associated reaction 
times cover the entire range from vp/R to h/T, and both 
the GOE description of resonances and the decomposi- 
tion ([3]) no longer apply. By the same token, precom- 
pound decay does not have the characteristic features of 
CN reactions: The average energy of the emitted par- 
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tides is larger than for CN decay; the emitted particles 
"remember" the incident channel(s) and emission is pref- 
erentially in the forward direction. Pre-equilibrium re- 
actions obviously call for a different approach, although 
the large number of configurations involved cannot be 
handled without statistical assumptions. Phenomeno - 
logical models designed for these reactions ( Blannl . [19751 ) 
were later followed by theories aiming at a quantum- 
statistical d escription. The latter were c ompared and 
analyzed by iKoning and AkkermansI ( 1990f ) . The theory 
of CN reactions reviewed in this paper is a closed theory 
based on the concept of quantum chaos as embodied in 
RMT and on strong empirical evidence reviewed in part 
I of this review. In its essential aspects CN scattering is 
universal and occurs likewise in the transmission of waves 
through disordered media. The theory of precompound 
reactions is specific to nuclei. It requires additional as- 
sumptions that go beyond chaos and RMT that cannot 
be tested directly. It is not reviewed here. 

Precompound processes occur in reactions induced by 
light particles (mass A < 6 or 8 or so) impinging on a 
target of mass A > 20 or so. A very different situation is 
that of reactions between "heavy ions", i.e., of two nuclei 
of mass A > 50 or 100 each. The energy of relative mo- 
tion is usually given per nucleon. An energy of 5 MeV 
per nucleon may then easily amount to a total kinetic 
energy of several 100 MeV. In the case of a grazing colli- 
sion, the reaction transports energy and angular momen- 
tum of relative motion into intrinsic excitations of either 
fragment. This is accompanied by the transfer of nucle- 
ons between both reaction partners. Excitation energies 
of several 10 MeV are easily reached in either fragment. 
But the two fragments basically keep their identity; the 
CN corresponding to complete fusion of both fragments, 
is typically not reached in a grazing collision. The theory 
of such processes uses concepts like friction and dissipa- 
tion developed in the theory of non-equilibrium processes 
and methods of quantum statistical mechanics. Natu- 
rally these arc inspired by RMT and chaos but, to the 
best of our knowledge, have never been strictly derived 
from such a basis. This is why in the present review we 
will not cover that area of the statistical theory of nuclear 
reactions either. In central collisions between two heavy 
ions, fusion is possible even though the density of CN res- 
onances is limited: Quasibound states of the shell model 
with too lar ge single-particle wid ths do not qualify as CN 
resonances ( Weidenmiilleil . [l963 ). Such fusion processes 
are mainly investigated via the gamma rays emitted in 
the decay of the highly excited CN; the resulting data do 
yield statistically relevant information on excited states 
of the CN; that information is discussed in Part I of this 
review. The literature on the subject is vast because 
the available energies have been much increased over the 
past 20 years, and because highly segmented gamma-ray 
arrays have become avail able. We confine ourselves to 
citing two early rev i ews ( Norenberg and Weidenmiilleil . 
Il980i: IWeidenmiilled . [l980h . 



III. RANDOM-MATRIX APPROACH TO QUANTUM 
CHAOTIC SCATTERING 

The statistical models reviewed in Section |TT1 although 
inspired by RMT, are not derived from a random-matrix 
description of CN resonances. The Bohr assumption 
is intuitively appealing. But it is not clear in which 
range of T /d it applies and with what accuracy. The 
Ericson model leads to interesting predictions that agree 
with experiment. However, completing the model for- 
mulated in Eq. ([T^ by an equation for the elastic case, 
i.e., for Saa, one finds that the resulting ^-matrix vio- 
lates unitarity. The same statement applies to the in- 
clusio n of direct reactions in the Hauser-Feshbach for- 
mula ( Kerman and Sevgenl Il976l ) . As for the Weisskopf 
estimate (jlip . the physical significance of the parameter 
F is not clear. A width parameter appears in the context 
of the ^-matrix autocorrelation function, see Eq. (|14p . 
We may also define F as the mean distance of the poles 
of S from the real axis. Are the two definitions identical? 
If not, which of the two (if any) is given by the Weis- 
skopf estimate? What is the accuracy of the estimate? 
What are the correction terms of next order? A theoret- 
ical approach based upon a random-matrix description 
of CN resonances that yields ^-matrix distributions and 
S'-matrix correlation functions within controlled approx- 
imations, is clearly called for. That approach is reviewed 
in the present and in the following two Sections. In the 
present Section we formulate the approach and, at the 
end of the Section, give a brief historical survey. The 
approach uses the average 5-matrix as input. Properties 
of and phcnomenological models for {S) arc reviewed in 
Section IIVI Results of the random-matrix approach are 
given in Section [V] The approach reviewed in this and 
the following Sections is not confined to CN reactions 
but applies in general to a random-matrix description of 
quantum chaotic scattering. 



A. Resonance Reactions 

A statistical theory based on RMT can be formulated 
only on the basis of a theory of resonance reactions. Only 
when the S'-matrix is written explicitly in terms of res- 
onance contributions can we implement the statistical 
properties of those resonances. We describe an approach 
that is based upon the coupling of N quasibound states 
to a number of channels and yields the S'-matrix directly 
in terms of the Hamiltonian governing those quasibound 
states. We show why this approach is very well suited 
for calculating GOE averages. 



1. Single Resonance 

Wc begin with the simplest case, a single resonance 
without any background scattering. The S-matrix has 
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Breit-Wigner form, 



and ( Mahaiix and Weidenmiilleil . llOGOfO 



Sab{E) — Sab — i 



lalb 



E-Eo + (z/2)r 



(18) 



The partial width amphtudcs ^a, lb, ■ ■ ■ give the proba- 
bihty amplitudes for decay of the resonance into channels 
a, b, .... Factorization of the numerator in the last term 
in Eq. (|18p is implied by quantum mechanics if the reso- 
nance is caused by a single quasibound state, see below. 
The matrix S{E) is obviously symmetric. We impose the 
unitarity condition ([ij for all energies E and find that all 
the partial width amplitudes must be real and that 



(19) 



The sum of the eigenphase shifts J2a of S increases 
by TT as i? increases from a value far below to a value 
far above Eg. To see this, we note that the form ([2]) 
for S implies detS' = cxp{2iJ2a^a)- To calculate det S* 
we use Eq. (jlSp and find in matrix notation det S = 
det(l - iii^ l{E - Eo + {i/2)r)) where T denotes the 
transpose. The last determinant is easily worked out and 
yields [E - Eq ~ {i/2)V)/{E - Eq + (V2)r). That term 
has magnitude one; the phase increases by 27r as E passes 
through the resonance. The resulting increase by tt of the 
sum of the eigenphase shifts of S over the resonance does 
not imply that one of the eigenphase shifts grows by tt 
while all others remain unchanged. On t he contrary, the 
non-crossing theorem of the eigenphases ( Weidenmiilleil . 
Il967l ) referred to above implies that generically all eigen- 
phases increase as E passes through the resonance, the 
average increase of every eigenphase being tt/A where A 
is the number of open channels. 

Eq. ([T5)) can be derived by considering a single qua- 
sibound state |(/)) with energy e which interacts with a 
number A of continuum states \xa{E)) (the range of the 
channel index a is A). The continuum wave functions are 
orthonormal, {xa{E)\xb{E')) = dabS{E - E') for all a,b 
and orthogonal to the quasibound state, {xa{E)\(l)) ~ 
which in turn is normalized, {4i\4i) = 1. In keeping with 
the absence of any background scattering in Eq. we 
neglect the dynamical coupling of the continuum states 
with each other and assume that the states \xa{E)) de- 
scribe free relative motion so that the scattering phase 
shift for each \xa{E)) vanishes. We focus attention 
on the coupling between the continuum states and the 
quasibound state mediated by the real matrix elements 
Wa{E). That coupling causes |0) to become a resonance. 
The Hamiltonian is 

H^e\4,){cp\+Y, I dE E\xa{E)){xa{E)\ 

+ J2 I dEWa{Em{Xa{E)\ + \Xa{E)){cp\) . 



(20) 

To determine t he s cat tering e igenst ates of 

H we write ( dPirad . Il958h . (iFanol . Il96lh . 



AE' ai^\E,E')\xc{E'))+b^^\EM 



(21) 

Here j^*) has an incoming wave in channel a only. We 
insert Eq. (|2T|) into the Schrodinger equation H\'^) = 
El"^) and multiply the result from the left with ((/)| and 
{XaiE)\. That gives a set of coupled linear equations for 
the coefficients a{E,E') and b{E). These are solved by 
imposing the boundary condition for |4'+(i?)). We define 

F{E) = Y.jdE' -^i^ =^-Y,^^Wm ■ (22) 

Here E'^ carries an infinitesimal positive increment. That 
corresponds to the boundary condition on 1^'"'"). The 
shift function A is defined in terms of a principal-value 
integral. We obtain 



b^''\E) 
ai^\E,E') 



Wa{E) 



E-e-F{E) ' 
6acS{E - E') 

1 Wa{E)Wc{E') 



E+ -E' {E-e- F{E)) 



(23) 



In the second of these equations we have again used the 
boundary condition: After insertion into Eq. (|2ip the 
delta function yields an incoming wave in channel a while 
the denominator 1/{E~^ ~ E') in the last term makes 
sure that all other contributions to |\I'+(£')) produce only 
outgoing waves at infinity. The amplitude of these waves 
can be easily worked out and is 



Sab{E) 



6ab — 2iTT 



WaiE)WbiE) 

E-e-A + {i/2)r 



(24) 



Here F = 2'k^^W^{E). We have used the super- 
script "res" on S to indicate the origin of Eq. ([M)) from 
a dynamical calculation employing a quasibound state. 
Eq. ([24]) agrees with Eq. ^ if we put 7a = y/2^Wa{E) 
and i?o = £ + A. The factorization of the numerator in 
the second of Eqs. (|23p causes the factorization of the res- 
onance numerator in Eq. ([M)) and holds for all resonances 
caused by a single quasibound state. 

To interpret Eqs. (gH) and ([Ml), we note that F{E) is 
the sum of the resolvents for free motion in the channels: 
The first factor Wc{E') is the amplitude for decay of the 
quasibound state into channel c at energy E' , the denom- 
inator l/(i?^ — E') is the propagator for channel c, and 
the second factor Wc{E') is the amplitude for returning 
to the resonance. We have to sum over all energies E' . 
The occurrence of F{E) = A — (z/2)F in the denominator 
of Eq. signals repeated decay of and return to the 
quasibound state. By that mechanism the quasibound 
state turns into a resonance. The real part A of F{E) 
corresponds to off-shell processes and gives the shift of 
the resonance energy e. While is due to the open 
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channels only, A receives contributions from both, open 
and closed channels. 

We turn to the general case where the background scat- 
tering is not negligible and leads to a background scat- 
tering matrix The S-matrix for a single resonance 
takes the form 



(0) 
ab 



lalb 



E-En 



(V2)r 



(25) 



The matrix S^^^ is unitary. This follows from the uni- 
tarity of S at energies far from the resonance energy 
Eq. Unitarity of S at all energies E implies that the 
complex partial width amplitudes obey J2c It'cP ~ ^■ 
The matrix S^^^ is also symmetric. This is implied 
by time-reversal invariance. Therefore, S^°^ possesses 
an eigenvalue decomposition of the form ([2]), S^J = 

X)c exp(2i(5c°^)C'j|!^ We use that relation and ma- 
trix notation to write Eq. (|25|) in the form 

Sab{E) = exp(i(5(o))5^"^(S)exp(i(5(o))(OW)^ 



ab 

(26) 

where the matrix S^^'^ in the cigcnchannel representation 
of has exactly the form of Eq. ((T8)) , 



SlTiE) 



Ocd 



Icld 



E-En 



(z/2)r 



(27) 



For the unitarity condition ([T]) to hold for all energies E, 

the partial width amplitudes 7c = exp(— i(5c°'') J2a ^ca^a 
must be real and must obey Eq. (|19p . It is normally as- 
sumed that 5*^"^ and the 7s are independent of energy. 
This approximation assumes that angular-momentum 
and Coulomb penetration factors and all relevant matrix 
elements depend smoothly on energy and is typically jus- 
tified over the width of a resonance although counterex- 
amples exist, see Section IVII. A. 21 The assumption does 
not apply for reactions at threshold energies which are 
not treated here. Analogous assumptions are made also 
when we later deal with many resonances; these are not 
mentioned explicitly again. It is easy to check that the 
sum of the eigenphases of S increases by it over the width 
of the resonance. 

As in the case of Eq. p8)) . it is possible to derive 
Eq. (|26p from a dynamical model. The Hamiltonian dif- 
fers from the one in Eq. (PH]) in two respects. First, the 
continuum states \xa{E)) are true scattering states with 
nonzero phase shifts. Second, the continuum states car- 
rying different channel indices interact with each other. 
Taking into account these phase shifts and continuum- 
continuum interactions first and neglecting the quasi- 
bound state yields the background matrix S''^\E). The 
actual calculation of that matrix involves a coupled- 
channels problem and may be rather involved. We show 
in Section IIVI how that problem is overcome. In the 
eigenchannel representation Eq. ^ of S*'"', the radial 
part of relative motion of the continuum wave functions 
exp{—iSc^'^)rtciE)) is real. (These functions are defined 



below Eq. ([2]). For neutral particles the radial part in 
any channel c is essentially given by a spherical Bessel 
function ji^{kcrc + Sc))- Therefore, the matrix elements 
Wa{E) coupling the quasibound state to the continuum 
wave functions are also real in the eigenchannel repre- 
sentation. A calculation analogous to the one leading 
to Eq. ^ then yields Eq. (gH) with 7^ = V2^Wa{E) 
and the resonance energy Eq given by e -I- A as be- 
for e. We do not give that calculat ion here but refer 
to ( Mahaux and Weidenmtillel . ll969D . In summary, the 
S'-matrix is given by Eq. (^5)) and S'^^^{E) has the form of 
Eq. except that the real matrix elements Wc{E) now 
refer to the eigenchannel representation of S*^"-* . Eq. ((26)) 
is very convenient as it clearly separates the effects of 
the continuum-continuum interaction and the effect of 
resonance scattering. 



2. N Resonances 

Except for the factorization of the numerator in the 
last term in Eq. (fTS)) . it is straightforward to guess the 
form of the S'-matrix for a single resonance. In the case 
of several or many overlapping resonances it is difficult to 
guess a form for S that fulfills the unitarity conditions ^ 
for all values of the resonance parameters (partial width 
amplitudes, resonance energies and total widths). There- 
fore, it is useful to derive that form from a dynamical 
model. We describe the model but skip the derivation 
because it is quite similar to the derivation leading to 

Eq. dm). 

We consider a set of A orthonormal continuum states 
\Xa{E)) and a set of N orthonormal quasibound states 
\(t>f_i). The continuum states interact with each other. 
As in Section IIII.A.ll that interaction gives rise to 
a smooth unitary and symmetric background scat- 
tering matrix S'^^'' for which wc can write S"^^' = 

^^Oic^ exp(2i(5^"')C'[°^ The quasibound states span an 
A^-dimensional Hilbcrt space. The unperturbed energies 
of and the interactions amongst the quasibound states are 
not specified in detail but arc represented jointly by the 
real and symmetric iV-dimensional matrix H, the pro- 
jection of the total Hamiltonian onto the TV-dimensional 
Hilbert space of quasibound states. In that space, we 
use an arbitrary basis of orthonormal states and write 
the matrix elements of H as iJ^^ with fi,i/ — I, . . . ^ N. 
As before we define the interaction between the qua- 
sibound states and the continuum states in the eigen- 
channel representation of the matrix S*'"' with scattering 
eigenfunctions eyi'p{~iS^^)na{E). Then the matrix ele- 
ments wjl^ (E) = Wa^ (E) of that interaction are real. 
The upper index zero indicates that we use the eigenchan- 
nel representation of S^^^ . The total S'-matrix again has 
the form of Eq. ([25)1 but the resonance part now differs 
and is given by 

S;r(£;) = Sab~2^7^Y,W^°JiE)iD-'),^w!^'^\E) . (28) 
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Here D{E) is a matrix in the space of quasibound states 
and given by 



(29) 



with 



F^,{E) ^Y./^^ 



, w^J'}{e')wS,\e') 



E+ -E' 

A<~°}-^nJ2w|,'J{E)W^'J{E) . (30) 



It is easy to check that S''''^{E) in Eq. (^5]) is unitary. Ac- 
cording to Eq. (j26p this statement imphes unitarity of the 
total S'-matrix. Also 5'''°^ is obviously symmetric, and so 
is, therefore, S{E). We note that the clear separation 
of the influence of the continuum-continuum interaction 
and that of the quasibound states expressed by Eq. ([26)) 
is not restricted to the case of a single resonance but holds 
in general. We also observe that Eqs. (^5]) to (PH)) apply 
for any strength of the coupling between resonances and 
channels. The equations are not restricted to the regime 
of weakly overlapping resonances and will serve as a basis 
for a statistical theory that applies for all values oiV/d. 

As for the case of a single resonance, F{E) given in 
Eq. ([50)1 describes virtual (real) decay of and return to 
the quasibound states via all channels by its real (imag- 
inary) parts, respectively. This gives rise to the real 
shift matrix A^*^) and the real width matrix with ele- 
ments 27r J2c {E)W^°^ (E) . Qualitatively speaking, 
the resonances are isolated (they overlap) if the average 
resonance spacing is large (small) compared to typical el- 
ements of the matrix F{E). That statement is quantified 
in Sections ITVl and IVJ] below. 

To define the analogue of the total width for a single 
resonance, we have to write the matrix 5*'°^ as a sum over 
poles in the complex energy plane. That is straightfor- 
ward only for isolated resonances. The Hermitean matrix 
H can be diagonalized by an orthogonal transformation 
O. We denote the eigenvalues of H by E^ and write for 

the transformed matrix elements W^l? = ^^OnijWl)^^. 
With a corresponding notation for F the transformed 
matrix D takes the form (£'+ — E^)5^v — F^^. The N 
eigenvalues E^ are coupled to each other by the non- 
diagonal elements of F. For isolated resonances, that 
coupling can be neglected, and only the diagonal ele- 
ments of F are retained. We omit A^°) in Eq. define 

r,. = 27rE.(W^ic^)' 
the form 



and obtain for isolated resonances 



SaTiE) — Sab 

-2^nY,W^°HE){E - E, + i^ /2)r .y'wj^f 



(31) 



This equation explicitly displays N resonances. For the 
resonance labeled /i, the partial width amplitude Wa^ 



gives the probability amplitude for decay into channel a, 
and the total width F^ is the sum of the partial widths 

|W^a?P- The scattering matrix (PT|) is unitary if we ne- 
glect the overlap of different resonances. The more gen- 
eral Eq. (pS)) obviously describes an S'-matrix with N 
resonances, too, but applies also outside the regime of 
isolated resonances. 

The S'-matrix in Eq. ([25]) possesses N poles in the 
complex energy plane. All poles of S"'^^ are located be- 
low the real energy axis. To see that we observe that 
the positions of the poles are given by the eigenvalues of 
Hfiu + F^i, . Let z be one such eigenvalue. It obeys 



(32) 



We multiply Eq. ([5^ by the complex conjugate eigen- 
function (5*^)*, sum over and take the imaginary 
part of the resulting equation. That yields the "Bell- 
Steinberger relation" 

-'^EiE^^A^^A^i'^Ei^^'i'^w . (33) 



This shows that indeed 3[z] < unless 'Yl,^l,^c^^'^ ti = 
7c = for all c. In the latter case we deal with a bound 
state embedded in the continuum with vanishing partial 
width amplitudes 7c. Then z is real and 5''°^ does not 
have a pole at E = z. 

The scattering matrix 3"°^ in Eq. ([^5)) can be rewrit- 
ten in terms of the i^-matrix, a matrix in chan- 
nel space. The ensuing relations and (|55|) are 
quite general and no t restricted to isolated reso n ances . 
They may be used ( Mahaux and Weidenmiillerl . Il969l ) 
to esta blish the relation between t he i?-matrix ap- 
proach ( Wiener and Eisenbudl Il947t) and the present 
framework. Using the orthogonal transformation that 
leads to Eq. ([3T|) we define 



Kab{E)=7:Y, 



E - E„ 



and find after a straightforward calculation 



SabiE) 



1 -iK 
1 + iK 



(34) 



(35) 



ab 



B. Stochastic Scattering Matrix 



The result of Section IIII.A.2I puts us into the position 
to introduce the stochastic description of the resonances 
in terms of the GOE. In a dynamical treatment of the res- 
onances based on the shell model we would express the 
matrix H in Eq. ([29)) in terms of the single-particle ener- 
gies and the matrix elements of the residual interaction. 
We replace such a dynamical treatment by a stochastic 
one and consider iJ as a member of the GOE. In other 
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words, in Eq. ([29|) we replace the actual Hamiltonian ma- 
trix H by the ensemble i?^*^^. That replacement gen- 
erates an ensemble of scattering matrices S each mem- 
ber of which is obtained by drawing H from the GOE. 
For the sake of completeness we recall the definition and 
some properties of the GOE and refer to Part I for fur- 
ther details. The independent elements H^^^ of the real 
and symmetric A^-dimensional matrix i/GOE uncor- 
related Gaussian-distributed random variables with zero 
mean values and a second moment given by 



^^GOE^GOE' 



"per 



\2 



(36) 



Here A is a parameter. The average spectrum of H^'^^ 
has semicircular shape. The radius of the semicircle is 
given by 2A. Near the center of the spectrum the av- 
erage level spacing d is given by d = nX/N. The GOE 
is invariant under orthogonal transformations of the un- 
derlying Hilbert space. The energy argument E of S 
is taken close to (i.e., a finite number of spacings away 
from) that center. All observables are calculated in the 
limit iV — >■ cxD. As pointed out in the Introduction, we 
expect the replacement of H by iJ^iOE j-^^ valid at 
bombarding energies of up to 10 or 20 MeV where the 
internal equilibration time of the CN is smaller than the 
decay time h/T. 

The statistical approach defined by the replacement 
H — > ijGOE gggjj^g contain a large number of parame- 
ters, i.e., the NA matrix elements Wa^ and the param- 
eter A of the GOE. To appreciate the simplifications due 
to a statistical treatment, it is instructive to visualize the 
effort required in case we would want to determine these 
parameters from a dynamical model such as the shell 
model. Using a shell-model basis of scattering states 
we would first have to calculate the background matrix 
5'(o)^ We would then have to calculate the matrix ele- 
ments wj^} = exp{—iSa){^a\V\(l)f_i,) in the eigenchannel 
representation of S*'"^ in terms of the residual interaction 
V by constructing the quasibound states 0^. For large 
N that effort would be huge. We now show that in the 
statistical approach, all this is much simplified by the 
invariance properties of the GOE. 

The statistical theory uses as input the elements (Sab) 
of the average S'-matrix (S) defined by the decomposi- 
tion ([3]). These elements must be given in terms of some 
dynamical calculation, or some suitable model. Options 
for a calculation of (S) are reviewed in Section IIVI The 
statistical theory aims at predicting moments and corre- 
lation functions (defined as ensemble averages) of (E) 
in terms of (S) in the limit N — > oo. The orthogonal 
invariance of the GOE implies that all such moments 
and correlation functions can depend only on orthogonal 
invariants of the parameters Wa^ . The only such invari- 
ants are the bilinear forms ^fj^Wa^wj^K But the S'- 
matrix is dimensionless and can depend only on dimen- 
sionless combinations of the parameters of the statistical 
approach. These are the real quantities Wa^wj^^ /X. 



Their number is A(A -|- l)/2, equal to the number of ele- 
ments of (S). But the elements of (S) are complex and it 
even seems that the model is overdetermined. We show 
in Section IIVI below that the moments and correlation 
functions of S actually depend, aside from overall phase 
factors, only on the magnitudes |(S'afc)| of the elements 
of the average scattering matrix. Thus, the statistical 
theory is well defined. We also show below that (S) de- 
termines not only the invariants Wa^wj^^ / X but also 

the orthogonal matrix O'"^ and the phase shifts Si^^ ap- 
pearing in Eq. ([26]) . In other words, knowledge of {S) 
suffices for a complete calculation of CN scattering pro- 
cesses. 

The quantities J^n^a^ ^jS'/^ ^'^^ effective pa- 
rameters of the statistical theory. We use this fact 
to simplify further the resonance part of the scatter- 
ing matrix in Eq. (|28p . The real and symmetric ma- 



trix J2fi Wa^wj^^ can be diagonalized by an orthogonal 
transformation O^^ in channel space. We define the new 
real matrix elements 



(0) 
bfi 



which obey 



J2 W^apM^^b = SabNvl 



(37) 



(38) 



To express Sat in terms of the VFa^s we define the unitary 
matrix 



Uab= O(°)cxp(z5(0))OC^ 



ab 



and write 



Here 



Sab{E) = ( C/5CN(i?)C/^ 



(39) 



(40) 



S^^{E) = Sab- 217: J2Wa^iE){D-')^,W,b{E) , (41) 



and 



where 



D^i, {E) — ES^y - H'^Jl 



TTCS _ TtGOE I TP 



(42) 



(43) 



is the non-Hermitean "effective Hamiltonian" that de- 
scribes the dynamics of the resonances. The non- 
Hermitean part of is due to the imaginary part of 
the matrix Ff^y{E) given by 



F^y{E) = J2 J di?' 



„ W ^aiE')WcAE') 
E+ - E' 



(44) 
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That matrix describes the effect of the coupling to the 
channels. 

The matrix S'^^iE) of Eq. (gl]) is the S'-matrix for CN 
scattering in its purest form: Absence of all continuum- 
continuum interactions and of all elastic scattering phase 
shifts, with real coupling matrix elements Wafj, obeying 
Eq. (1221). Thus, S'^^ is the central object of study of 
the statistical theory. The ensemble of matrices S''^^ (E) 
constitutes a matrix-valued random process. (For fixed 
E, every element of S™ is a random variable. Because 
of the dependence on energy, a generalization of the con- 
cept of random variable is called for, and one speaks of 
a random process.) The moments of S'^^ depend on the 
dimensionless parameters 



(45) 



A 



see Eq. (pSl) , where we have taken the average GOE level 
spacing d at the center of the GOE spectrum. The form 
of Eq. (l45l) is reminiscent of Fermi's golden rule. The 
S'-matrix correlation functions depend, in addition, on 
energy differences given in units of d, see Section IV] 

The unitary matrix U in Eq. (l39l) contains all of the 
non-statistical effects that connect S'^^ with the physical 
channels: The transformation to the eigenchanncl repre- 
sentation, the eigenphase shifts of S^^\ and the transfor- 
mation that diagonalizes the bilinear form of the W^^'^s, 
see Eq. (|37p . Statistical assumptions are made only with 
respect to the matrix if^OE ^j^g^^ describes the interac- 
tion among quasibound states. We do not impose any 
statistical requirements on the matrix elements Wp,c{E), 
on , or on U. This is in keeping with the evidence on 
chaos in nuclei reviewed in Part I. It also corresponds to 
chaotic features of quantum dots and microwave billiards. 

In many cases it is possible to simplify the matrix 
Ffi,u{E). As in Eq. (pOl) we decompose Ff^^{E) into its 
real and imaginary parts. Often the coupling matrix 
elements Wfj.a{E) are smooth functions of E (so that 
d [d\nWf^c{E)/dE] <C 1). Then the principal-value inte- 
gral is small, the real shift matrix can be neglected, and 
'(i/2)r^i, where the width matrix F^^ is given 



F. 



by 



T^,^2nJ2w^,ciE)W,,iE) . 



(46) 



The assumption d [dlnWf^dE) /dE] < 1 is not justi- 
fied near the threshold of a channel, a say, where the 
Wfj,aS with fi = 1,...,N depend strongly upon en- 
ergy. A better approximation is obtained when we write 
the principal- value integral giving the real part of 
in th e approximate form ( Mahaux and Weidenmiilleil . 
I1969D J2c W^^cAcWc^ where the A^s are some channel- 
dependent constants. Then F^i, « tt J^c W^/ic(Ac — «)Wci/- 
That form differs from the first approximation by the 
replacement of the channel propagator —i by (Ac — i). 
lAgassi et all ([l9li) have shown that such a replacement 
only modifies the dependence of the transmission coef- 
ficients Tc on the coupling strengths Xc and leaves the 



theory otherwise unchanged. But the coefficients Tc serve 
as phcnomenological input parameters anyway, and the 
coupling strengths Xc do not appear explicitly anywhere 
in the final expressions for the moments or correlation 
fimctions of S*^^. Therefore, we use in the sequel the 
simple approximation « — « X^c ^^c^ci/- With this 
approximation, the effective Hamiltonian in Eq. flSl) 
takes the form 



rrcff 



(47) 



Eqs. (|40|41|42|1 and ^ together with Eq. ^ are the 
basic equations of the statistical model as used in the 
present paper. They are used throughout except for Sec- 
tion 1vIlX21 Isobaric analogue resonances occur for pro- 
ton energies below or at the Coulomb barrier where the 
energy dependence of the matrix elements Wc^{E) can- 
not be neglected. 

In some cases it is useful to write S'~^^ in the diagonal 
representation of H'^'~'^. We diagonalize every realiza- 
tion of ijCiOE -^j^j^ g^j^ orthogonal matrix O. The eigen- 
values Ef^, /i = 1, . . . , follow the Wigner-Dyson distri- 
bution, see Section I.II.D. For A^ — > oo the elements of O 
arc uncorrelated Gaussian-distributed random variables 
and so arc, therefore, the transformed matrix elements 
W^a = Y^bOfi„W„a- These are not correlated with the 
-B^s, have zero mean values and have second moments 
given by 

{Wf^aW.b) = v%b5^. , (48) 

see Eq. fSSl) . It is assumed that the v'^s are indepen- 
dent of energy. In that representation the matrix S'^^ of 
Eq. (HH) is given by 



S^i'iE) = 6ab - 2inY,Wat,iE)iD-^)^uW,biE) (49) 



where 



(50) 



As done for S"'"" in Eq. {351), in Eq. gj) can also 
be written in terms of the if-matrix. That form was, 
in fact, the starting point of some work on the statis- 
tical theory reviewed in Section IIII.CI The forms (|4T1) 
and (|35p . although mathematically equivalent, differ in 
one essential aspect. The S'-matrix in Eq. ^^\\ depends 
explicitly on H'^'^^ . Ensemble averages of moments of 
as given by Eq. flTl) are calculated by integrating 
over the Gaussian-distributed matrix elements of ff^^^. 
The orthogonal invariance of the GOE makes it possible 
to do the calculation analytically for some of these aver- 
ages, see Section IVl In contrast, the AT-matrix depends 
on the eigenvalues and eigenvectors of H'^'-'^. Thus, the 
orthogonal invariance of the GOE is not manifest in the 
form (1211) , and ensemble averages require a separate inte- 
gration over the distribution of eigenvalues and of eigen- 
vectors of the GOE. The calculation turns out to be pro- 
hibitively difficult. That is why the AT-matrix formalism 
has only been used for numerical simulations. 
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We have derived Eqs. (|40I41I42I) and (gT]) in the 

framewo rk of the shell-model appr o ach t o nuclear re- 
actions ( Mahaux and Weidenmiilleil . Il969| ). That ap- 
proach directly yields the dependence of the scattering 
matrix on the Hamiltonian for the quasibou nd states . 
Feshbach's u nified theory of nuclear reactions ( Feshbachl . 
Il958l Il962l I1964D yields similar but more formal ex- 
pressions written in terms of the projection operators 
onto the closed and the open channels. To actually im- 
plement an RMT approach into these expressions, and 
to work out averages from the resulting formulas, one 
has to write Feshbach's ex pressions explicitly in term s 
of nuclear matrix element s (iLemmer and Shakin . 1964[) . 
Such an approach yields ( Lewenkopf and Weidenmiilleil . 
Il99l[ ) formulas quite similar in structure and content 
to Eqs. (j40l41l42p and (gT]). The removal of direct re- 
action contributions by the matrix U in Eq. (j40p and 
the reduction of the 5 -matrix to the canonical form 
i n Eg (gT |) goes back to Engelbrecht and Weidenmiillerl 
(|1973[ ) and lNishioka and Weidenmaileil (|l985[ l 

Averages of observables are theoretically worked out 
as averages over the ensemble of scattering matrices, i.e., 
over iJ^^^, and are denoted by the same angular brack- 
ets as used to indicate energy averages. The theoretical 
result (A) for an observable A is compared with the ex- 
perimental running average over energy of the same ob- 
servable measured for a specific nucleus, i.e., for a specific 
nuclear Hamiltonian. It is stipulated that that nuclear 
Hamiltonian is a member of the GOE. Ergodicity would 
then guarantee the equality of both averages, see Sec- 
tion I.H.C.3. Ergodicity holds true for A if the energy- 
correlation function of A goes to zero for large values 
of the argument. In the case of the scattering matrix 
S™ , that property is analytically fully established for 
all observables formed from S'~^^ only in the domain of 
strongly overlapping resonances. There is no reason to 
doubt, however, that ergodicity applies also in the regime 
of isolated and weakly overlapping resonances. 

As discussed in I.II.C.2 the local spectral fluctuation 
measures of the GOE do not depend on the Gaussian 
distribution assumed for the matrix elements H^^^ in 
Eq. p6p but apply for a wide class of random-matrix en- 
sembles. That property is referred to as universality. The 
fluctuation properties of the 5'-matrix are studied over 
energy intervals that are measured in units of the mean 
level spacing of the resonances. The relevant observables 
(S-matrix correlation functions) are, thus, likewise lo- 
cal fluctuation measur es. The proof of universality of 
such m easures given by iHackenbroich and Weidenmiillerl 
( 19951 ) applies also to these observables. 

The theoretical framework of Eqs. (|40l41l42p and (|T7)) 
is quite flexible and allows for extensions of the the- 
ory. These may account for violation of isospin sym- 
metry, of parity, or for tests of time-reversal invari- 
ance in CN reactions. Such applications of the statis- 
tical theory are reviewed in Sections IVII.Al IVII.Bl and 
IVII.Cl Violation of isospin symmetry or of parity is 
treated by replacing ijGiOE (jTT)) by a block ma- 



trix ( Rosenzweig and Ported . [r960[) as done in Eq. (1.30) 
and in Eq. ([M]) below. Each diagonal block refers to 
states with the same isospin (or parity) quantum num- 
bers while the off-diagonal blocks contain the matrix el- 
ements of the symmetry-breaking interaction, see Sec- 
tion I.III.D.l. Time-reversal invariance is broken when 
jjGOE j-gpiaced by i^GOE _|_ where a is a real pa- 
rameter and ^ is a real and antisymmetric random ma- 
trix, see Eqs. (1.32) and (1.33) and Section IVlI. CI below. 

The statistical theory for S formulated so far applies 
to every set of flxed quantum numbers (total spin, par- 
ity) of the CN. In keeping with Eq. ([5]) we assume that 
S'-matrices referring to different quantum numbers are 
statistically uncorrelated. That assumption is not as 
innocent as it may look. Indeed, let us imagine that 
in Eq. is determined not from the GOE but 
from the shell model. That calculation would yield the 
matrices H^^, (and, thereby, the S'-matrices) simultane- 
ously for all conserved quantum numbers. A change of 
the residual interaction of the shell model would cause 
all these matrices to change simultaneously. Different 
realizations of the GOE may be thought of as corre- 
sponding to different choices of the residual interac- 
tion. Therefore, the matrices i/GJO^ appearing in S'- 
matrices carrying differen t quantum numbers arc ex- 
pected to be correlated (jMulhall et al I I2000D . Such 
correlations among Hamiltonian matrices referring to 
stat es with different quantum nunibers do indeed ex- 
ist ( Papenbrock and Weidenmiilleil . l2007t ) and were dis- 
cussed in I.V.B.5. To what extent is the assumption for- 
mulated in Eq. (O invalidated by the existence of such 
correlations? As remarked in Section lll.Ci Eq. ([5]) implies 
that CN cross sections are symmetric about 90 degrees 
in the cm. system. The available experimental evidence 
supports that prediction and we expect, therefore, that 
deviations from Eq. ([5]) are not signiflcant. We return to 
that point in Section IV.GI 



C. History 



The theory developed in Sections IIII.AI and IIII.BI may 
appear quite natural. However, it took several decades 
to arrive at that formulation. Since the phenomeno- 
logical models reviewed in Sections III. CI and III.DI led 
to predictions that were in good agreement with ex- 
periment, a deeper understanding of these models was 
called for f rom the outset. By way o f justific a tion o f 
their work, Brink an d Step hen! ( 19631 ): lEricsonI (I960); 
iHauser and Feshbach (1952.1 themselves had referred to 
the statistical properties of isolated resonances (which 
were then supposed and are now known to agree with 
GOE predictions). But the actual derivation of the 
Hauser-Feshbach formula and of Ericson fluctuations 
from a GOE model for CN resonances posed a severe 
challenge, especially if the aim was a comprehensive sta- 
tistical theory based on the GOE that would apply for 
all values of the parameter T /d. That resulted in a large 
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number of publications in the years 1950 to 1985. Here 
we can give only an outline of the main developments. 

The development of the statistical theory of CN reac- 
tions depended on the availability of a suitable theoret- 
ical framework to describe resonance reactions. Such a 
framework is needed to formulate statistical assumptions 
on the resonance parameters. But to a large extent the 
availability of such a framework depended on the devel- 
opment of nuclear-structure theory. Early attempts to 
formulate statistical models mirror the development of 
the dynamical theory of nuclei. 

The i?-matrix theory by I Wiener and EisenbudI (I1947D 
was formulated at a time when virtually nothing was 
known about nuclear structure and is, by necessity, a 
very formal theory: Except for the short range of the 
nucleon-nucleon interaction, it docs not refer to any spe- 
cific feature of the nuclear Hamiltonian. Resonances are 
constructed as follows. The nucleus is thought to be en- 
closed by fictitious boundaries (one for each two-body 
fragmentation). On these boundaries fictitious bound- 
ary conditions arc imposed. As a result, the spectrum 
of the nuclear Hamiltonian in the "internal region" (the 
domain enclosed by the boundaries) is discrete. The 
eigenvalues and eigcnfunctions depend on numerous pa- 
rameters (distance of the boundaries from the center of 
mass, values of the boundary conditions). Green's theo- 
rem is used to connect these discrete levels with the chan- 
nels, and the levels become CN resonances. Even today 
the resulting form of the scattering matrix is extremely 
useful for the analysis of experimental data containing 
several partly overlapping resonances ("multi-level R- 
matrix fit"). However, few-level approximations to the 
S'-matrix are not automatically unitary. Moreover, the 
form of the S'-matrix is rather unwieldy and depends 
explicitly on the parameters just mentioned. For these 
reasons, the i?-matrix has never been used as the starting 
point of a statistical approach. 

iLane and LvnnI ( 19571 ) calculated the average cross sec- 
tion for neutron capture reactions in the regime of iso- 
lated resonances. The calculation was possible without 
referring to a comprehensive theory of resonance reac- 
tions. Indeed, in the regime F <C d each resonance is 
independently described by a Breit-Wigner formula as 
in Eqs. or (^5)) . If one assumes that the distribution 
of resonance parameters of isolated resonances follows the 
GOE, the partial width amplitudes in each channel are 
Gaussian-distributed random variables. For each reso- 
nance the total width F^ is related to the partial widths 
Tafj, by the sum rule (jl9p . The resonance energies do 
not enter the calculation of the energy-averaged cross 
section. The task consists in calculating the ensemble 
average of Fa^Fb^/F where a and b are two specific chan- 
nels (here: the neutron and the gamma decay channel). 
This can be done using the GOE and as input values for 
the averages of the partial widths, the strength functions 
{^ati)/d, see Section FV.Al A correction factor deduced 
from i?-matrix th eory accounted fo r the effect of weak 
resonance overlap. iMoldaueil ( 196ll ) extended the result 



of Lane and Lynn pcrturbatively to weakly overlapping 
resonances. 

The first attempt t o construc t a comprehensive statis- 
tical t heory is due to iMoldaueil (|l96ll . Il963l 1 1964 Il969l 
Il975l . Il976l Il980l) (see also references to further work 
therein). In order to go beyond the regime of isolated 
resonances, Moldauer used an expansion of the scatter- 
ing matrix in terms of its poles in the complex energy 
pl ane. This form of the S-mat rix had been proposed 
bv lHumblet and Rosenfeldl (|l96l[ ). In contrast to the R- 
matrix theory, the Humblet-Rosenfeld theory is a dy- 
namical theory: The parameters of the pole expansion 
are, in principle, completely determined by the Hamil- 
tonian, and there are no arbitrary parameters. Unfor- 
tunately an explicit analytical connection between the 
Hamiltonian and the positions of the poles and the val- 
ues of the residues, is not known. Thus, the Humblet- 
Rosenfeld theory is effectively a formal theory like the R- 
matrix theory. Moreover, the following diflaculty arises. 
For isolated resonances there is a one-to-one correspon- 
dence between a resonance and a pole of the scattering 
matrix, see Eq. ((T8|) . The parameters of the pole possess 
a direct physical interpretation as partial width ampli- 
tudes and as energy and total width of the resonance. 
In the general case of many overlapping resonances the 
Humblet-Rosenfeld expansion uses the same parameters 
(residues and locations of the poles of the S'-matrix plus 
a smooth background matrix). But the simple interpre- 
tation valid for isolated resonances does not apply. More- 
over, the theory is not manifestly unitary. For isolated 
resonances the unitarity constraint leads to the simple 
sum rule (|19p , while it imposes complicated relations be- 
tween the pole parameters for overlapping resonances. 
These have never been untangled. Therefore, the choice 
of statistical assumptions for the pole parameters was 
far from obvious. At some point, Moldauer used the K- 
matrix formulation (|35p of the S-matrix, assumed that 
the energies i?^ and matrix elements Wafj. in Eq. ((34)) 
obeyed GOE statistics, and determined the distribution 
of pole parameters via a numerical simulation. His re- 
lentless efforts met with limited success but kept interest 
in the problem. Some of his results are reviewed in Sec- 
tion |VT1 

Fesh bach's u n ified theory of nuclear reac- 
tions (iFeshbachl Il958l I1962L I1964D was the first 



theory of GN reactions that expressed the scattering 
matrix in terms of the Hamiltonian of the system with- 
out the help of arbitrary parameters. It uses projection 
operators onto the spaces of open and of closed channels. 
The projection of the nuclear Hamiltonian onto the 
space of closed channels defines a self-adjoint operator 
with a discrete spectrum. The bound states of that 
operator generate the r esonances in the fuU prob lem. 
The theory was used by iLemmer and ShakinI (I1964D for 
a first calculation of elastic neutron scattering on ^^N 
using the nuclear shell model. Quasibound shell-model 
states appeared as neutron resonances. Although not 
related to the statistical theory, that work demonstrated 
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the possibihty to account for CN resonances in terms 
of a dynamical approach using the nuclear s hell m odel. 
Later work by the MIT group ( Kawai et all Il973t ) used 
Feshbach's theory to formulate an extension of the 
Hauser-Feshbach formula that accounts for the presence 
of direct reactions, see Section IILEI The approach 
uses plausible assumptions but is not based upon a 
random-matrix approach. 

The general theory of n uclear resonance reactions 
develo ped in the 1960s by iMahaux and Weidenmiillerl 
(I1969D was based on the nuclear shell model and work 
by iDirad (I1958D and iFanol (|l96lD . It was more explicit 
than Feshbach's theory of nuclear reactions. That was a 
natural consequence of the fact that nuclear-structure 
theory had made significant progress since the early 
1950s. Resonances were shown to be mainly due to qua- 
sibound states of the shell model. The approach ex- 
presses the scattering matrix in terms of the Hamilto- 
nian governing the dynamics of the quasibound states as 
in Ec^. ([28| . That was the starting point of later devel- 
opments. The theory is connected with the i?-matrix 
theory of I Wiener and E iscnbud (1947,) through relations 
of the form of Eq. ((35l) . 

Indep endently of that develop- 
ment, lEngclbrecht and Weidenmiilleil ( 1973f ) showed 
that there exists a unitary transformation which diago- 
nalizes (S). The transformation acts on channel space 
and, therefore, leaves the statistical p r opertie s of the 
resonances unchanged ( Hofmann et all Il975bl ). As a 
consequence, the transformation reduces the problem 
of calculating CN cross sections in the presence of 
direct reactions to the problem without direct reactions 
(i.e., for a diagonal (S)). That was an important 
simplification, see Eq. pO|) . 

The conspicuous lack of a comprehensive t heory of CN 
reacti ons with predictive power m otivated iTepel et all 
(|1974| ) and lHofmann et al. r (|l975allbl ) to perform numer- 
ical simulations with the aim of establishing fit formulas 
for average CN cross sections valid for all values of F/d. 
The authors used the transformation mentioned in the 
previous paragraph and focused attention on 5*^^ . They 
proved that aside from overall phase factors, the distri- 
bution of ^-matrix elements depends only on the trans- 
mission coefficients. That fact simplified the construc- 
tion of the fit formulas. For the numerical simulations, 
the authors used the if -matrix form psp of the stochas- 
tic scattering matrix (|4ip . The A"~matrix parameters 
were determined by the eigenvalues and eigcnfunctions of 
GOE matrices. Repeated random drawings from a Gaus- 
sian distribution of the elements of GOE matrices led to 
statistically meaningful averages and determined the pa- 
rameters in the fit formulas, see Section IV.D.1[ Similar 
formu las were subsequently also developed by Moldaueil 

dHzi. 

The first use of the explicit dependence of Eq. p9)) on 
the Hamilto nian H in the calcu lation of CN processes 
was made bv lAgassi et al\ ( 1975[) . CN scattering theory 
was extended so as to include precompound reactions. 



The authors used a stochastic model for the matrix H 
in Eq. ([^^ which allowed for the existence of classes of 
shell-model states of increasing complexity. Averages of 
cross sections were calculated in terms of an asymptotic 
expansion valid for F 3> d. The Ericson regime was ob- 
tained as a special case of very strong coupling between 
classes so that the internal equilibration time of the CN 
becomes small compared to the decay time Ti/T. As a 
result, the Hauser-Feshbach formula and predictions of 
the Ericson model were for the first time derived from a 
microscopic statistical theory. That was possible because 
resonance spacings were assumed to be constant (i.e., not 
to follow GOE predictions). Later work showed that for 
F ^ d that simplification led to correct results. 

In 1984, a c o nnectio n was estab- 
lished (jVerbaarschot et all Il984al) between statistical 
nuclear theory and quantum field theory. More precisely, 
field-theoretical concepts used in condensed-matter 
theory and in statistical mechanics were applied to the 
RMT description of both, nuclear spectra and nuclear 
reactions. In a very concrete sense, RMT as applied 
to nuclei became part of the statistical mechanics of 
many-body systems. Specifically, generating functionals 
familiar from condensed-matter theory were used to 
describe fluctuations both of nuclear spectra and of 
^-matrix elements in a common framework. In the limit 
iV ^ 1, the evaluation of the generating functionals 
permitted a clear separation of average properties of 
observables and of their fluctuations. The stochastic 
scattering matrix was deflned in the form given by 
Eq. (|41[) and was written as a suitable derivative of 
a generating functional Z. As a flrst application, the 
^-matrix correlation func tion was calcu l ated w ith the 
help of the replica trick ( Weidenmiilleil . Il984l ). That 
yields an asymptotic expansion in powers of d/F, see 
Section IV.BI Later, Z was ex pressed in te rms of 
Efet ov's supersymmetry a pproac h (lEfetovl . Il983l ) . That 
led (|Verbaarschot et all Il984bl Il985[) to the exact 
expression for the iS-matrix correlation function given 
in Section NT\ Eqs. (|40I41I42|) and g7|) and the 
supersymmetry approach have since been established 
as the main tools to study chaotic scattering based on 
RMT. 



IV. AVERAGE 5-MATRIX 

The average S-matrix serves as input for the statis- 
tical theory of nuclear reactions. In the present Section 
we display central properties of {S) and review ways to 
determine {S) phenomcnologically. 



A. Calculation of {S'^^ {E)) 

We focus attention on the only stochastic element in 
the scattering matrix, i.e., on the CN part (S™), see 
Eq. (pij) . The calculation of averages of S™ (and of 
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powers of S'~^^) is simpler than that of terms involving 
both S™ and {S'^^y because as shown in Eq. §^ aU 
poles of S'~^^{E) lie below the real E-axis. 

We calculate {S'^'^{E)) by replacing the ensemble av- 
erage by a running average over energy with a Lo rentzian 
weigh t function of width / centered at energy Eq ( Brownl . 
Il959l ). We choose the center of the GOE spectrum 
iJo = as the center of the weight function and have 



{S™iE)) = - / dE S^'^iE) 

J-oo 



E^ + {l/4)P 



(51) 



The average in Eq. (jSip extends over very many reso- 
nances so that I ^ d with d the average level spacing in 
the center of the GOE spectrum. The integration con- 
tour can be closed in the upper half of the energy plane. 
We obtain 



(52) 



We use Eqs. (|49|) and (HO]) and expand S™{il) in a Born 
series with respect to the imaginary part of , 

oo 

(5™(i?)),b = 5,b + 2^[(-z^W)"]ab . (53) 

n=l 



Here W is a matrix in channel space given by 



fib • 



(54) 



Since / ^ d, it is legitimate to neglect the fluctuations 
of the Efj_ in the denominator in Eq. i.e., to as- 

sume for the Efj^ a picket-fence model with fixed nearest- 
neighbor spacing d. The sum over fi amounts to averag- 
ing the product Waf^Wf^b- We have (1/iV) Wa^^Wf^t = 
(l/N) WaiiW^b- We use Eq. ([38]), change the summa- 
tion into an energy integration, and find 



Wab 



-ITTV 



(55) 



Inserting the result into Eq. (|53p and using the defini- 
tion p5)) . we find 



1 



— Oab 



(56) 



Eq. ((56)) relates the average scattering matrix to the pa- 
rameters Xa of the statistical theory. The result ([55)1 
is also obtained in the limit iV — >■ oo via the replica 
trick (see Section [V.B.2|) and via the supersymmetry ap- 
proach (see Section IV. C[) . The average ^-matrix given 
by Eq. ([56]) is real because we have chosen E in the center 
of the GOE spectrum, E = 0. For a different choice of E, 
i^aa) would not be real. The additional phase caused 
by such a choice is not part of the matrix U constructed 
in Section UlLBl 

Eq. ([55)) shows that (S™) is always diagonal. But 
that does not imply that (S) is diagonal, too. In fact. 



Eq. ((40)) demonstrates that (5*) is diagonal only if the 
matrix U is diagonal. That matrix embodies the effect of 
couplings between channels. That is why a non-diagonal 
form of (S) is taken to be synonymous with the presence 
of "direct reactions" (i.e., reactions that proceed without 
intermediate CN formation). 

The average of a product of ^-matrix elements can 
be worked out in the same way. Let {ai,bi} with i — 
1, . . . , fc denote k arbitrary pairs of channels, and con- 
sider the average of the product of the matrix elements 
S™. {Ei) taken at energies Ei, E2, ■ ■ ■ , Ek- Again using 
a Lorentzian averaging function, closing the contour in 
the upper half of the complex energy plane, and using 
Eq. dSl]) we find 



(57) 



More generally, we have for any function f{S™) that is 
analytic in the upper half of the complex energy plane 
the relation 



CN\ 



(58) 



This is an important result. It shows that S'-matrix ele- 
ments taken at the same or at different energies are un- 
correlated. The result ((57)) does not hold for products 
involving both and (5''-'^)* since these have poles 
on both sides of the real energy axis, and closing the 
integration contour does not yield a simple expression. 
The calculation of such terms poses the main technical 
difficulty in the statistical theory. 



B. Physical Interpretation of (S) 

The decomposition ([3]) of the scattering matrix into an 
average part and a fluctuating part applies likewise to the 
CN part S'CN of S. The discussion in Section [iLFj of the 
physical significance of that decomposition in terms of 
time scales for the CN reactions then applies to S™ , too, 
as do the limitations of the statistical theory established 
there. 

Eq. ((56)) shows that the average CN S'-matrix is sub- 
unitary, I (5*™) I < 1. The equality sign holds only when 
Xa vanishes, i.e., when there is no coupling between chan- 
nel a and the N resonances. We emphasize that the lack 
of unitarity of (5^^) is due to taking the average over 
N resonances (and, since {S'~^^) is diagonal, not to in- 
elastic scattering processes that would deplete the elastic 
channel) . To see how averaging reduces the magnitude of 
S'~^^ we consider the simplest case of a single channel cou- 
pled to N resonances. There are no inelastic processes 
by definition. The scattering amplitude has the form 
S'~"^{E) = exp(2iS{E)) and can be viewed as a point on 
the unit circle of the complex plane. As E increases, 
6{E) increases by tt over the width of every resonance, 
and S™ {E) moves counter-clockwise once around the 
origin on the unit circle. The average of S'~^^{E) over 
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N resonances must then lie in the interior of the unit 
circle, i.e., be subunitary. These points are b eautifully 
illuminated in ( Friedman and Weisskopi Il955[ ) . 

It is remarkable that (S^^) as given in Eq. ([55]) de- 
pends only on the coupling coefficient Xa and not on the 
other coupling coefficients xt with b ^ a. That shows 
that (S^^) taken all by itself describes the loss of prob- 
ability amplitude in channel a due to CN formation. In- 
formation on the manner in which the CN eventually de- 
cays back into the various channels, is not contained in 
i^aa)- The information as to where the lost probability 
eventually reappears (in quantum mechanics, probability 
is conserved) is supplied by the fluctuating part of . 
Thus, CN scattering theory can b e viewed as a spec ial 
case of quantum transport theory ( Agassi et a/.l . Il975!). 

In Eq. ([S]); the transmission coefficients were defined 
under the assumption that (S) is diagonal. That assump- 
tion is met by (S'^'^). We calculate Ta from Eq. ^ using 
Eq. dSSl) for (S-CN) and find 



Ta 



(59) 



(For a non-diagonal (S) the connection with the trans- 
mission coefficients is given in Section lIV.Cp . The 
form (j59p for the transmission coefficients applies when 
we choose the energy E in the center of the GOE spec- 
trum, E = 0. For a different choice of E, the form of Ta 
differs. That difference is irrelevant in practice because 
the transmission coefficients serve as phenomenological 
input parameters of the theory anyway. 

In the theory of resonance reactions formulated in Sec- 
tions nil. A. 21 and IIII.BI the coupling between levels and 
channels seems to increase monotonically with increas- 
ing strength of the coupling matrix elements Wafj, or, 
equivalently, of the parameters Xa defined in Eq. (|45|) . 
However, Eq. ((5^ shows that this is not the case. With 
increasing Xa, the transmission coefficient Ta increases 
monotonically until it reaches its maximum value = 1 
at Xa = 1. Thereafter Ta decreases monotonically with 
increasing Xa and tends towards zero as Xq — > oo. Simi- 
larly, Eq. dsn]) shows that {S^^) decreases from unity to 
zero as Xa approaches unity from zero. As Xa increases 
further, (S™) approaches minus one. We also note that 
under the substitution Xa — I/xq, Ta is invariant while 
(5™) changes sign, {S^!'){l/xa) - -(5™)(^a). Ah 
this shows that there is a maximum value of the cou- 
pling strength, Xa = 1, beyond which a further increase 
of Xa effectively reduces the coupling between channel a 
and the resonances. 

To make that fact physically plausible we consider 
th e matrix D in Eq. (H^. We follow the work 
of ISokolov and Zelevinsk^^l988L I1989L I1992L I1997D . For 
large values of the Xa (all a), the width matrix F^^ 
(Eq. (gni)) dominates F^OE jjoS ^ ^ There- 

fore, one chooses a basis in which F^^ is diagonal. 
Eqs. ([55]) show that F^i, possesses A non-zero eigenval- 
ues 2ttNv'^ with a = 1, . . . , A. The associated orthonor- 
mal eigenvectors are given by (l/-\/]Vtjf)VFap. Com- 



pleting in some arbitrary fashion these A eigenvectors 
to a set of N orthonormal vectors, denoting the result- 
ing orthogonal matrix by G and transforming the matrix 
D = GDG'^, we obtain 



D., 



{E + inNvl)8^, 



(60) 



where H = GH'^^^G^ and where the first A diagonal 
entries ttNv'^ are given in terms of the A non-vanishing 
eigenvalues 2'kNv'^ with a ~ 1 , . . . , A of the width ma- 
trix V^i, while = for ^ > A. Concerning these first 

A resonances, we consider H as & small perturbation. 
That is justified as long as the decay widths 27riVw^ of 
the resonant states /i = 1 , . . . , A are large compared with 
their spreading widths 27r([if^i,]^)/c? due to mixing with 
the other resonances. (We recall that the concept of the 
spreading width was introduced in the context of door- 
way states in Section I.II.G. The entire discussion there 
applies likewise to the present case). With d = itX/N 
and with the help of Eq. ([5^ that implies ^ d or 
Xa ^ I for all a. In zeroth order in H we then deal 
with A poles located at —iirNv^, /i = 1, . . . , A and with 

— A poles located on the real axis. Taking into ac- 
count H we obtain a "cloud" of A^ — A poles located 
below but close to the real axis and A poles very far from 
that axis. Each of the latter contributes only a smooth 
phase to the ^-matrix (smooth over a typical distance 
d). The additional phase causes the minus sign in the 
relation {S^^){l/Xa) = -{S^^)ixa). Aside from that 
phase, the A far-away poles do not affect the CN scat- 
tering process. On the other hand, CN scattering due 
to the cloud of A^ — A narrow resonances is, for A fixed 
and N ^ oo, indistinguishable from that of A^ narrow 
CN resonances. That is why, aside from overall phase 
factors, all moments and correlation functions of S*^^ do 
not depend on the phase shifts due to the A far-away 
poles and can be expressed in terms of the transmission 
coefficients Ta- The picture of the distribution of poles 
that emerges from this qualitative discussion is quantita- 
tively confirmed in Section lV.FI As all the coefficients Xa 
increase monotonically from zero, the A^ poles of S move 
away from the real energy axis into the complex plane. 
As the Xa increase beyond Xa ^ 1 (all a), the "cloud" of 
N poles begins to separate into a small cloud containing 
A poles that move ever further away from the real axis, 
and a big cloud of A^ — A poles that moves back toward 
the real axis, see Fig. [5j 

According to Eq. ([5]), the transmission coefficients 
measure the unitarity deficit of the average S'-matrix. As 
mentioned in Section [II. C[ they also measure the proba- 
bility of CN formation. It is now clear why the T^s are 
better suited for that purpose than the coupling coeffi- 
cients Xa- 



C. Satchler's Transmission Matrix 



In Section IIV.BI we have used S™ (which is diago- 
nal on average) to express (S^^) and Ta in terms of the 
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coupling coefficients Xa, see eqs. ([56)1 and ((59)) . We now 
address the question how to determine the transmission 
coefficients if (S) is not diagonal. That is the generic 
case. In principle, the question can simply be answered 
with the help of Eq. (|40]): To find the TaS, take the ma- 
trix S™ = U'fSU*. By definition, the average of 
is diagonal, and Eq. ([8]) yields the transmission coeffi- 
cients. But historically the question came up before the 
form (|40)) with U defined in Eq. (|39)) was known. It 
arose in conjunction with the experimental discovery of 
"direct reactions" mentioned in Section III.Ei and with 
the theoretical treatment of such reactions which are de- 
scribed by non-diagonal average scattering matrices, see 
Section HfcH 

A general measure of the unitarity d eficit of (S) is gi ven 
by Satchler's transmission matrix P ( Satchlerl . Il963[) . a 
matrix in channel space defined in Eq. (flT)) . Since P is 
Hermitean, it can be diagonalized by a unitary matrix U. 
We denote the eigenvalues by Pa with a = 1 , . . . , A and 
have 



where < Pa 



< 1 



(61) 



Applying the transformation U to the right-hand side 
of Eq. ([T7| and using U*U'^ = 1a we are led to con- 
sider the matrices A = W {S)U* and A* = U^{S}IJ. 
Eq. (EI]) implies that {AA*)ab = 5ab{l - Pa) = iA*A)ab 
which in turn shows that the symmetric matrices A 
and A * are normal and can be diago n alized simultane- 
ously ( Engelbrecht and Weidenmiilleil . 11973) . The ma- 
trix product AA* being diagonal already, we conclude 
that A and A* are diagonal, too (except for accidental 
degeneracies of the pa which we do not consider). Thus, 



iUHS)U*)ab - SabV^^aCMMa) 



(62) 



where the 4>aS are real. Defining Uab — UabCyi'p{i4>b) we 
conclude that the average of the matrix U"^ SU* is real 
and diagonal, {If^ SU*)ab = SabV^ ~ Pa- That shows 
that the transmission coefficients are given by the eigen- 
values Pa of Satchler's transmission matrix, Ta = Pa, o, = 
1, . . . , A. Moreover, the function in Eq. (j^H) is in- 
variant under an arbitrary unitary transformation in the 
space of channels. Therefore, the matrix = W SU* 
is identical with the scattering matrix of pure CN scat- 
tering in the absence of direct reactions. In other words, 
the matrix U reduces the CN scattering problem in the 
presence of direct reactions to the pure CN scattering 
problem without direct reactions. 

The transformation U introduced 

by lEngelbrecht and Weidenmiilleil ( 19731 ) and de- 
fined in terms of S atchler's transmission matrix is a 
purely formal device. iNishioka and Weidenmiilleil ( 19851 ) 
have shown how U is defined physically. We do not 
retrace the steps of the argument here because we have 
used it to construct the matrix U in Eq. (p9)) . The 
present construction of the matrix U via Satchler's 
transmission matrix yields exactly the same matrix U 
as in Eq. ([551) if the energy E is chosen in the center 



of the GOE spectrum, so that E = 0. Otherwise, the 
matrix U obtained from Eqs. ((CT|) and ([5^ contains an 
additional phase (the phase of the average S-matrix, see 
the remark below Eq. ((SS]) ). 



D. Optical Model and Strength Function. Direct Reactions 

The average S'-matrix serves as input for the statisti- 
cal theory and must be determined phenomenologically. 
It is intuitively clear that the fast part (S) of the reac- 
tion amplitude S can involve only few degrees of freedom. 
Therefore, models that do not take into account the full 
complexity of the many-body problem suffice for an ac- 
curate determination of (S). These are the optical model 
of elastic scattering and the coupled-channels approach. 
The latter describes fast inelastic processes that are re- 
ferred to as direct reactions. We do not review these 
models in detail but confine ourselves to those aspects 
that are essential for the understanding of the statistical 
theory. 

In the optical model ( Hodgsoiil . ll963[ ). elastic scattering 
of nucleons is described in terms of a radial Schrodinger 
equation with a central potential V{r) + iW{r). The bulk 
of the real part V{r) is due to the shell-model potential. 
(The shell model is the subject of Section I. IV. A). The 
imaginary part W{r) describes CN formation by nucleon 
absorption. Via a dispersion relation, the imaginary part 
W{r) also modifies the real (shell-model) part V{r). Be- 
cause of the presence of W{r), the scattering amplitude 
calculated from the optical model is subunitary and is 
identified with the average scattering amplitude in the 
nucleon channel of the statistical theory. Obvious gener- 
alizations apply to the scattering of composite particles. 
Fig. d] shows results of an optical-model calculation. 

As mentioned in Section fllLCi the CN was originally 
considered a black box. The novel asp e ct of the opti- 
cal model introduced by iFeshbach et al. 1 (Elll was the 
partial transparency of the target nucleus. That modi- 
fied Bohr's picture of the CN as a system of strongly in- 
teracting particles: Nucleons in the nucleus had a finite 
mean free path. The nuclear shell model had indicated 
such behavior already for the ground state and low-lying 
excited states. Now that feature was extended to reso- 
nances (states above neutron threshold). The influence of 
shell structure on the CN cross section is manifest in the 
neutron strength function. In analogy to the spreading 
width introduced in Section I.II.G, the neutron s-wave 
strength function s{E) is defined as s{E) = 27r(r^a)/d. 
Here is the partial width for s-wave neutron emis- 
sion of resonance ^. With F^a ~ SttVF^^ and Eq. ((45)) 
that yields s{E) = Axa and, for Xa <. 1, s{E) = Ta, 
see Eq. ([5^ . The strength function s{E) measures the 
intensity with which a group of resonances is coupled to 
the s-wave neutron channel. The strength function can 
both be measured (by performing a running average over 
a number of resonances) and be calculated directly from 
the optical model. We recall that the radius of the real 
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part of the optical-model potential (i.e., essentially the 
shell-model potential) increases with A. The strength 
function displays maxima versus A whenever a single- 
particle s-wave state is about to be pulled into the po- 
tential well. 

The phenomenological description of direct reactions 
is based upon an extension of the optical model for 
elastic scattering. For a set of channels, the radial 
equations describing elastic scattering with the help of 
the optical model, are coupled. The term which cou- 
ples two channels is obtained as the matrix element 
of the interaction between projectile and target sand- 
wiched between the two channel wave functions, see 



the reviews bv lAusternl (Il970l) : iGlendenninS (|2004[ ): 
iKoning and Delarochd (|2003l ). Direct reactions are most 
important for strongly coupled channels. Often these are 
channels where the target nucleus is in different states 
of collective excitation. Nucleon transfer between pro- 
jectile and target may also lead to strong interchannel 
coupling. In many cases, the solution of the coupled- 
channels problem can be simplified by using the Born 
approximation. The result is the distorted-wave Born 
approximation: The plane waves in entrance and exit 
channel are distorted by the optical potentials, the tran- 
sition between both channels is calculated perturbatively 
to first order. 

Direct reactions induce correlations between 
parti a.l-width arn p litude s relating to differ ent chan- 
nels (jHafner et ai Il967l: iKawai et all Il973l) . In the 
present context, that is seen as follows. Eq. pS)) 
shows that for a ^ b, the partial-width amplitudes 
V^nWfj^a and ^/2nWfj,b are uncorrelated random vari- 
ables. But in the presence of direct reactions, the 
matrix U in Eq. ([101 is not diagonal, and the partial- 
width amplitudes of the full S matrix in channel 
a are given by ^f,Uab\/27rWf_,b- For two channels 
a b connected by a direct reaction these are, in 
general, correlated. Indeed, using Eq. we find that 
(7^a7;6>=27rEcf^acC/fe>cV0. 



V. RESULTS OF THE STATISTICAL THEORY 



According to Eqs. (|40l41l42p . S{E) is a matrix-valued 
random process. The goal of the statistical theory con- 
sists in finding the joint probability distribution of S{E) 
and S* [E) . Because of Eq. (|40)) it would actually suffice 
to determine the joint probability distribution of S'^^^ (E) 
and {S'-^^{E))*. That distribution should be given in 
terms of the transmission coefficients Tc, of energy differ- 
ences in units of the mean GOE level spacing d, and 
overall phase factors. We are far from that goal be- 
cause averaging over the N{N + l)/2 random variables 
of ijGOE ii^g^^ appear in the denominator of S'~^^ (see 
Eqs. (|41l42p and d!?])) turns out to be extremely difficult 
in general. We possess only partial information, mainly 
on low moments and correlation functions of S'~^^ (E) and 
{S™{E))*. That information suffices for a comparison 



with the available experimental data and is now reviewed, 
together with the various methods that have been used 
to calculate the answers. We do not go into the full com- 
plexity of some of the calculations as these are technically 
quite demanding, and refer the reader to the original lit- 
erature. 

In the calculation of moments and correlation func- 
tions, a simplification arises because the coupling matrix 
elements obey Eq. ((38|) . It follows that all moments and 
correlation functions of S™{E) and {S'^^{E))* vanish 
unless the channel indices are pairwise equal. For in- 
stance, using the unitary transformation of Eq. (j40p to 
calculate (|S'obp), we find that only terms of the form 
{\S™\^) and {S^^{S^3)*) give non-vanishing contribu- 
tions. 



A. Isolated Resonances (F <C d) 

Isolated resonances occur in two limiting cases, Xa <C 1 
and Xo ^ 1 for all a. The two cases differ only in the 
phases of the average 5'-matrix elements and lead to iden- 
tical expressions for the CN cross section. Therefore, we 
consider only the case Xa I for all a. We use the 
form (j49l50p of S™ and the statistical properties listed 
above that equation. We replace the ensemble average 
by an energy average and use contour integration. That 
is possible for F <C d because the resonances are iso- 
lated. We define F^^ = 2ttW^^ and F^ = EcTmc- The 
partial widths F^a sltc squares of Gaussian-distributed 
random variables and follow the Porter-Thomas distri- 
bution given in Eq. (1. 18). According to Eq. (1. 17) that 
distribution depends on a single parameter, the aver- 



age w idth (Tua)- We find (Bethe ^ l937 : Lane and Lvnnl . 



19571) and (iMoldaueil . 119611) to (jMoldaueH. \197^ 



27r / F^gF^fc 
~d 



(63) 



The average on the right-hand side is over the Porter- 
Thomas distribution of the F^gS. Since F^ = X^a^Mo, 
the variables in numerator and denominator are not sta- 
tistically independent. We use that for Xa ^ 1 we have 
Ta ~ 4:Xa = {2tt / d) {T f^a) ■ Thus, it is possible to express 
the right-hand side of Eq. (including the Porter- 
Thomas distribution of the F^a) completely in terms of 
the transmission coefficients Ta, as required by the sta- 
tistical theory. We write Eq. (|63|) in the form 



\S, 



(CN fl)|2\ 



(64) 



Here Waf, is referred to as the "width fluctuation correc- 
tion" (to the Hauser-Feshbach formula) and is defined 
as 



r^aF^fcV (F^ 



(65) 



Eq. shows that for the single-channel case we have 
W = 1. For the general case (several open channels) 
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Wn/) h as been calculated numerically by Lane and LvnnI 
(|1957[ ) and especially by iReffo et al\ ([19761). The val- 
ues range from one to three. This shows that for iso- 
lated resonances, the Hauser-Feshbach formula is only 
approximately valid. For the special case of a large 
number of channels with very small transmission coef- 
ficients each (so that A ^ 1 but X^c-^c ^ 1) numer- 
ator and denominator in the first term on the right- 
hand side of Eq. (|65p become uncorrelated, we have 



« (r^ar^6>/(r^a)(r,.f,) = l + 2(5a&, and the Hauser- 
Feshbach formula does apply with an elastic enhance- 
ment factor of three. That is why it is sometimes stated 
that the elastic enhancement factor increases from the 
value 2 for r ^ d (see Section IV.Bj) to the value 3 for 
r <^ d although the last value actually applies only in a 
special situation. 

Eq. holds for completely isolated resonances. A 
correction taking into account weak re sonance overlap 
and u si ng fl-matrix theo ry was given bv lLane and LvnnI 
(|1957[ ). iMoldaued (|l96ll) went beyond Eq. by us- 
ing a perturbative expansion of the S'-matrix in powers 
of the non-diagonal elements of the width matrix. The 
expressions become soon very cumbers ome. Therefore, 
iMoldaueil (fl96ll Il963l [l96l I1969L [l975h (see also further 
references therein) tried to go beyond the limit T <^ dhy 
usi ng the pole expansion o f the S '-matrix as formulated 
bv lHumblet and Rosenfeldl ( 196l[) . But the pole param- 
eters (locations of poles and values of the residues) are 
linked by unitarity in a complicated way and, thus, not 
statistically independent. That difficulty has never been 
resolved. We return to some aspects of Moldauer's work 
in Section EE] 

Eq. predicts average cross sections for elastic and 
inelastic scattering as functions of the transmission co- 
efficients. It would be of interest to calculate pertur- 
batively also the cross-section autocorrelation function 
^|^(CN fl)(£;^)|2|^(CN ^)(^E2)\^). That seems not to have 
been done yet. Likewise, the quantity (5™ ^{S^^^ ")*) 
seems not to have been worked out explicitly for T <^ d. 
As remarked in the introduction to Section|Vl that quan- 
tity is needed for the average CN cross section in the 
presence of direct reactions. The reason for the neglect 
is probably that F <C d is strictly realized only for a single 
open channel. 



B. Ericson Regime (F :» d) 



The Ericson model reviewed in Section FlLDI leads to in- 
teresting predictions that agree with experiment but has 
the shortcomings listed at the beginning of Section IIIII 
While it is probably not possible to derive the Ericson 
model as such from the RMT approach developed in 
Section UTT] we now show that it is indeed possible to 
derive both, the Hauser-Feshbach formula and all re- 
sults of the Ericson model, from that approach. That 
is done with the help of an asymptotic expansion in 
powers of d/T. Two methods have been used, a dia- 



gramr natic approach (lAgassi et all Il975l) and the replica 
trick (|Weidenmulleij . Il984[) . These are reviewed in turn. 
Combining the results of both, one obtains a complete 
theory of CN reactions in the regime T ^ d. 



1. Diagrammatic Expansion 

We use the form (|49l50p of 5"^^ and the statistical 
properties that come with that representation. We note 



that with F^^{E) 



■ J2c W^cWc^, we have 



(66) 



The last equation defines /. With (5F^y = F^^, — 6^^f, 
the matrix D is written as = (E — — if)S^i, — 

SFgi, = b^^^5^y — SFfj^iy. In the approach of lAgassi et all 
(I1975D . each of the two S'-matrix elements in the product 
Sab'iEi)iScdiE2))* is expanded in a Born series with re- 
spect to SFf,^. To calculate {S^i^ {Ei){S^]^ {E2)y} , one 
first performs the ensemble average over the Ws. That 
average is taken separately for each term of the double 
Born scries. The matrix SF^^ is bilinear in the Gaussian- 
distributed Ws. Each term of the double Born series 
therefore contains a product of Ws. The ensemble av- 
erage over such a product is taken by Wick contraction: 
The Ws are grouped in pairs, each pair is replaced by 
its ensemble average, see Eq. pS)) . All possible ways of 
pairing the Ws must be taken into account, including 
pairs where one member stems from the Born series of 
S'^^ and the other, from that of (S™)*. The number of 
ways of pairing the Ws increases dramatically with the 
order of the double Born series. Among these, only those 
which contribute to leading order in d/T, are kept. These 
are identified with the help of two rules, (i) Contrac- 
tion patterns are neglected that yield terms of the form 

Eu(-Di°V'((-Dl°^)*)"' withni > 2and/orn2 > 2. That 
is because when the summation over fi is changed into an 
energy integration and the integrals are done via contour 
integration, the result vanishes, (ii) The second rule is il- 
lustrated by the following example. We compare two sim- 
ple contraction patterns both of which occur as parts in 
the Born series and yield a term proportional to (d')^^)^. 

The first one occurs in the term D''°''[{WiWD''°^)'^]f,^ 
and involves the contraction of the first with the third 
and of the second with the fourth factor W. From 
Eq. (gHD the resuh is Sf^piD^'^'^f J2ai^l)^ ■ The second 
pattern occurs in the term Dl^\{WiWD^'^^)^]p_i, and in- 
volves the contraction of the first with the fourth, of the 
second with the third, of the fifth with the eighth, and 
of the sixth with the seventh factor W. The result is 
S^^AD'i?^)Hi-^^^/d)Eaivl)Y ■ The first (second) pat- 
tern yields a single (double) sum over channels. Thus, 
the first pattern is small of order d/F compared to the 
second and is neglected. Denoting the contraction of a 
pair of Ws with an overbar, one finds that only those 
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contractions survive in leading order in d/V where no 
two contraction lines intersect (nested contributions). 

In evaluating these patterns we have replaced the sum- 
mation over eigenvalues E"^ by an integration over en- 
ergy. That yields the factors l/d in the second pattern. 
Doing so corresponds to the neglect of correlations be- 
tween eigenvalues (the actual eigenvalue distribution is 
replaced by one with constant spacings). Wigner-Dyson 
eigenvalue correlations have a typical range given by d, 
and for d ^ r one expects that the neglect is justified. 
With the help of that same approximation the remaining 
terms in the double Born series can be resummed. With 
e = E2 — El that yields 



Tan 



(67) 



For £ = that coincides with the Hauscr-Fcshbach for- 
mula ^ and for e ^ agrees with Ericson's predic- 
tion if we identify in the latter the factors Air'^v'^/d 
with the transmission coefficients Ta , see the text follow- 
ing Eq. Moreover, in the framework of the diagram- 
matic approach, the Weisskopf estimate (fTTj) is seen to 
yield the exact expression for the width of the S-matrix 
correlation function in the Ericson regime. Eq. (j67p im- 
plies that the elastic (a = b) CN cross section is enhanced 
over the inelastic one by a factor of two ( "elasic enhance- 
ment factor" ) . This result is beautifully confir med exper- 
imentally, see Fig. 121 In ( Agassi et al\ . \l97di ) the terms 
of next order in d/T were also calculated. We return to 
that point in Section rV.B.2l 

The diagrammatic expansion can also be used to cal- 
culate higher-order ^'-ma trix correlation functions. It is 
found ( Agassi et aLl . ll975[ ) that to leading order in d/T, 
such correlations can be expressed completely in terms 
of correlations of pairs of S'-matrix elements. That im- 
plies that in the Ericson regime, the S'-matrix elements 
possess a Gaussian distribution. 

In summary, the diagrammatic approach shows that 
for T ^ d, the S'-matrix elements arc Gaussian- 
distributed random processes, with first and second mo- 
ments given by Eqs. (j56p and (|67p . respectively, and with 
a correlation width given by the Weisskopf estimate ([TT|) . 
All of this agrees with the predictions of the Ericson 
model. We conclude that within the diagrammatic ap- 
proach (i.e., under the neglect of eigenvalue correlations), 
the distribution of the S-matrix elements is completely 
known and is determined in terms of the average S'- 
matrix elements (j56p via the transmission coefficients ([5]) . 



2. Replica Trick 

Calculation of the average of |S^^^p is difficult be- 
cause i^GiO^ appears in the denominator of S'~^^ . A 
way to overcome that p roblem consists in th e use of a 
generating functional Z (jWeidenmiilleil . Il984l) . The ob- 
servable (here: IS^jf'p) is given in terms of a suitable 



derivative of Z, and iJ^^E appears in Z as the argu- 
ment of an exponential. That fact would greatly simplify 
the integration over the random variables (the matrix 
elements of ff°°^) were it not for the need to normal- 
ize Z. The normalization problem is overcome by the 
replica trick originally developed in condensed-matter 
physics ( Edwards and Andersonl . [T975I ). 

We introduce N real integration variables tpf^, fi — 
I, . . . , N and define the generating functional 

ZiE, J) ^[lif ^^m) exp I (z/2) ^,E4'^ } 

X exp I (i/2) ^ V^(-i?M<> + «7r^ W'^cVKci.)^'.. [ • 

^ fil/ C ^ 

(68) 

Here the matrix elements W include the "source terms" 
J and are defined as 



The CN scattering matrix is given by 



d 



(69) 



(70) 



j=o 



Unfortunately, little has been gained because calculating 
{Sab) or (•S™(5,T)*> from Eq. dTO]) involves averaging 
the logarithmic derivative of Z or the product of two 
logarithmic derivatives of Z and is next to impossible to 
do. The difficulty is that at J = 0, Z is not and cannot 
easily be normalized to unity. The problem is overcome 
by using the identity 



InZ = lim -(Z" - 1) 

ri-!-0 n 



(71) 



For integer values of n, the average of (or of a prod- 
uct of such terms) can be calculated, and Eq. ((7T|) is then 
used to calculate (In Z) or of ((In Z)-^). Instead of In Z we 
use n replicas of Z to calculate averages, hence the name 
of the method. The calculation must be done analyti- 
cally as otherwise the limit n — >■ cannot be taken. That 
is possible only for positive integer values of n. However, 
Eq. ([7T|) is strictly valid only when n is not restricted to 
integer values. Otherwise one may miss non-zero contri- 
butions that happen to vanish for all positive integer n. 
That is why the method is not guaranteed to be exact 
and is referred to as a "trick" . Later investigations have 
shown, however, that when used for an asymptotic ex- 
pansion as is done below, the replica trick gives correct 
answers ( Verbaarschot and Zirnbaueij . [l985l) . 

A detailed d es cripti on of the calcula- 
tion (IWeidenmliUed . Il984l) exceeds the frame of 
this review. We only sketch the essential steps. In 
the case of (Z"), we deal with n x N real integration 
variables i/'^ with k = 1, . . . ,n and fj, = 1, . . . ,N. With 
the help of Eq. (|36p . averaging over i/GOE yjci^jg ^ 
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quartic term in the f/'^s. That term can be written as 
the trace of the square of the symmetric real n x n 
matrix Akk' — J2fj. V'^V'^i • The form of A reflects the 
orthogonal invariance of the GOE. The bilinear term 
in A is removed by means of a Hubbard-Stratonovich 
transformation, 



exp 



n 

k<k' 



^Trace (A^) 



+ CJO 



dcTfcfc' exp 



—Trace (a^) 
4 



X exp 



iATracc((TA) 



(72) 



In Eq. ((72|) we have introduced a set of new integra- 
tion variables which appear in the form of the sym- 
metric real n-dimensional matrix Ok^i . After removal 
of the quartic term, the remaining integral over the 
V'^s is Gaussian and can be done. The integral over 
CTfcfc' is performed with the help of the saddle-point ap- 
proximation. For iV 3> 1 that approximation is excel- 
lent. For {S^^), it yields an asymptotic expansion in 
inverse powers of N with the right-hand side of Eq. ([56|) 
as the leading term. Proceeding in the same way for 

(5^™ ^\Ei){sf^ ^\E2))*) one finds instead of a sin- 
gle saddle point a continuum of s addle points. In a sem- 
inal paper, ISchafer and Wegneil (ll980D have shown in a 
different context how to deal with that manifold. That 
is done with the help of a suitable parametrization of 
the matrix a^k' which displays both, the massive inte- 
gration variables (that describe integration points out- 
side the saddle-point manifold) and the Goldstone mode 
(which defines the coordinates within the saddle-point 
manifold). After these steps, one finds that the energy 
difference E2 — Ei = e appears in the integrand only in 
the combination X^c-^c + {2iTTe/d). That shows that the 
correlation width F of the correlation function is given 
by the Weisskopf estimate pT|) . The integrals over the 
saddle-point manifold cannot be done exactly. For the 

calculation of (5^™ "^(£^i)(5^d^ "^(-^2))*), the replica 
trick can only be used to calculate an asymptotic expan- 
sion in inverse powers of Tc- The leading term agrees 
with Eq. (|67l). T erms of higher orde r also agree with 
the result given bv lAgassi et al\ ( 19751 ). Moreover, these 
terms can be used to confirm unitarity for each order of 
the expansion and to show that the leading term gives re- 
liable answers for ^^Tc > 10 or so, see also Section lV.CI 
In the case of direct reactions, we use Eq. PO)) and 
obtain from Eq. (p7)) . the definition P7|) . and with Ta = 
Pa as a result Eq. (ITCl) . That shows tha t to leading order 
in d/T, the result of lKawai et all (|l973l ) agrees with that 
of the statistical theory. 



3. Summary 

The replica trick yields an asymptotic expansion in in- 
verse powers of ^ J, Tc for the S'-matrix correlation func- 



tion. The leading term is given by Eq. ((67|) . That form 
implies that the correlation width agrees exactly with the 
Weisskopf estimate ([TT|). The diagrammatic approach 
neglects eigenvalue correlations of the GOE Hamiltonian 
in favor of a constant-spacing model but leads likewise 
to Eq. (p7|) . That shows that for F ^ d such correla- 
tions are indeed negligible, and that the results of the 
diagrammatic approach are trustworthy. Investigating 
higher-order correlation functions within the diagram- 
matic approach (which seems prohibitively difficulty in 
the framework of the replica trick) , one finds that the S'- 
matrix elements are Gaussian-distributed random vari- 
ables with first and second moments given by Eqs. (j56p 
and (p7|) . respectively. We see that combining the replica 
trick and the diagrammatic approach we obtain a com- 
plete theoretical understanding of the distribution of the 
elements of S''-^^ in the Ericson regime. The results agree 
with predictions of the Ericson mo del and, in th e pres- 
ence of direct reactions, with those of lKawai aLl(tl973i') . 



C. S-Matrix Correlation Function 

The replica trick yields only an asymptotic expansion 
but not the full S'-matrix correlation function. That 
functi on caii be obtai ned exactly using "supersy mmc- 
try" (lEfetovl . [1983; Vcrbaarschot and Zirnbaueii [l985D . 
We briefly motivate the use of that method. 

The use of the generating functional defined in Eq. 
would greatly simplify if it were possible to normalize Z 
so that Z{E,0) = 1. Then in Eq. ^ we would have 
[d\nZ{E,J)/dJ]\j=o = [dZ{E,J)/dJ]\j^o- Averaging 
S would amount to averaging Z (and not In Z) , and it 
would not be necessary to use the replica trick. Since 
in Z the random variables appear in the exponent, the 
calculation of (Z) would be straightforward. 

That goal is achieved by defining the normalized gen- 
erating functional as the product of two factors. The 
first factor is Z^ with Z as defined in Eq. §8^. The 
second factor has the same form as Z except that the in- 
tegration variables anticommute. Integration with anti- 
comm uting variables i s a well-known mathematical tech- 
nique ( BerezirJ . Il986l ). In the present context, we use 
the following property. The normalization factor of a 
Gaussian integral involving N commuting real integra- 
tion variables ipf^, fj, = 1, . . . , N and a symmetric matrix 
A is given by 



+00 



Yl / dtp^, exp <^ {i/2) Vm^m^V'iv 



[det(A/(2i7r))] 



-1/2 



(73) 



The same integral with the commuting variables ip^ re- 
placed by anticommuting variables ^'^d x*^ is given 
by 



[det(A/(2i7r))] . 

(74) 
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Combining two factors of the form of the left-hand side 
of Eq. ([75]) with one factor of the form of the left- 
hand side of Eq. (f74|) we obtain a normalized generat- 
ing functional involving both commuting and anticom- 
muting integration variables. The term "supersymme- 
try" commonly used in that context is somewhat inap- 
prop riate. It s tems from rclativistic quantum field the- 
ory ( Wess and Zum ino. 1974a.bj) where fcrmionic (i.e., 
anticommuting) and bosonic (i.e., commuting) fields are 
connected by a "supersymmetry" . That specific symme- 
try does not occur in the present context. 

In a formal sense, the calculation of {S™) and of the 
correlation function runs in parallel to the one for the 
replica trick. In content, the steps diff'er because of the 
simultaneous use of commuting and anticommuting vari- 
ables. The steps arc: Averaging the generating functional 
over the GOE, replacing the quartic terms in the inte- 
gration variables generated that way with the help of a 
Hubbard-Stratonovich transformation, execution of the 
remaining Gaussian integrals over the original integration 
variables, use of the saddle-point approximation for the 
remaining integrals over the supcrmatrix a. That ma- 
trix is introduced via the Hubbard-Stratonovich trans- 
formation and has both commuting and anticommuting 
elements. For ^ 1 the saddle-point approximation 
is excellent and yields for {S'~'^) an asymptotic expan- 
sion in inverse powers of N , the leading term being given 
by the right-hand side of Eq. ([55)1 when E is taken in 
the center of the GOE spectrum. When the same for- 
malism is used for the S'-matrix correlation function in- 
stead of S itself, the saddle-point changes into a saddle- 
point manifold. With the help of a suitable paramctriza- 
tion (jSchafer and Wegneii Il980l ) of (J, the integration 
over that manifold can be done exactly. The result is 
again valid t o leading order in (\/N). W ith e = E2 — Ei, 
one obtains ( Verbaarschot et aLl . ll985l ) 



Jo Jo 

1 r ITTS 

X- ^(Ai,A2,A) exp j - _(Ai + A2 + 2A) 



(1 - TeA) 



(i-KreAi)i/2(i + r^A2)i/2 

X>/af)cd(Al, A2, A) . 



(75) 



The factor /i(Ai,A2,A) is an integration measure and 
given by 



m(Ai,A2, A) = 



(l-A)A|Ai -A2I 



nti[((l + A0A0^/^(A + A.)2] 



(76) 



while 



>/a6cd(Al, A2, A) — {Sac^bd + SadSbc)TaTi, 



A,(l + A.) 



2A(1 - A) 



- (l + T,AO(l + TfcA,) (l-r,A)(l-rbA) 



X, 2A 



2 



TaX,) 1 - TaX 

A,- 2A 



TA) 1-T,A 



(77) 



describes the dependence of the correlation function on 
those channels which appear explicitly on the left-hand 
side of Eq. ([75l). 

Eqs. ([75|) to ([77]) give the S'-matrix correlation func- 
tion in closed form, i.e., in terms of an integral represen- 
tation. It does not seem possible to perform the remain- 
ing integrations analytically for an arbitrary number of 
channels and for arbitrary values of the transmission co- 
efficients. In the way it is written, Eq. (|75p is not suited 
very well for a numerical evaluation because there seem to 
be singularities as the integration variables tend to zero. 
Moreover, the exponential function oscillates strongly for 
large values of the As. These difficulties are overcome by 
choosing another set of i ntegr ation variables. Details are 
given in ( Verbaarschoti Il986[ ). For the case of unitary 
symmetry, for mulas corresponding t o Eqs. ([75]) to ([77)1 
were given m (|Fvodorov et al\ . \200^ . 

We turn to the physical content of Eqs. (|75p to ()77p . 
The unitarity condition ([l} must also hold after averaging 
for S™ . Using a Ward identity one finds t hat unitarity 
is indeed obeyed ( Verbaarschot et all . Il985l ) . Except for 
the overall phase factors of the average S-matriccs ap- 
pearing on the right-hand side of Eq. (|77p . the S-matrix 
correlation function depends only on the transmission co- 
efficients, as expected. Eqs. ([75]) to ([77)) apply over the 
entire GOE spectrum if d is taken to be the average GOE 
level spacing at E = {l/2){Ei + E2). That stationarity 
property enhances confidence in the result. When E is 
chosen in the center of the GOE spectrum, E = 0, then 
(S™) is real, see Eq. (|56p . and the complex conjugate 
sign in the last term in Eq. ([77]) is redundant. Writing 
the product over channels in Eq. (|75p as the exponential 
of a logarithm, expanding the latter in powers of Ai, A2, 
and A, and collecting terms, one finds that e appears only 
in the combination [{2iTTe/d) + J2c -^c]- This fact was also 
established in the framework of the replica trick. It shows 
that F as defined by the Weisskopf estimate ([TT|) defines 
the scale for the dependence of the correlation function 
on s. Put differently, the energy difference e appears uni- 
versally in the dimensionless form e/F (and not e/d), and 
for {a, b} = {c, d} the correlation function ([75)1 has the 
form /(I + ie/T; Tq, T^; Ti , T2, . . . , Ta). The right-hand 
side of Eq. ([77)) is the sum of two terms. These cor- 



respond, respectively, to {si'','' ''\e,) (5^^ "'(^2))*) 

and to {Sia^ ^\Ei) {si^^ ^\E2))*) and are exactly the 
terms expected, see the introduction to Section [V) 

The limiting cases (F < d and F » d) of Eqs. ([75)) 
to (I7 7J) were studied and compared with previous re- 
sults ( Verbaarschoti Il986[) . The case F ^ d is obtained 
by expanding the result in Eqs. ()75p to ([77)) in powers 
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of the transmission coefficients. The result agrees with 
Eq. Conversely, expanding the result in Eqs. ([75]) to 
([77]) in inverse powers of one generates the same se- 

ries as obtained from the replica trick, see Section lV.B.2l 
That shows that the Hauser-Feshbach formula (|9]), the 
extension ([T6|) of that formula to the case of direct re- 
actions, and the Ericson result (|67| all are the leading 
terms in an asymptotic expansion of the relevant expres- 
sions in powers of d/V. Moreover, in the limit F ^ d 
the Weisskopf estimate pT|) gives the exact expression 
for the correlation width. All of these facts were also ob- 
tained previously with the help of the replica trick, see 
Section rV.B. 21 A more detailed investigatio n using su- 
persymmetry ( Davis and Boosel . Il988l [l989( ) has shown 
that the elastic S'-matrix elements Saa possess the Gaus- 
sian distribution assumed by Ericson only for strong ab- 
sorption [Ta ~ 1). If that condition is not met, devi- 
ations from the Gaussian are caused by unitarity, see 
Section lEB 

Eqs. ((75|) to ([77)) contain the central result of the statis- 
tical theory. The S'-matrix correlation function is given 
analytically for all values of F/d. These equations su- 
persede both the perturbative approach (Section IV. A[) 
and the replica trick (Section IV.B.2P because the results 
of both these approaches turn out to be special cases 
of the general result. It would be highly desirable to 
extend the supersymmetry approach to the calculation 
of moments and correlation functions involving higher 
powers of S^"^ and {S'^^)* than the first. For the cor- 
relation function (|77|) . the supermatrix a has dimension 
eight. For the cross-section correlation function, that 
matrix has dimension sixteen. Integration over all ma- 
trix elements must be done analytically. While that is 
feasible for the function ([77| . it seems beyond reach for 
the cross-section correlation. That is why the results of 
the diagrammatic approach (Section IV.B.ip are needed 
to complete the theory in the Ericson regime: They show 
that the S'-matrix elements have a Gaussian distribution. 
There is no analytical information about cross-section 
fluctuations versus energy outside the Ericson regime. 

We note that in both limits F <g; c? and F ^ d the fluc- 
tuation properties of the S'-matrix are completely deter- 
mined by the Gaussian distribution of the eigenvectors of 
£jGOE g^jj^ infiepgndent of the fluctuation properties 
of the eigenvalues of that matrix. Indeed, for F ^ d the 
value of (ISfjJ^^p) in Eqs. ([Mj) and ([HS]) is determined en- 
tirely by the Porter-Thomas distribution of the partial 
widths (which in turn is a consequence of the Gaussian 
distribution of the eigenvectors of 77^*^^). And for F ^ d 
(Eric son limit) the asymptotic expansion ([Agassi et all . 
I1975D based on a picket-fence model for the eigenvalues 
gives the same result as the calculation that takes fully 
into account the eigenvalue correlations of _ff*^°^. A sig- 
nificant dependence of the S-matrix fluctuations on the 
distribution of the eigenvalues of i^ci^E ^^^^^ thus, occur 
only in the intermediate domain F « rf. 

A dynamical model for an interacting fermionic many- 
body system coupled to a number of open channlcs 



was investigated by ICelardo et al. I ([20071 iooi). Spin- 
less fermions are distributed over a number of single- 
particle states with Poissonian level statistics and inter- 
act via a two-body interaction of the EG0E(2) type, see 
Section I.V.A. The Hilbert space is spanned by Slater 
determinants labeled ji ~ 1,...,N. These are cou- 
pled to the open channels labeled a,b, . . . by amplitudes 
A" which are assumed to have a Gaussian distribution 
with mean value zero and a second moment given by 
{A^j^Al) = SabSfj.uY'/N. The parameters 7" determine 
the strength of the coupling to the channels. One focus 
of these papers is on the way the S-matrix fiuctuations 
change when the intrinsic dynamics is changed from reg- 
ular (Poisson statistics) to chaotic (Wigner-Dyson statis- 
tics) by increasing the strength of the two-body interac- 
tion. Many of the features discussed in t his review are 
i llustrated b y the numerical simulations bv lCelardo et al\ 
([2007[ . [20081 ). The assumption that the amplitudes 
have a Gaussian distribution with a second moment that 
is proportional to the unit matrix in the space of Slater 
determinants, puts a constraint on the calculations as 
it implies that the distribution of the is invariant 
under orthogonal transformations in the space of Slater 
determinants. Put differently, together with the intrin- 
sic Hamiltonian Hfj^i, all Hamiltonian matrices obtained 
from Hf^i, by orthogonal transformations yield the same 
distribution for the S-matrix elements. This induced or- 
thogonal invariance of H implies that the eigenvectors 
of H arc Gaussian-distributed random variables, com- 
pletely independently of the detailed form of H, and that 
only the distribution of the eigenvalues of H depends on 
whether the intrinsic dynamics is regular or chaotic. The 
statistical assumption on the amplitudes A° is physically 
justified when many channels are open and when the 
channel wave functions describe states that are chaotic 
themselves. But the assumption allows for a partial test 
only of the transition from regular to chaotic motion. 



Eqs. ([75[) t o (1771) have been derived in the limit — > 00 
and A fixed. [Lehmann et all ( 1995bf ) pointed out that in 
the extreme Ericso n regime y^„Te ^ 1, a nother limit 
(first considered bv [Sokolov and Zelevinskvl (fl988. 198^ 
Il992[ )') may be more appropriate: A and N tend jointly 
to 00 while the ratio m = A/N is kept fixed. The au- 
thors calculated the S-matrix correlation function in that 
limit. Using supersymmetry. they find corrections of or- 
der m to the saddle-point equation. These modify the 
S'-matrix correlation function but keep the correlation 
width (given by the Weisskopf estimate) unchanged. The 
modifications are due to the fact that in the limit con- 
sidered, the range of the correlation function becomes 
comparable to the range of the GOE spectrum, and the 
universality that characterizes the regime F <C A (where 
2 A is the radius of the GOE semicircle) is lost. Further 
details are given in Section [V.FI 



28 



Decay in Time of the Compound Nucleus 

Rather than studying the dependence of the correla- 
tion function (|75p on the energy difference e, we inves- 
tigate the Fourier transform of that function. We con- 
fine ourselves to pure CN scattering (so that the ma- 
trix U in Eq. (|40)) is the unit matrix). We first show 
that for a very short wave packet incident in channel a, 
the time dependence of the CN decay feeding channel 
b is given by the Fourier t r ansfo r m of the corr e lation 
function (1751) (iDittes et all . Il992l: iHarnev et all . Il992l: 
iLvuboshitd . 1978al lb(). We then use the explicit form of 
that function to determine the time dependence of CN 
decay. 

Let / dE cxp{-iEt) g{E) Yia{E)) with / AE \g{E)\'^ 
= 1 be a normalized wave packet incident in channel a. 
The functions \^ a{E)) are solutions of the Schrddinger 
equation for a CN system with A channels and ^ 1 
quasibound states subject to the boundary condition that 
there is an incident wave with unit flux in channel a only. 
We put h = \. The wave packet is short in time if g{E) is 
very broad, i.e., covers many resonances. The outgoing 
flux in any channel h is asymptotically (large distance) 
given by ^ dEi ] dE^ exp[i(£;2 - Ei)t] g{Ei)g*{E2) 
^ab (-^1 ) i^ab (-^2 ) ) * ■ We introduce new integration vari- 
ables E = {l/2){Ei + E2) Si-ad e = E2 - Ei and ob- 
tain J dE ] de exp[iei] g{E - {l/2)e)g*{E -f- (l/2)e) 
SabiE-{l/2)e){S^^{E + {l/2)e))*. Under the assump- 
tion that g{E) changes very slowly over a scale of order d 
or r (whichever is bigger), we may write that expression 
as / de eM^et] J dE \g{E)\' S^^^ {E - {l/2)e) {S^i^{E + 
{l/2)e)r = Jds exp[*ei] {Sg^ {E - {1/2) s){S^i' {E + 
(l/2)£))*). We use the decomposition ^ and obtain a 
sum of two contributions. The first, 6{t)\{S^^)\'^ , shows 
that the average 5'-matrix describes the fast part of the 
reaction as claimed in Section IIVI The second has the 
form 

Pabit) = |deexpM(5i;;N«)(i?-(l/2)£) 



{l/2)e)r 



(78) 



That shows that the Fourier transform (FT) of the cor- 
relation function ([75]) gives the time dependence of the 
CN decay. More precisely, the FT is the derivative of 
the function which describes the decay in ti me of the CN 
in channel b within the statistical theorv ([Pittes et all . 
I1992I: iHarnev et all Il992l: ILvuboshitd . ll978al lbD The 
form of t he matrix D i mplies that Pab(t) = for t < 
as expected ( Dittes et all . Il992l : IHarnev et aU[l993) . 

The function Pab{t) cannot be worked out analytically 
in i ts full generality. Some limiting cases are of inter- 
est (jPittes ami992l:lHa7nev et ai.l . ll992l ). For a single 
open channel, the decay is not exponential. Rather, it 
is asymptotically (t — 00) given by the power-law de- 
pendence t~^/'^. For A channels weakly coupled to the 
resonances (all Tq <C 1), the decay has asymptotically the 
form t~^~^/'^. In both cases, the characteristic mean de- 
cay time is given by the average of the coupling strengths 



of the resonances to the channels. Exponential decay 
given by exp[— Ft] with F given by the Weisskopf esti- 
mate pT|) is realized only in the Ericson regime F 3> d. 
The deviations from the exponential decay law are due 
to the Porter-Thomas distribution of the partial widths, 
see Section LILD.l. The decay in time of an isolated 
resonance is exponential, of course, except for very small 
and very large times. But superposing many such reso- 
nances with different widths causes deviations from the 
exponential distribution. Such deviations ha ve been ob- 
serve d experimentally, see Section [VLAi In ( Hart et'oH . 
I2009D . a somewhat different conclusion was drawn. The 
decay in time of excitations in a chaotic microwave cav- 
ity was investigated. The cavity mimicks a chaotic quan- 
tum system. Therefore, the general RMT results should 
apply. The authors find that in the limit t 00 , 
the ensemble avera ge considered in ( Dittes et all . Il992l : 
IHarnev et all . Il992l ) differs from the behavior of (every) 
single realization. For the latter, the longest-living res- 
onances finally dominate, these are well separated in en- 
ergy, and the decay, therefore, becomes eventually ex- 
ponential. The transition from algebraic to exponential 
deca y follows a univers al law if time is properly normal- 
ized (iHart et a/.l. l2009l) . T he apparent discrepancy w ith 



ipp a 

the result of (iDittes et all . Il992t iHarnev et all. I1992D is 



due (fOttl . I2OO9D to a finite-size effect: iHart et all (|2009l ) 
considered the case o f a finite number N of resonances , 
while the result of ( Dittes et all . 119921 IHarnev et all . 
I1992D holds in the limit A^ — > 00. The turning point 
in time where in ( Hart et I2009D power-law decay 
changes into exponential decay increases with N and is 
moved to infinity for N 00. 

Another approach to tim e delay uses t he Wigner- 
Smith time-delay matrix ( Smithl . Il960allbl : IWigneil . 
ll955bD 



Qab 



de 



xisi':'' «)(i? + (l/2)e))*) 



(79) 



r=0 



That matrix is the matrix of average time delays. In- 
deed, it is easy to see that Qaa = (27r)~^ / dt t Pab{t)- 
The eigenvalues of the Hermitean matrix Q are called 
the proper delay times. For A channels with Tc = 
1 for all c, the joint probability distribution and the 
dens ity of the proper delay t i mes have b e en w orked 
out (|Brouwer et a/.l . ll995l . [19971 ). ISavin et all i^200lh gen- 
eralized that approach to A equivalent channels with 
Tc ^ 1. Analytical formulas valid for an arbi- 
trary number of channels and arbitra ry values of the 
transm ission coefficients are given in ( Lehmann et all . 
Il995bl ). For time-reversal non-invariant s ystems 
these issues were treated iii ( Fvodorov et all . Il997bl : 
iFvodorov and Sommersl. 1997 ). The single-c hannel case 
was studied in ( Ossipov and Fvodorovl l2005l) . 
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D. Distribution of S-IVIatrix Elements 

Except for the Ericson regime, correlation functions 
that relate to physical obscrvables and involve higher 
powers of S than in Eq. (|75p . are not known. Is the 
situation better when we ask for the distribution of S- 
matrix elements all taken at the same energy? As we 
shall see, the answer lies strangely between yes and no. 



1. Fit Formulas 

Fit formulas for the second moments of 5* based 
upon an RMT simulation were developed prior to the 
derivation of the general result in Eqs. (|75p to ((77|) 
and, in principle, have been superseded by that de- 
velopment. We give these formulas here b ecause they 
are st ill fr equently used i n app lications. iTepel et al\ 
(|1974[ ) and lHofman~ et al. 1 (Il975ai rbl) used an ansatz for 

(I'^ab^^ ^■'P) that was inspired by the Hauser-Feshbach 
formula. It reads 



{\S™?) 



VgVb 



(1 + (W, - l)5ab) 



(80) 



Unitarity relates the expressions on the right-hand side 
of Eq. (|80|) with the transmission coefficients Ta- The 
resulting equations possess unique solutions for the VaS 
provided the W^s are known. That leaves the latter as 
the only parameters to be determined by fits to a nu- 
merical simulation. For numerous sets of transmission 
coefficients, the simulation was done using the if-matrix 
form dSni) of S'^'^, taking in Eq. ^ the E^,s as eigen- 
values and determining the wi^^s in terms of the eigen- 
functions of a GOE matrix, with a strength defined by 
Ta- The result is a fit formula for Wa, 



W„ - 1 = 



0.3+1.5(T„/^^Te; 



1 + {Ta) 



Here T is the arithmetic mean of the T^s. 



(81) 

The other non- 
vanishing bilinear form {Sj^i^ ^\^bb^ *'))*^ -^ag simi- 
larly assumed to factorize for a ^ b. Fit formu- 
las for th at function based on si mulations are likewise 
given m (iHofmann et al\ . Il975bt ). Numerical evalua- 
tion (jVerbaarschotl Il986l ) of Eqs. ([75]) to ([771) showed 
good agreement with Eqs. (|80)) and (|8T|) within the statis- 



tical errors ex pected. Similar f ormulas were subsequently 



ticai errors ex pectea. ai r unar l o 
developed bv lMoldaueil (|l976[ ). 



By definition, these have average value zero and, thus, 
transmission coefficients Tq = 1 in all channels. More- 
over, the limit N oo implies the (unrealistic) limit 
of infinitely many channels. It is perhaps for these rea- 
sons that the orthogonal circular ensemble has not found 
notable applications in nuclear physics. 

The distribution of 5-matrix elements is analytically 
accessible for systems with absorption. Absorption oc- 
curs, for instance, in microwave resonators where it is 
due to Ohmic losses. Absorption is described by intro- 
ducing ficticious channels and associated transmission co- 
efficients. If the resulting total widths of the levels are 
dominated by absorption, they become statistically inde- 
pendent of the partial width amplitudes of the physical 
channels. The calculation of the di stribution of 5'-matrix 
elements is then m uch simplified ( Fvodorov et a^.l . [20051 : 
ISavin et flZ.l . l200fih . 

The approach may be used to obtain partial informa- 
tion on the distribution of S'-matrix elements in CN reac- 
tions. One must focus attention on a distinct pair (a, b) of 
channels and assume that many channels are open (these 
may all have small transmission coefficients so that we do 
not necessarily work in the Ericson regime). All chan- 
nels different from (a, b) then play the same role as the 
ficticious channel s in the case of absorp tion, and the re- 
sults o btained bv lFvodorov et ail ( 20051 ) and lSavin et al\ 
( 20061 ) may be used to determine the distribution of Sab- 
The method obviously does not yield the joint distribu- 
tion function of all elements of the scattering matrix. We 
are not aware of applications of that approach to CN re- 
actions and do not reproduce the relevant formulas here, 
see, however. Section IV. El 



3. Exact Results for Low Moments 

The supersymmetry approach can be used to calcu- 
late the third and fourth moments of the scattering ma- 
trix (jPavis and Boosel Il988l Il989l ). The need to intro- 
duce supermatrices of dimension larger than eight (un- 
avoidable if one wishes to calculate higher-order correla- 
tion functions) is circumvented by writing these moments 
as higher-order derivatives of the very same generating 
functional that is used to calculate the result (|75|). The 
resulting analytical formulas are valid for all values of 
r/d. For more details, see Section fV.EI In view of the 
complexity of the calculations, an extension of that ap- 
proach to higher moments than the fourth seems very 
difficult. 



2. Many Open Channels 



4. Maximum-Entropy Approach 



iDvsonI (|l962a[) has defined the orthogonal ensemble of 
unitary symmetric matrices S of dimension N (the "cir- 
cular orthogonal ensemble") and has studied the distri- 
bution of its eigenvalues for N ^ oo. The members of 
that ensemble may be interpreted as scattering matrices. 



The approach deve loped by the Mexican 
group (|Mello et all Il985t) is based on an appealing 
idea. If S'^^ is determined by an RMT approach (as 
done in Section IIII.B|) . then itself should be as 

random as is consistent with basic properties of that 
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matrix. These properties are unitarity, symmetry, and 
the property Eq. ([57)1 to which the authors refer as 
"analyticity-ergodicity" . In addition, it is required that 
i^aa) should have the value given by the ^-matrix S°^^ 
of the optical model. With these requirements used as 
constraints, expressions for the distribution F{S) of 
are then obtained from either an analytical approach or 
a variational principle. The results agree. In the latter 
case, the probability density for S*^^ is determined by 
maximizing the entropy — J F{S)\nF{S)dii{S) under 
the said constraints. Here fi{S) is the Haar measure 
for unitary and symmetric matrices, see Section I.II.B. 
With A the number of channels, 1a the unit matrix in A 
dimensions, and V a A-dependent normalization factor, 
the result is 



CN^ ^ 1 [det(lA-(^°P^)*g°P*)](^+i)/^ 
^ ^ V |det(lA- (5°pt)*5CN)|A+i 



(82) 



The function F is the most likely distribution function 
for S™ under the constraints mentioned. Do results cal- 
culated from F{S) agree with those based on the RMT 
approach in Section IlII. BP The answer is a (co nditional) 
yes. F irst, the distribution ((82)) can be derived (jBrouweii 
from the stochastic S'-matrix in Eq. (|4ip under 
the assumption that the Hamiltonian is a member of the 
"Lorcntzian ensemble" (rather than of the GOE). The 
Lorcntzian ensemble and the GOE have the same eigen- 
vector distribution and th e same level-cor relation func- 
tions in the large N limit ( Brouweri Il995[ ). It is, there- 
fore, extremely likely that the distribution ([5^ also holds 
for the stochastic S'-matrix in Eq. (|4T|) with the Hamilto- 
nian taken from the GOE. That view is supported further 
by the following facts. For strong absorption {{Sab) = 
for all {a, b}) Eq. §^ reduces to F{S™) = constant. In 
other words, the distribution of 5-matrix elements is de- 
termined entirely by the Haar measure. In that limit 
the ensemble (15^ agrees, the refore, with Dyson's cir- 
cular ensemble ( Dvsonl . Il962a[) . In the Ericson regime, 
Eq. (|82|) yields the Hauser-Feshb ach formula with an 
elastic enhancement factor of two ( Friedman and Mellol . 
Il985l ). Moreover, for A = 2 and e = Eq. ((5^ 
agrees (|Verbaarschotl Il986| ) with Eqs. ^ to Fur- 
ther support comes from resu lts for the unitary case 
([Fvodorov and Sommer [T997^ . 



All these facts make it seem highly probable that 
Eq. (|82p correctly describes the distribution of S'-matrix 
elements for chaotic scattering. Unfortunately, in the 
general case of several open channels the expression on 
the right-hand side of Eq. is so unwieldy that it has 
not been possible so far to evaluate it. That puts us into 
the strange situation that we do seem to know the distri- 
bution of S-matrix elements without being able to use it. 
It must also be remembered that the maximum-entropy 
approach docs not yield information on S-matrix corre- 
lation functions. 



E. Cross-Section Fluctuations 

For the analysis or prediction of cross-section fluc- 
tuations, the theoretical results reviewed so far in this 
Section do not suffice. While the third and fourth mo- 
ments o£jS|^j^E}_at_^xed_energy__E are known analyti- 
cally (|Davis and Boosel Il988l Il989l ). information on the 
corresponding correlation functions of S^ {E) doe s not ex - 
ist. The problem was addressed bv lDietz et dl (|2009c[) . 
The authors used the available analytical results and in- 
formation obtained numerically and/or experimentally 
from microwave billiards (sec Section IVI.Ap to investi- 
gate cross-section fluctuations, and to identify the range 
of parameters where predictions can safely be made. 

We use Eqs. (|4T]) to (gS]) for S^^. The cross-section 
autocorrelation function Cab(e) is defined as 



Cab{e) 



{\S^r{E + e/2)\'\S^^^{E-s/2)\') 



(83) 



Using Eq. we write S'^^{E) = {S™) + S'^™ ^\ We 
use the fact (see Eq. ([56)) ) that {S™) is real and obtain 



Cab{s)^2Sab^\{S^r)\'^[Cil\e)] 
+ (^Sf > 5R[(S(™ ^>{E + £/2)|S(?^ ^He - s/2)\')] 
+ (^aT) SftK^r «)*(i?-e/2)|Sr«)(i? + s/2)p)] 



+ \{\slT'\E + e/2r\Si:f;'->{E~sm 

1 2 > 

q(CN fl)|2 



|2| c(CN fl). 



(84) 



Here C^j^ (e) is the S-matrix autocorrelation function in 
Eq. (|75)) taken at c = a, d = b. Eq. (|M)) shows that 
there is a substantial difference between the elastic case 
(a = b) and the inelastic one {a b), caused by the fact 
that {S^'^) is diagonal. We address the inelastic case 
first. In the Ericson regime, S^"^^ is Gaussian, and 
the autocorrelation function of IS*-^^**-*!^ is. therefore, 
given by the square of the S-matrix correlation func- 
tion in Eq. (|75)) . However, that relation cannot be ex- 
pected to hold much outside the Ericson regime because 
we must expect the cross-section fluctuations (in units of 
the average cross section) to increase significantly as F/d 
dec reases. Thi s expec tation is quantitatively confirmed 
in (jPietz et al. I. l2009d) : The autocorrelation function of 
|^(CN fl)|2 least approximately, given by the square 

of the S-matrix correlation function in Eq. ()75p when- 
ever S^'^^ posseses a bivariate Gaussian distribution, 
and that is essentially the case when F > c? or so. For the 
elastic case (a — b) the situation is difficult even in the 
Ericson regime unless |(S^^)| <C 1. Indeed, whenever 

that constraint is violated, the distribution of S^q^ 
cannot be Gaussian: Combined with the decomposition 
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, the constraint < 1 im- 



(CN fi) 



to 



CCN _ /qCN\ I c(CN fl) 
'->aa — Waa / "r '^aa 

phed by unitarity forces the distribution of S'l 
be skewed. The distortion of the Gaussian distribution 
grows with decreasing T/d and is strongest when the cou- 
phng to the channels becomes very smaU. (Then S^a 
is dominated by the first term on the right-hand side 
of Eq. (|41]l .) The terms in the second and third hne of 
Eq. ((M)) vanish only if the distribution of is Gaus- 

sian so that the phase of si^^ is uniformly distributed 
in the in terval |0,27 r |. Th at condition is found to be 
violated ( Dietz et a/J . 120093 ) already when T/d < 4, and 
a full evaluation of all terms on the right-hand side of 
Eq. ([M)) is then necessary. 

This is possib l e with the help of the results of 
iDavis and Bo'^ (|l988l . I1989D who calculated analyti- 
cally the functions 



FrA^) = {[S^JiE + e/2)Y[S^,,{E-s/2)Y) , 
F!.l\e) = {S^JiE + e/2) [S!,{E - e/2)]') . (85) 

We note that for £ 7^ 0, these functions differ from the ex- 
pressions appearing on the right-hand side of Eq. ([84| . 
Howe ver, numerical simu lations and experimental data 
show ( Dietz et aLl . 120093 ) that for all values of F/d the 
last curly bracket in Eq. (|84)) (denoted by C^^^(e)) is well 
approximated by cil\o) F^l\e) / F^^\o) , and that simi- 
larly we have afiKSi^"" ''^*{E+e/2)\si'^'' ^)(£;-e/2)|2)] « 
(3) 

Faa (s). The last relation holds with good accuracy ex- 
cept for the regime F « d of weakly overlapping reso- 
nances. With these results, Eq. (|84| takes the form 



p(3)/ 



CifiO) ,(4) 



(86) 



All terms in Eq. (|86| are known analytically. Expres- 
sions useful fo r a numerical comp utation are given in the 
Appendix of (jPietz et all l2009cl) . Thus, from a practi- 
cal point of view theoretical expressions for cross-section 
fluctuations are available for all values of T/d except for 
the elastic case where the relation (|86p docs not hold for 
F«d. 



F. Poles of the S-Matrix 

Resonances correspond to poles of the 5'-matrix. The 
distribution of the poles of the stochastic S-matrix de- 
fined in Eq. (|4ip has, therefore, attracted theoretical at- 
tention from the beginning. Obvious questions are: How 
is the correlation width (|lll) related to the distance of the 



poles from the real axis? Is it possible to verify quan- 
titatively the picture drawn of the pole distribution in 
Section IEb]/ 



It was mentioned in S ectionHlFCl that lMoldaue^ (Il96ll 
I1963L 11964 I1969L I1975L Il976l Il980l) based his approach 



to CN scattering on the pole expansion of the S'-matrix. 
He seems to have been the first aut hor to determin e 
the distribution of poles numerically ( Moldaueil Il964l) . 
He made a number of important discoveries that stim- 
ulated later research, (i) There exists a gap separating 
the poles from the real axis, (ii) For strong coupling to 
the channel s, some poles oc cur far away from the real 
energy axis (lMoldaueilll975D. (ii i ) A "sum rule" for res- 
onance reactions (see ( Moldaueil Il969( ) and further ref- 
erences therein) relates the transmission coefficients and 
the mean distance of the poles from the real axis. 

Later work by various authors has led to a deeper 
understanding of these results. We begin with the 
"Moldauer-Simonius sum rule fo r resonance reac tions" 
and follow the derivation due to ISimoniu i (|1974D . We 
assume that the unitary S'-matrix defined in Eq. (|4ip 
has N simple poles £^ only. (Coincidence of two poles is 
considered fortuitous and actually excluded by quadratic 
repulsion of poles, see below). Then 



detSCN=exp{2j0}[]^|-|^ 

H/ On 



(87) 



The denominator on the right-hand side represents the N 
poles. The form of the numerator follows from the unitar- 
ity of S™: det[(S"^^)*] is inverse to det[S'CN]_ xhe only 
energy dependence is due to the poles of S'^^ , and the 
phase 4> is, therefore, constant. With = — (1/2)F^ 
we have 



Indet S™-2i</) = ^ln 1 



E-£,. 



To average Eq. (|88p over energy, we use a Lorentzian 
averaging function with width /, see Section IIV.AI and 
expand the logarithm in powers of Fp/(£^ — £^). Contour 
integration shows that only the linear term gives a non- 
vanishing contribution. With / ^ F^ for all fi we find 



In 1 - i- 



(89) 



E-£f,J/ E-^{£f,) + iI 
This yields 

(IndetS™) + (ln(detS™)*) 

" ^^2^ (£;-g?(£:^))2 + /2 - d ' 

We use Eq. ((58]) (i.e., the equality of (In dot S"*^^) and 
of lndet(S'^'^)) and the definition (g]) and obt ain the 
Moldauer-Simonius sum rule 



(91) 
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Eq. (|9T|) implies (F^) > T where T is the correlation 
width given by the Weisskopf estimate (|lip . Moreover, 
(r^) diverges whenever a single (or several) transmis- 
sion coefficients approach unity. The divergence is caused 
by the fact that one (or several) pole(s) of S™ is (are) 
shifted far below the real i?-axis, see Section IIV.BI and 
the text below. In deriving Eq. (|9l]) we have assumed 
that the averaging interval / is large compared to 
for all /i. It is not clear from the derivation whether 
the Moldauer-Simonius sum rule applies as one of the 
Ts approach es unity. However, wo r k on the unitary 
case (GUE) (iFvodo rov and Somrneri. 1 19961) and on the 
single-channel case for the GOE (|Sommers et all Il999l ) 
has shown that the distribution of the poles acquires a 
tail that causes the divergence as Tc 1- 

The locations of the poles of S are given 

by the eigenvalues of the effecti ve Ha milto - 

. ISokolov and Zelevinskvl (|1988L Il989l I1992D 
in an effort to extend statistical spec- 
(as based on the properties of g^Q^) to 
resona nces. Numerical work of iKleinwachter and Rotted 
(I1985D liad amplified Moldauer's observation that for 
large coupling to the channels, one or several poles is 
(are) located far from the real energy axis. Sokolov and 
Zelevinsky gave a semiquantitative analytical explana- 
tion of that observation summarized in Section IIV.BI 
and showed that at the same time, the majority of poles 
moves back towards the real axis ( "trapped states" ) . For 
the single-channel case they derived the distribution of 
poles in the complex energy plane. They connected the 
existence of one or several poles far from the real energy 
axis with the phenomenon of supcrradiance ( Dickd . 
1 1954) in quantum optics. They showed that the non- 
Hermitean part of causes quadratic repulsion of the 
poles in the complex plane, and that finding the density 
of poles of S in the complex energy plane is equivalent 
to the reconstruction of the two-dimensional charge 
density from a given electrostatic field. 

To determine analytically the joint probability density 
of the poles in the complex energy plane for the ^-matrix 
defined in Eq. (j4T|) . two approaches have been taken. If 
the limit — > oo is taken with A fixed, terms of order 
m = A/N do not contribute to the saddle-point equa- 
tion, neither in the replica approach (see Section IV. B.2|) 
nor in the supersymmetry approach (see Section FV.Cp . A 
mo del different from but re l ated to Eqs. (1411) i ntroduced 
by ISokolov and Zelevinskvl (Il988l . Il989l Il992[) makes it 
possible to overcome that limitation, and to discuss the 
pole distribution in the framework of the saddle-point 
approximation. The number A of channels is assumed 
to be large and to go with N to infinity while the ratio 
m = A/N is held fixed. The parameters Wafi describing 
the coupling of level /U with channel a are taken to be 
Gaussian-distributed random variables with mean value 
zero and common second moment 



TT 



(92) 



semblc of S'-matrices is invariant with respect to orthog- 
onal transformations of the channel space. The dimen- 
sionless strength 7 of the coupling is the only parame- 
ter. In contrast to the use of ijGOE -j^ -^j^g 
invariance of the distribution of the Ws under orthogo- 
nal transformations of the channels cannot be deduced 
from quantum cha os and is, th erefore, somewhat arbi- 
trary. In (pLehmann et al\ . \\QQba^ \ it was shown, however, 
that the pole distribution obtained from the model ((92)) 
is quite similar to the one where the Ws are fixed. The 
limit A — ^ 00 corresponds to the Ericson regime. How- 
ever, by choosing 7 <C 1, one can approach the limit of 
weakly overlapping resonances. 

The model ([92]) was used for an extensive dis- 
cussi on of the distri bution of poles of the iS-matrix 
in (iHaake et al. ^ 19921) (where th e replica trick was used) 
and in (jLehmann et al\ . Il995a[ ) (where supersymmetry 
was applied). In both cases, the saddle-point equations 
differ from those obtained for fixed channel number by 
terms of order m = A/N. These equations are used to 
determine the average pole distribution in the complex 
energy plane with the help of an electrostatic analogy 
similar to the one mentioned above. Analytic expres- 
sions are obtained for the boundary curve separating the 
area of non-vanishing pole density from the empty one. 
Typical results are shown below. 

As mentioned above, Moldauer discovered a gap sep- 
aratin g the poles of S from the real energy axis. 
Later. iGaspard and Ricd ( 1989al |bl) deduced the existence 
of a gap in the framework of the semiclassical approx- 
imation for chaotic systems with few degrees of free- 
dom. That suggests that the gap is a universal feature 
of chaotic scattering. For the model of Eq. ([^ . the gap 
was shown to exist and the gap parameters were worked 
out an alytically in ( Haake et al\ . Il992l : iLehmann et al\ . 
Il995a[) . In the center of the GOE semicircle, the width 
Fgap/2 of the gap se parating the cloud of p oles and the 



real axis is given by ( Lehmann et a/.l . ri995at ) 



F, 



gap 



^±A- 
27r 1 



47 



(93) 



The A channels are all equivalent, and the resulting cn- 



Using Ea. (j59|) and Xa = 7 we can write the right- 
hand side as (d/(27r)) Tr. This is the Weisskopf es- 
timate. In ( Lehmann et al\ . 19953), the S-matrix corre- 
lation width for the model of Eq. ([^^ was also worked out 
and found to coincide with the one found in the frame- 
work of Eqs. (|4T|) (fixed number of channels), i.e., with 
the Weisskopf estimate. For the model of Eq. ([U^ . the 
equality of gap width and Weisskopf estimate holds un- 
less m approaches unity. Then, the correlation width be- 
comes comparable to and is modified by the range of the 
GOE spectrum. Put differently and positively, the result 
implies that S-matrix fluctuations are universal as long 
as there is a clear separation of the two energy scales, the 
gap width and the range of the GOE spectrum. 

The joint probability density of the poles of S in 
the complex plane can be obtained from the above- 
mentioned algebraic equations. With ^{E)/\ = 
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'^{E)/X = y, Fig.[5]shows the distribution below the real 
i5-axis. The single cloud seen for 7 = 0.2 and for 7 = 1 
splits into two as 7 is increased further. The separation 
begins to develop at 7 = 1. From the Moldauer-Simonius 
sum rule, we would expect that many poles move to —ioo 
as 7 approaches unity. That is not seen in the figure. We 
ascribe the discrepancy to the assumption (used in the 
derivation of the sum rule) that the entire energy depen- 
dence of S™ is due to the poles, see Eq. (gT]). That 
assumption fails when the widths of the resonances be- 
come comparable with the range of the GOE spectrum. 



the distribution of poles of the scattering matrix in 
the complex energy plane employed a microwave res- 
onator and the method of harmonic analysis, see also 
Sectio n IVLBl The res ults agree with theoretical predic- 
tions ( Sommers et al\ . [l99^ . 




FIG. 5 Distribution of the poles of the 5'-matrix in the com- 
plex energy plane. The distribution is shown for several val- 
ues of the strength 7 of the average coupling to the chan- 
nels. The abscissa spans th e entire range of the spectrum. 
From iLehmann et al\ (|l995al ) . 



The second approach (iFvodorov and Koruzhenkol 
119991 : iFvodorov and Sommersl . l2003l ) takes the limit N 
00 for fixed channel number A and arbitrary values of the 
transmission coefficients Tc ("almost Hermitian matri- 
ces"). For the unitary case (GUE Hamiltonian) it yields 
the complete joint probability density of the eigenvalues 
of the effective Hamiltonian in Eq. Corresponding 
results for the GOE are not known except for the single- 
channel case A = 1 (ISokolov and Zelevinskvl . [l989l: lUllahl . 
I1969D . 

The first experimental study ( Kuhl et all . l2008l) of 



In summary: The distribution of the poles of S in 
the complex energy plane gives valuable insight into the 
scattering mechanism, even though that information is 
not sufficient to construct the scattering amplitude(s). 
The Moldauer-Simonius sum rule shows that in general 
the average distance of the poles from the real axis is 
bigger than would be concluded from the Weisskopf esti- 
mate pT|) . With increasing coupling to the channels, up 
to A poles are moved ever further away from the real axis, 
a fact related to superradiance in quantum optics. An ex- 
ample for the actual distribution of poles in the regime of 
many strongly coupled channels {N — 00, A — > 00 and 
m = A/TV fixed) is shown in Fig[S] A gap separates the 
cloud of poles from the real i5-axis. The width of the 
gap is given by the Weisskopf estimate. The regime of 
almost Hermitean matrices has also been worked out. 



G. Correlations of ^-Matrix Elements carrying different 
Quantum Numbers 

We return to the question raised in the last paragraph 
of Section fill. Bl Is the assumption ([5]) justified that S- 
matrix elements carrying different quantum numbers are 
uncorrelated? 

While RMT per se is obviously not in a position to 
give an answer to that question, a realistic large-scale 
shell-model calculation would. We have in mind a cal- 
culation using the two-body random ensemble (TBRE) 
of the nuclear shell model. The TBRE was introduced in 
Section I.V.B where details and references to the original 
papers may be found. In the TBRE, several neutrons and 
protons occupy the single-particle states of a major shell 
of the shell-model and interact via a two-body interac- 
tion. The matrix elements of that interaction are taken 
to be Gaussian-distributed random variables. Spin is a 
good quantum number, and the need to couple nucleon 
angular momenta and spins to good total spin creates 
considerable complexity. Matrix elements of S carrying 
different spin quantum numbers would have to be calcu- 
lated numerically for different realizations of the TBRE, 
and their correlations worked out. For reasons given be- 
low, the dimensions of the underlying shell-model spaces 
would have to be very large. We are not aware of any 
such calculation. 

H owever, in the conclusions 

of ( Papenbrock and Weidenmiillerl . |2007[ ) the au- 
thors argue that for sufficiently large shell-model 
spaces, the correlations between S'-matrix elements 
carrying different quantum numbers might be very 
weak. The arguments are based on t he st udy of a 
model ( Papenbrock and Weidenmiillerl . |2007| ) that is 
conceptually close to but technically simpler than the 
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TBRE and, therefore, analytically accessible. In the 
model parity is the only quantum number, and m spin- 
less fcrmions occupy £i {£2) degenerate single-particle 
states with positive (negative) parity, respectively. The 
fermions interact via a parity-conserving two-body 
interaction with random Gaussian-distributed uncor- 
related two-body matrix elements. The nth moments 
A^„(±) of the Hamiltonian H of the model (defined as 
normalized traces of i/" with positive integer n) can 
be worked out analytically for the many-body states of 
both positive (-I-) and negative (— ) parity. The case of 
large matrix dimension is attained in the "dilute limit" 
defined by £i,£2,"i — > 00, m/li 0, — >■ 0. It is 
shown that Mni+) — Mn{^) for all n up to a maximum 
value that is bounded from above by m but tends to 00 
in the dilute limit, and that the two moments differ ever 
more strongly when n grows beyond that bound. That 
result shows that the spectra of states with positive and 
negative parity are strongly correlated. In particular, the 
shapes of the two average spectra are extremely similar. 
At the same time, the result suggests that the local 
spectral fiuctuations of the two spectra are uncorrelated. 
Indeed, in their work on the use of moments for nuclear 
spectroscopy, French and collaborators concluded that 
such fluctuations are determ ined by the very highest 
moments of the Hamiltonian ( Brodv et a/.l . [ToSlI ) . 

The results just stated apply in the dilute 
limit only. They do not contra di ct ea rlier find- 
ings ( Papenbrock and Weidenmiilleil . [2003) for the 
TBRE on correlations between spectra carrying different 
quantum numbers. These calculations involved Hilbert 
spaces of small dimension only. 

Assuming that a result similar to the one just stated 
holds in the limit of large matrix dimension for the TBRE 
and observing that the fluctuation properties of the S*- 
matrix are caused by the local spectral fluctuation prop- 
erties of the underlying Hamiltonian (see Eqs. (|28|29p ). 
we conclude that within the framework of the nuclear 
shell model, S'-matrix elements carrying different quan- 
tum numbers are likely to be uncorrelated, in agreement 
with Eq. ©. 



VI. TESTS AND APPLICATIONS OF THE STATISTICAL 
THEORY 

The statistical theory reviewed in Sections IIVI and |V] 
has been tested thoroughly. Moreover, it has found nu- 
merous applications, both within the realm of nuclear 
physics and beyond. In this Section we review some re- 
cent such tests and applications. 



A. Isolated and Weakly Overlapping CN resonances 

In the regime of isolated resonances, thorough tests 
of the statistical theory of nuclear reactions (F <C d) 
were undertaken already many years ago. Especially 



for ne utron reson ances there exists a comprehensive re- 
view ( Lvnnl . [l968l ). In the regime of weakly overlapping 
resonances, tests have so far not been performed in nu- 
clei. Here we do not summarize the early works but 
rather focus attention on recent data and tests of the 
theory. These have bec ome possible in microwave bil- 
liards ( Dietz et al\ . \206§t] . Such devices simulate the CN 
and its resonances or, for that matter, any other chaotic 
quantum scattering system. Indeed, in sufficiently flat 
microwave resonators and for sufficiently low values of 
the radio frequency (rf ) — so-called quantum billiards — 
only one vertical mode of the electric field is excited, and 
the Helmholtz equation is mathematically equivalent to 
the Schr odinger equation for a two-dimensional quantum 
billiard (iGraf et all Il99ll: iRichteil . '199* 'So e t aLl. 11995 ; 
Sridhad . [T99ir IStockmann and Stein . .1990; .StockmannI 



2OOOI ). If the classical dynamics (free motion within the 



microwave resonator and elastic scattering by its bound- 
ary) is chaotic, the statistical properties of the eigen- 
values and eigenfunctions of the closed r esonator in the 
quan tum case follow RMT predictions ( Bohigas et all 
1 19841) . and the scattering of rf amplitudes by the res- 
onator corresponds to quantum chaotic scattering. 

rf power rf power / 
in out / / 




FIG. 6 Flat microwave resonator (left-hand side) as a model 
for the compound nucleus (right-hand side). The height 
d = 0.84 cm of the flat resonator makes it a chaotic quan- 
tum billiard up to a frequency of 18.75 GHz. 

The left-hand side of Fig. [5] shows a typical quantum 
billiard realized in the form of a flat microwave resonator. 
For the measurement of the spectrum, rf power is cou- 
pled via an antenna labeled 1 into the resonator, thereby 
exciting an electric field mode within the resonator, and 
the reflected output signal at the same antenna (or the 
transmitted one at the antenna labeled 2) is determined 
in magnitude and phase in relation to the input signal. 
Hence, the resonator is an open scattering system where 
the antennas act as single scattering channels. The scat- 
tering process is analogous to that of a CN reaction as 
indicated schematically on the right-hand side of Fig. [51 
The incident channel A + a consists of a target nucleus 
A bombarded by a projectile a leading to a compound 
nucleus which eventually decays after some time into the 
channel with the residual nucleus B and the outgoing 
particle b. (We disregard angular momentum and spin). 
Attaching more antennas to the resonator (or dissipat- 
ing microwave power in its walls) corresponds to more 
open channels C + c, D + d, ... of the compound nu- 
cleus. The resonator in Fig.lHlhas the shape of a so-called 
Bunimovich stadium billiard which is known to be fully 
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chaotic in the classical limit ( Bunimovicbl . Il985l ). 
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FIG. 7 Part of the transmission spectrum of a superconduct- 
ing microwave billiard. We note the extremely high resolution 
of the resonances. Two singlets (a) and (c) and a doublet (b) 
of resonances are magnified in the upper part. For a ll three 
cases 7i-matrix reso nance formulas (iBeck et all |2003| ) based 
on expressions f rom iLane and Thomas (|l958l ) were fitted to 
the data. From (|Dembowski et a/.ll2005l ). 



Here and are the real and imaginary parts of 

the eigenvalues of an effective Hamiltonian ifoff = H ~ 
i TT WW^ for the microwave billiard in vifhich the Ws de- 
note the coupling of the resonator states to the antenna 
states and to the walls of the billiard, see Section llII.A.2l 
With Wa°J — >■ Tq^^, Eq. ([M]) is completely equivalent to 
Eq. (|3ip which describes the coupling of isolated nuclear 
quasi-bound states to the channels. 

As discussed below Eq. (|5T|) the partial width am- 

1/2 

plitude TfJa is the probability amplitude for decay of 
resonance fi into channel a, and = F^a is the 
total width. In an experiment with the superconduct- 
ing quantum billiard (shown on the left-hand side of 
Fig. ini) coupled to three antennas c = 1,2,3 (corre- 
sponding to a CN reaction with three open channels), 
the complete set of resonance parameters (resonance en- 
ergies, partial widths, to t al wid ths) for 950 resonances 
was measured ( Alt et all . Il995f ). The total widths F^ 
and the partial widths F^2 are found to fluctuate ran- 
domly about a slow secular variation with fi, i.e., with 
frequency. The GOE predicts a Gaussian distribution 
for the decay amplitudes vjl^ or, equivalently, a dis- 
tribution wit h one degree of freedom (th e Porter-Thomas 
distribution ( Porter and ThomasllT956l ). see Eq. (F18)). 



The distribution of F^2 in Fig. [5] exhibits this behavior 
impressively. 



The resonator in Fig.lHlwas made of niobium and oper- 
ated in a superconducting mode. This strongly increases 
the quality factor of the resonator and yields very high 
resolution. The resonator was used for both, the study of 
spectral prop erties, i.e., the sta tistics of the resonances 
in the cavity ( Graf "oil. 1991), an d a measurement of 



their decay widths ( Alt et al . 19951) . The latter are pro- 



portional to the square of the billiard eigenfunctions at 
the locations of the antennas. Figure [7] shows a transmis- 
sion spectrum. Spectra of that quality can be typically 
obtained in superconducting billiards in the regime of 
isolated and weakly overlapping resonances. The mea- 
sured ratio of Pout,b, the rf power signal transmitted into 
antenna b, and of Pin, a, the incoming rf power signal at 
antenna a, is shown on a semilogarithmic plot. The ra- 
tio is equal to \Sab{f)\'^ where Sab{f) with a,b — 1,2 are 
the elements of the complex-valued frequency-dependent 
2x2 scattering matrix S. More generally, measurements 
of the modulus and phase of the outgoing and the in- 
coming signals performed with a network vector analyzer 
determine magnitude and phase of all elements Sab{f) of 
S. Such detailed information is not usually available for 
other chaotic scattering systems where in general one can 
only measure intensities. For a number of isolated res- 
onances labeled /x with fi = I, . . . , N without any back- 
ground scattering the S'-matrix is a sum of Breit-Wigner 
terms. 



Sab — Sab — i 
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(94) 



FIG. 8 The experimental distribution of the the partial 
widths r,i2 in units of their mean value (r^2). The 
solid hue correspond s to a Porter-Thomas distribution. 
From (|Alt et alivm^ . 



As a further test of the statistical theory, the autocor- 
relation function of the S'-matrix of the cavity, defined 
as Ce(e) = (Se,(/) S:^{e + /)) - \{S,,{f))\\ was deter- 
mined for channels c = 1,2, and 3. The average is taken 
over frequency. For c = 2 the result is plotted as cir- 
cles in the upper part of Fig. HI The shaded band is a 
measure of the experimental uncertainty of C2(e). The 
shape of the shaded band differs markedly but not unex- 
pectedly ( Lewenkopf and Weidenmiilleil . [l99ll ) from that 
of a Lorentzian with width (F^) shown as a solid line in 
the upper part of Fig. [51 The contribution of each of 
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the over 900 individual resonances in Fig.|9]is, of course, 
Lorentzian in shape with width F^. However, different 
resonances have different widths, and the average over all 
resonances is not a Lorentzian. This happens only when 
the total widths fluctuate strongly with Thus, 
both the number of open channels and the absorption 
in the walls of the microwave billiard (also contained 
in F^) must be small. Both conditions can be satis- 
fied in experiments with superconducting microwave res- 
onators, but generally not in experiments at room tem- 
peratu re. (In the microwave experiment of lDoron et al\ 
( 1990l ). for instance, absorption was strong and conse- 
quently the_data^_did_jiot display a non-Lorentzian line 
shape (iLewenkopf et al. I, I1991D .) Clearly, the observed 
non-Lorentzian shape is a quantum phenomenon: In the 
semiclassical approximation, i.e., for many open chan- 
nels, we have purely exponential decay. The result dis- 
played in the upper part of Fig.lHlis in quantitative agree- 
ment with the statistical theory, i.e., with Eqs. ([75]) to 

dZT]). 




t {1/<T>) 

FIG. 9 Upper part: The experimental autocorrelation func- 
tion |C2(e)P (circles within the shaded band of errors), 
the prediction of the statistical theory (dashed line), and a 
Lorentzian (solid line). Lower part: Fourier coefficients of the 
autocorrelation function (dots) with errors indicated by the 
shaded band together with the prediction of the statistical 
theory (dashed line) and an exponential (Fourier transform 
of a Lorentzian). The non-exponential decay in time of the 
Fourier transform of th e ^'-m atrix autocorrelation functions 
is clearly visible. From (|Alt et a/.. . ,19951 ). 



The Fourier transform (FT) of the autocorrelation 
function C2 (e) shown in the lower part of Fig. |9] decays 
non-exponentially as a function of time. This feature, 
too, reflects the non-Lorentzian shape of the autocorre- 
lation function shown in the upper part of Fig. [9] and 
is in agreement with the statistical model. We refer to 
the discussion in the paragraph following Eq. ([78| . One 
may say that in the experiment non-exponential decay 
of a quantum system with chaotic dynamics has been 
"observed" for the first time. 

A thorough experimental investigation of chaotic scat- 
tering in microwave billiards has recently also been per- 
formed in the regime of weakly overlapping resonances 
(F < d), and the results have been analyzed with the 
help of the formulas of Section IV. CI A normally con- 
ducting fiat microwave res onator mad e of c opper was 
used as a quantum billiard ( Dietz et~al\ . l2008l ) . The res- 



onator shown i n the inset of Fig. fTOlhas the shape of a 
tilted stadium ( Primack and Smilanskvlll994[ ). The sta- 
dium has fully chaotic classical dynamics. The resonator 
carried two antennas, and the complex elements of the 
symmetric scattering matrix Satif) were measured ver- 
sus frequency. Figure [TU] shows examples of the measured 
transmission (|S'i2p) and reflection (IS'np) intensities as 
functions of resonance frequency. We note the strong 
fluctuations of both quantities with frequency. 




FIG. 10 Transmitted (upper part) and reflected (lower part) 
intensity versus frequency between 9.0 and 9.5 GHz. The reso- 
nances overlap and create a fluctuation pattern. The shape of 
the two-dimensional quantum billiard used in the experiment 
is shown in the inset. Points 1 and 2 in dicate the positions of 
the antennas. From (|Dietz et all\206^ . 

As in the case of isolated resonances the S'-matrix au- 
tocorrelation function Cab = (Sabif) S'*^(/+£)) - 1 (S'ab) P 
for a, & = 1 or 2 was computed from the data. The upper 
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part of Fig. [TT] shows that the values of the scattering 
matrix Sabif) are correlated with a correlation width F 
of order several MHz. The values of Sab{f) measured 
at M equidistant frequencies with step width A have 
been Fourier-transformed. The complex Fourier coefS- 
cients are denoted by Sabit) with t > 0. We use the 
discrete time interval t = d/AIA elapsed after excita- 
tion of the billiard resonator instead of the Fourier index 
k. The Fourier transform Cab{t) of Ca6(e) has Fourier 
coefficients \Sab{t)\'^ and any two coe fficients are uncor- 
rclated random variables ( Ericsonlll96d) . The lower part 
of Fig. [TT] shows data for log;^o (^ab{t)- We note that the 
Fourier coefficients scatter over more than five orders of 
magnitude. The cutoff a.t t = 800 ns is due to noise. The 
decay in time is strikingly nonexponential, i.e., powerlike, 
as for isolated resonances (Fig. [5]). 




4 6 
e(MHz) 




FIG. 11 Upper part: The autocorrelation function \Cab{£)\^ 
for values of a and b as indicated. The values of |Cab(e)P were 
calculated from the measured S-matrix elements (points) and 
from the fit of Eq. (|75|) of the statistical theory to the data 
(solid line). Lower part: Fourier coefficients Cab{t) of the 
autocorrelation function (points) and the Fourier transform of 
the fit of Cabie) (Eq. ((75])) to the data (solid line). The non- 
exponential decay in time in the region of weakly over l appin g 
resonances (F/d ~ 0.2) is striking. From (|Dietz et a/.ll200"^ . 



The solid lines in Fig. [TT] result from fitting Eqs. ([75)1 



to ([77)) of the statistical theory to the data points. The 
parameters in Eqs. ([75]) to ([77)1 for the S'-matrix autocor- 
relation function Cab(e) are the transmission coefficients 
Ta with a = 1, 2 for the open antenna channels, Tc with 
c = 3,4,... for ad ditional fictiti o us ch annels modeling 
Ohmic absorption ( Schafer et al\ . l2003l ) in the walls of 
the normally conducting resonator, and the average level 
spacing d. For a = 1 , 2 the transmission coefficients were 
calculated from Ta — 1 — \{Saa)\'^ with data on Saa as 
input. The averag e level spacing d was c alculated from 
the Weyl formula (|Baltes and Hili [1976[) . The product 
in Eq. (|75|) over the ficticious channels was replaced by 
an exponential function of the sum Tabs of the transmis- 
sion coefficients of these channels. This is a good ap- 
proximation when all Tc <C 1 and left Tabs as the only 
free parameter in the fit. In order to allow for secular 
variations of Tabs the experimental data were analyzed in 
1 GHz intervals. It was found that the sum Ti -I- T2 -I- Tabs 
increases from the value 0.11 in the interval 3-4 GHz to 
the value 1.15 in the interval 9-10 GHz. (The resulting 
increase of Tabs is consistent with conductance properties 
of copper.) Using these numbers and the Weisskopf es- 
timate F « [d/27r] X^c^c d) correlation 
width F one finds that T/d increases from 0.02 to 0.2 over 
the same frequency range. Thus, the statistical theory 
of chaotic scattering was indeed experimentally tested 
in the regime of weakly overlapping resonances. The 
Fourier coefficients turn out to be uncorrelated, Gaussian 
distributed random variables. That fact and the large 
number of such coefficients (2400 per frequency interval) 
made it possible to assess the quality of the agreement 
between data and the fits in terms of a Goodness-of- Fit 
(GOF) test, with excellent results (lDietzIiraI|,[200i). 



The experiment of lDietz et al. I (I2OO8D has produced in- 
teresting results also on the distribution of moduli and 
phases of S'-matrix elements (see Section IV.D|1 , and on 
the elastic enhancement factor, see Sections IV. Al and 
IV.B.ll Concerning the first point, we expect theoretically 
that for F ^ d the distribution of S'-matrix elements is 
non-Gaussian. This is because unitarity constrains the 
distribution of S'-matrix elements. The constraints are 
strongest for F <C d, see Section [V. El Fig. [12] shows that 
the distribution of the real part of Sn is strongly peaked 
near unity, especially for the lower frequency interval 
from 5 to 6 GHz in which mostly isolated resonances are 
found. But even in the regime of weakly overlapping reso- 
nances, i.e., from 9 to 10 GHz the distributions of the real 
and imaginary parts of Sn deviate from Gaussians (solid 
lines). The distributions of the phases (rightmost pan- 
els in Fig. [T2)) are peaked. However, the valley between 
the two peaks fills up as T/d increases. For T/d ^ 1, 
i.e., in the Ericson regime, the phases are expected to be 
uniformly distributed and the S-matrix elements to be 
Gaussian distributed. The elastic enhancement factor, 

defined as I^ = {{\Sf,\'') {\S^^\'')Y'^ / {\Sf^\') , is deter- 
mined from the data as a function of / in two ways. One 
may either use the autocorrelation functions (Fig. Ill|) or 
the widths of the distributions of the imaginary parts 
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of the scattering matrix. Both results agree and yield a 
smooth decrease of W with / from W « 3.5 ± 0.7 for 
4 < / < 5 GHz toW^ 2.0 ± 0.7 for 9 < / < 10 GHz, 
the errors being finite-range-of-data errors. The calcu- 
lation of the enhancement factor using Eqs. ((75)) to ([77)) 
gives the values W = 2.8 and 2.2 for the respective fre- 
quency intervals. Wc return to the determination of W 
once more in Section IVlI.CI below. 







FIG. 12 From left to right: Histograms for the scaled dis- 
tributions of the real and imaginary parts of the reflection 
amplitude Sn and the real part and the phase of the trans- 
mission amplitude 5*12, respectively. The data were taken in 
the two frequency intervals 5-6 GHz (upper panels) and 9- 
10 GHz (lower panels). The scaling factors are given in each 
panel. The solid lines are best fits to Gaussian distributions. 
From ()Dietz et aU 120081 '). 

In summary, we have reviewed in this Section how a 
wealth of very precise experimental data on chaotic quan- 
tum billiards that mimic the CN and its resonances, can 
be used for a stringent test of the theory of quantum 
chaotic scattering. Measurements of the phases of the 
scattering matrix elements provide valuable additional 
information that is usu ally not accessib le in CN reac- 
tions. In the work of lAh et al\ (|1995D on the statis- 
tics of partial widths of isolated resonances (F ^ d), 
all tests applied — Porter-Thomas distribution, lack of 
correlations between partial widths in different scatter- 
ing channels, lack of correlations between partial widths 
and resonance frequencies, non-Lorentzian decay of the 
S-matrix autocorrelation function, and the related non- 
exponential time decay of the chaotic quantum system — 
are in perfect agreement with GOE predictions for the 
statistics of eigenfunc tions and eigenvalu es. In a contin- 
uation of this work by jPietz et aLl (l2008l ) into the regime 
of weakly overlapping resonances (F « d), the distribu- 
tions of S'-matrix elements are found to be non-Gaussian 
while the Fourier coefficients of these S'-matrix elements 
do have an approximately Gaussian distribution. These 
data were used for a highly sensitive test of the statis- 
tical theory reviewed in Sections IIVI and |V| In particu- 
lar, the predicted non-exponential decay in time of iso- 
lated and weakly overlapping resonances and the values 
of the elastic enhancement factors are confirmed. The 
evidence for non-exponential decay in time, obtained by 
Fourier-transforming measured S'-matrix elements into 
time space, is still indirect. A direct measurement of the 
decay time of an excited nucleus might become possible 
at high-power laser facilities such as the National Igni- 
tion Facility (NIF) where all the nuclear resonances (or 



subsets of them) might be excited simultaneously by a 
short laser pulse (jMoses et all |2009() . 



B. Strongly Overlapping CN resonances: Ericson 
Fluctuations 

In the 1960s and 1970s, the newly discovered phe- 
nomenon of Ericson fluctuations formed a central part 
of research in nuclear reactions, see Section III.Dl Pro- 
tons, deuterons and light ions up to oxygen or so, 
but also fast neutrons have been used as projectiles 
in various C N reactions. The field h a s bee n re- 
viewed earl v (lEricson and Maver-Kucku5 Il966[ ) and 



again later ( Richtei . 1974 ). Both articles show that 



all of Ericson's predictions were confirmed experimen- 
tally. We do not reiterate here what has been known 
for many years about cross-section autocorrelation func- 
tions, their Lorentzian shape, the mean coherence width 
F, the vanishing of cross-correlation functions between 
cross sections in different reaction channels in the ab- 
sence of direct reaction contributions, Hanbury Brown- 
Twiss behavior of cross-correlation functions of cross sec- 
tions measured at different scattering angles, probability 
distributions of randomly fluctuating cross sections, and 
about the interplay between direct and CN reaction pro- 
cesses in general. From the few experiments performed 
lately in nuclear physics, two are chosen to exemplify over 
and above what has been stated in Section II. D why Eric- 
son fluctuations are now commonly viewed as a paradigm 
for chaotic behavior of a quantum system. 

Recently, a number of mainly neutron-induced CN 
reactions on medium-heavy nuclei has been studied in 
the Ericson regime, primarily in o rder to deduce nu- 
clear level densities from the data ( Grimed ( 2002) ) and 
Refs. therein). A striking example, the excitation 
function for the reaction^^Si( n, iPo4-i)^^Al, is shown in 
Fig. [13] ()Bateman et all )1997)) . Decay of the CN ^^Si 
populates the ground state and the first excited state of 
the final odd-odd nucleus ^^Al. The two proton channels 
Po+i are not resolved. The high-resolution measurement 
of the CN cross section reveals significant fluctuations 
with energy. The peaks and minima of the excitation 
function are not caused by individual, more or less iso- 
lated resonances, but instead result from the construc- 
tive (or destructive) superposition of many overlapping 
CN resonances. The amplitudes of the resonances are 
random variables. Therefore, the curve connecting the 
measured points in Fig. [13] has the curious feature of be- 
ing both reproducible and random ([Wci d cnmiillcr, . ,1990.) . 
A measurement of the same reaction in the same energy 
interval with the same energy resolution will reproduce 
Fig. [13] Nonetheless, the energy dependence of the cross 
section in the figure displays the features of a random 
process. 

The second example for the role of Ericson fluc- 
tuations in nuclei is from the field of giant reso- 
nances where the question of direct versus statisti- 
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FIG. 13 Fluctuating excitation function of the reaction 
^*Si(7i,po+i)^**Al in the regime of strongly overlapping res- 
onances (r » d) for the CN ^^Si. The cross sections of 
proton exit channels po and pi leading to the ground and 
first e xcited states in '^^Al, respectively, were not resolved. 
From (|Bateman et aZ.I . [l997l ) . 



cal decay plays a central ro l e (iBortignon et all Il998l : 
iHarakeh and van der Woudd . l200lD . The giant reso- 
nance is a "distinct" and "simple" mode of excitation 
of the nuclear ground state whose amplitude is usually 
spread over many "complicated" states for which the dis- 
tinct mode acts as a "doorway". The strength function 
is typ ically of Breit-Wigner shape ( Bohr and Mottelsonl . 
Il969l ). Doorway states were reviewed in Section I.II.G. 
To describe the approach, we have to distinguish several 
contributions to the total wid th Fp of the gian t reso - 



nancc. In good approximation ( Goeke and Spethl . Il982l ) 
To can be written as Fq = AF + Fl' + F-^. Here AF 
stands for the Landau damping of the giant resonance. 
(Electric dipole excitation of the ground state in the first 
step produces a coherent superposition of one particle- 
one hole (Ip-lh) states that have different single-particle 
energies.) The Ip-lh states in turn couple to the contin- 
uum and acquire an escape width F^ which gives rise to 
a direct decay contribution into dominant hole states of 
the daughter nucleus. Finally, F^ describes the spreading 
width resulting from mixing the Ip-lh states with more 
complex 2p-2h and further np-nh configurations until an 
equilibrated compound nucleus is reached. A primary 
goal of all giant-resonance high-resolution decay exper- 
iments is, thus, to determine the relative contributions 
of the widths F^ and F-'' to the total width Fq. The ex- 
perimental signature of direct nuclear decay of a giant 



resonance in a nucleus with mass number A is the en- 
hanced population of hole states in the daughter nucleus 
A — 1. Statistical decay can be identified either by a 
comparison of the measured decay spectrum with pre- 
dictions of the Hauser-Feshbach formula (fTO)) or by an 
Ericson-fluctuation analysis of the fine structure in the 
decay spectrum measured with high resolution. Here we 
address the second possibility. 

Giant-resonance spectroscopy of ^°Ca has recently 
been performed through exclusive electroexcitation ex- 
periments of the type '*°Ca(e, e'x) where x stands for 
either protons or alpha particles detected in c oinci- 
dence with th e scattered electron ( Carter et al\ . l200ll : 
iDiesener et a/.l . [2001a .b). The description of the {e,e'x) 
reaction is based on the one-photon exchange mecha- 
nism. The incident electron is inelastically scattered on 
'^'^Ca, and a virtual photon with energy and momen- 
tum q is transferred to the ^''Ca nucleus exciting it into 
the giant-resonance region. The excited nucleus propa- 
gates in time and finally decays into a residual nucleus 
(■^^K, ^^Ar) by emitting particle x. It is assumed that 
excitation and decay can be treated independently. The 
upper part of Fig. [14] shows the double differential cross 
section for electrons that are scattered inelastically on 
■^'^Ca and measured in coincidence with protons that leave 
the residual nucleus '^^K in its ground state. The elec- 
trons were detected at an angle = 22°. The protons 
were measured under various angles. The total yield (in- 
tegrated over all angles) was determined. The range of 
excitation energies in ^°Ca is E^ « 10-27 MeV. For 
Ex between 10 and 15 MeV, a number of isolated states 
is observed but the most prominent excitation is clearly 
the MeV-wide peak at E^ « 19 MeV. It is the electric 
giant dipole resonance in ^"^Ca. Superposed upon the 
peak is considerable fine structure due to Ericson fiuc- 
tuations of the underlying overlapping CN resonances. 
Such fluctuations have been seen already some time 
ago ( Diener et aLl . 119731 ) in the reaction '^^K(p, 7o)^°Ca. 
The original spectrum has been smoothed with a Gaus- 
sian of FWHM = 800 keV (continuous solid line). The 
middle panel shows the ratio of the actual cross section 
and the smoothed one. The data fluctuate around the 
value unity. The autocorrelation function shown in the 
lower part was computed from the fluctuating cross sec- 
tion within 16 MeV < Er, < 23 McV (dashed lines in the 
middle part of the figure). The Lorentzian predicted by 
Ericson, C(£) = C(0) {T^y /{{T^y+e^), was fitted to the 
experimental points (open circles). This determines the 
spreading width with a value between about 15 to 30 keV , 
depending on the meth od of averaging ( Carter et all 



I2OOII : iDiener et aLl . [19731 ). The scatter of the points re- 
sults from the finite range of the data ( Richterj Il974l ). 



The value of C(e) at e = 0, i.e. the normalized vari- 
ance C(0) = {a'^) / {(t)"^ — 1 is related to the direct part 

= + for the reaction feeding the pro- 

ton decay channel (with S = {S)+S^) through the equa- 
tion C(0) = {1/nN) (1 — j/^). Here, N corresponds to the 
effective number of spin channels contributing to the re- 
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FIG. 14 Upper part: Double differential cross section of the 
■"'Ca(e, e'po) reaction at an electron energy Eq = 183.5 MeV. 
Middle part: Ratio of the measured cross section and the 
smoothed c ross section. Lowe r part: Autocorrelation func- 
tion. From ([Carter et aZ.L 120011 ). 



action, and n describes a damping factor due to the finite 
experimental energy resolution. The detailed analysis of 
the fluctuating cross sections in the ^°Ca{e, e'p)^^K re- 
action yields as fraction ya of the direct cross section a 
value between about 85-95 % for the po and pi decay into 
the ground and first excited states in ^^K, respectively. 
From the point of view of the shell model, these states are 



rations, respectively. Furthermore, the ratio of the cross 
sections for these two decay channels is close to the ratio 
of the single-nucleon transfer spectroscop i c fact ors from 



19761). Thus, 



the 40Ca(d,3 He)39K reaction (iDoll et all 

the analysis of Ericson fluctuations ( Carter et al\ . l200l[ ) 
has shown that in ^°Ca, the escape width F^ for direct 
proton emission of the electric giant dipole resonance 
feeding low-lying states of '^^K, is considerably larger 
than the spreading width F'''. We note that for the nu- 
clear giant dipole resonance, the ratio T^/T'^ strongly in- 
creases with mass number so that F^ dominates in heavy 
nuclei. 

We return once more to the analogy between a flat 
chaotic microwave resonator — a quantum billiard — and 
a CN. The data on billiards described so far relate to the 
cases of isolated or weakly overlapping resonances. An 
extension of the measurements with the chaotic tilted sta- 
dium billiard (inset of Fig. [TOl) into the Eri cson regime 
(F > d) was reported bv lDietz oil (|2009bf ). The abso- 
lute squares of the strongly fluctuating matrix elements 
S12 and 5*11, respectively, taken at high excitation fre- 
quencies, are plotted versus frequency in Fig. 1151 The 
autocorrelation functions and their Fourier transforms 
are shown in Fig. [1^1 The data show that already for a 
value of r/D sa 1.06 the system decays exponentially in 
time. Again, the decay pattern and the autocorrelation 
functions of the S'-matrix elements in Fig. [16] are very 
well described by the statistical theory. Furthermore, as 
shown in Fig. [17] for S12, the S'-matrix elements have 
a Gaussian distribution, and the phases are uniforml y 
distributed in the interval {0, 27r} (jPietz et all l2009bh . 
We note that the doorway-state phenomen on can also be 
mode lled in terms of a microwave billiard ( Aberg et all . 
l2008i). 

iKuhl et al. I (I2OO8D have taken a different approach to 
chaotic scattering in microwave cavities in the Ericson 
regime. They used the method of harmonic analysis to 
determine the locations of the poles of S in the complex 
energy plane. The resulting widt h distribution was com - 
pared with theoretical results bv ISommers et al\ (Il999[ ). 

These examples from recent experiments together 
with t he many others summarized in the earfier re - 
vicws ([Ericson and Maver-KuckukL Il966l : iRichted . Il974[ ) 
should suffice to demonstrate the importance of Ericson 
fluctuations in nuclei. The pheno menon occurs, however, 
also in other quantum systems. [Bliimel and Smilanskvl 
( 1988f ) analyzed numerically the effect of irregular classi- 
cal scattering on the corresponding quantum-mechanical 
scattering matrix. Using semiclassical arguments they 
showed that the fluctuations of the S'-matrix and 
the cross section arc consistent with Ericson fl uctua- 
tions. Wcidcnmiillcr (1990) and [Sorathia et'oH ([2009f ) 
compared universal conductance fluctuations of meso- 
scopic systems in the metallic regime with Ericson fluctu- 
ations of CN cross sections and pointed out the common 
stochastic features of the resonances in both phenomena 
but also the differences. These arise because the con- 
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FIG. 15 Same as in Fig. [10] but in the Er icson regime of 
overlapping resonances. From (|SchafeJ , [20091 ). 



ductance is a sum over many channels, and because the 
length of the mesoscopic system can be varied contin- 
uously. [Main and Wunneil ( 1992[ ) reported exact quan- 
tum calculations for the photoionization cross sections 
of the hydrogen atom in crossed magnetic and electric 
fields and found strong Ericson fluctuations as a char- 
acteristic feature of cha otic scattering. A featur e article 
in physical chemistry bv lReid and Reisleil ( 1996[ ) empha- 
sized the manifestation of interfering overlapping reso- 
nances and Ericson fluctuations in the unimolecular dis- 
sociation reaction of NO2 molecules into NO + O final 
states. They treat this process as resonance scattering 
within the forr nalism of random-matrix th eory. In a pi- 
oneering work, IStania and Waltheil ( 20051 ) reported the 
first experimental observation of Ericson fluctuations in 
atomic and molecular systems. Quantum chaotic scatter- 
ing was studied in *^Rb in strong crossed magnetic and 
electric fields in an energy regime beyond the ionization 
threshold. The impressive experimental results of pho- 
toexcitation cross sections were later supported by ex- 
act n umerical calculations ( Madrofiero and Buchleitnerl . 
|2005[ ). Very recently, Ericson fiuctuations have also been 
predicted for the inelastic electron scattering cross sec- 
tion on a helium ato m near the double ionization thresh- 
old ( Xu et at l2008[ ). The universality of Ericson fluctu- 
ations in quantum chaotic scattering is, thus, very well 
established. 



VII. VIOLATION OF SYMMETRY OR INVARIANCE 



As mentioned towards the end of Section IIII.Bl the 
theoretical framework of Eqs. (|40l41l42p and ([47]) is quite 
flexible and allows for extensions of the theory that de- 
scribe violation of isospin symmetry, of parity, or of time- 
reversal invariance in CN reactions. These are reviewed 
in turn. In the case of both isospin and parity violation, 
the statistical theory is taken for granted and used as a 
means to analyze the data. The same is true for tests 
of time-reversal invariance in nuclear reactions. But in 
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FIG. 16 Autocorrelation functions (upper part) and Fourier 
coefficients (lower part) with a fit of Eqs. (jTSp to (|77|l of the 
statistical theory. Same as in Fig. [11] but for the Ericson 
regime. The exponential decay in time in the region of over- 
lappi ng resonances {T/D ~ 1.06) is striking. From (|SchafeJ . 
I2OO9I ). 



experiments on induced violation of time-reversal sym- 
metry in microwave billiards, the underlying theoretical 
framework has been thoroughly tested. 



A. Isospin 

Violation of isospin symmetry is treated by replacing 
£jGO E ij^^ ,^ Hamfltonia n H with block struc- 

ture ( Rosenzweig and Porterj . ll960l ) as done in Eq. (1.30). 
We simplify the presentation and consider only two 
classes of states with different isospins Ti and T2 (where 
typically we have T2 = Ti + 1) containing TVi (N2) ele- 
ments, respectively. The generalization to more classes 
with different isospins is straightforward but is not 
needed in practice. The matrix representation of H has 
the form 



(1) 



V 



(2) 



(95) 
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FIG. 17 Upper part: Histogram for the distribution of the 
real part of the matrix element Si2- Lower part: Histogram 
for the distribution of the phase of S12. T he data are ta ken 
in the Ericson regime {T/D « 1.06). From (|Schafeil[2009l ). 



Here (ij'^') are the Hamiltonian matrices acting 

on states that carry the isospin quantum numbers Ti 
and T2, respectively; each is taken from the GOE. The 
isospin-breaking interaction is represented by the rect- 
angular matrix V residing in the non-diagonal blocks. 
The elements of V are due to the Coulomb interaction 
between protons and to other isospin- violating effects, 
see Section I.III.D.l, and are typically small compared 
to those of i/^^-* and H^^\ A similar model is also used 
to describe parity mixing, with and if^^) 

now de- 
noting the Hamiltonian matrices of states of positive and 
negative parity, respectively, and V the induced parity- 
violating nucleon-nucleon interaction. 

The density of states with isospin Ti is usually consid- 
erably larger than that of states with isospin T2 — Ti + 1. 
Isospin mixing has been extensively investigated for two 
cases: (i) In the Ericson regime. For both isospin classes 
the resonances overlap strongly, (ii) For isobaric ana- 
logue resonances. Resonances with isospin T2 are well 
separated. Because of isospin mixing, each acts as a door- 
way for the weakly overlapping resonances with isospin 
Ti. We deal with both cases in turn. 



1. Ericson Regime 

Isospin violation in the Ericson regime (strongly over- 
lapping resonances for both values Ti and T2 of isospin) 
was extensively reviewed in ( Harnev et all . Il986l ) . Since 
then, essential new developments have not occurred ei- 
ther in theory or in experiment. We confine ourselves 
here to a brief summary. The only recent experimen- 



tal information on isospin mixing in very highly excited 
CN is from measurements of 7-ray spectra in heavy- 
ion fusion reactions which we address below. We first 
sketch the modifications of the general framework of Sec- 
tions IIII.Bl IIVI and |V] that are required if two values of 
isospin contribute to the reaction. We then turn to a 
summary of the main experimental results. 

In order to model isospin violation theoretically in the 
Ericson regime one must, in addition to Eq. (|95p for 
the Hamiltonian, also specify the isospin properties of 
the channels. The physical channels labeled (at) carry 
the quantum number t. It is given by the projection of 
total isospin onto some axis and equals half the differ- 
ence of neutron and proton numbers in both fragments 
(projectile plus target). The background matrix S**-*^^ 
that describes scattering in the absence of resonances 
and is given by Eq. ^ as S'-"^ UU"^ , is assumed 
to be diagonal with respect to the physical channels, 



^atM' 



Sab^w 6xp{2z(5^j''}. In other words, in Eq. 



we put both O^^-* and O'^^ equal to unit matrices. Then 
the physical ^-matrix differs from the matrix S"^^ in 
Eq. (|iT|) only by phase factors; these are suppressed in 
what follows. Our assumption neglects direct transitions 
between "mirror" chan nels that are related by neutron- 
proton symmetry; in (|Harnev et all Il986[) it is shown 
that in the Ericson regi me that neglect is ir relevant. 

Although the theory ( Harnev et aLl . [T986[ ) is more gen- 
eral, we focus attention here on the main mechanism of 
isospin mixing. It is due to "internal mixing" induced 
by the matrix elements of V in Eq. ([M)) . Because of 
charge effects, isospin mixing also occurs in the channels 
but is negligible. The coupling matrix elements in 
Eq. (j41[) carry the additional isospin quantum number T 
and so do the transmission coefficients labeled Tqt (we 
deviate here from our standard notation to distinguish 
the transmission coefficients from the isospin quantum 
number). In the physical channels the transmission co- 
efficients T^rp are not diagonal and arc given by project- 
ing (with the help of angular-momentum coupling co- 
efficients {aT\at)) the transmission coefficient t^t onto 
t,t' by r*^ = {aT\at){aT\at')TaT- The autocorrelation 
function of is then given by 

mn 

+ SadSbc X! '^aln ^mnrU ■ (96) 
mn 

Here H„i„ is a 2 x 2 matrix in isospin space. With e = 
E2 — El, the inverse of H is given by 



H" 



JVi+ z + 2iTTe/di 



JV2 + z + 2iTr£/di 



(97) 



The average level spacing for the states with isospin Tm, 
m ~ 1, 2 is denoted by dm, and J^m — J2a '''aT,„ ■ Eqs. 
and ([57)1 generalize Eq. (pT)) for the Ericson regime to the 
case of isospin mixing. The strength of isospin mixing is 
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characterized by a single dimensionless parameter z given 
by 



47r' 



(98) 



The parameter z bears an obvious close analogy to the 
spreading width introduced in Section I.II.G. The ma- 
trix becomes diagonal if z = 0. In the full theory not 
reviewed here, the parameter z comprises both internal 
isospin mixing (via the Coulomb interaction) and exter- 
nal mixing (via the reaction channels). The theory con- 
tains also both, the predictions from the static criterion 
and the dynamic criterion for isospin symmetry breaking: 
The effect is large when either the Coulomb mixing ma- 
trix elements are of the order of the mean level spacing 
(so that z « 1) or when the spreading widths zd\ [zd'i) 
are comparable to the decay widths 2'nNi/di (27rA/'2/d2, 

respec tively). 

In ( Harnev et al\ . Il986l ). the theory was applied to 
a number of experimental examples of isospin mixing 
from which z was determined. The examples can be 
divided into four classes. In the first class, the aver- 
age cross section of a reaction forbidden by an isospin 
selection rule is compared to the average cross section 
of an isospin-allowed reaction. The resulting suppres- 
sion factor / is related to the mixing parameter z and is 
a measure of the average mixing probability of the CN 
levels with isospin T2 with those that have isospin Ti. 
Nuclear reactions in this first category use self-conjugate 
target nuclei (i.e., nuclei with equal neutron and pro- 
ton numbers) and are of the type (c?, a), (d, d'), (a, a'), 
(^Li, a), .... For CN excitation energies close to neu- 
tron threshold one finds suppression factors / around 
0.3-0.5. These decrease to values of only a few percent 
at high excitation energies suggesting that at those en- 
ergies, isospin symmetry is restored. The detection of 
isospin-forbidden dipole radiation from alpha and heavy- 
ion capture reactions supports that observation. While 
the former yields suppression factors / « 0.15 for com- 
pound nuclei in the sc?-shell nuclei at about 14 MeV ex- 
citation energy, in the fusion reaction ^•^C -|- ^^O pop- 
ulating the CN 2^ Si at = 34 MeV, / was found to 
be 5 . 0.05. S i nce th e early experiments ( Harakeh et al\ . 
llQSd ISnoveil Il984l ) several more heavy-ion fusi on reac- 
tions were studied. iKicihska-Habior et al\ ( 19941 ) e.g. in- 
vestigated isospin mixing in the CN ^^Al and ^^S i up to 
excitation energies E,j; sa 63 MeV. Di Pietro et al\ (2001) 



m 24Mg up to « 47 M eV. iKicihska-Habioil (|2005[ ) 
in and ^^Ar, and latelv IWojcik et al\ \2mm also in 
'''^Ti and ^"Zn at E^^ w 50 MeV. (The isospin mixing pa- 
rameter determined in these heavy-ion capture 7-ray 
reactions is directly re lated to the parameter / intro- 
duced in (iHarnev et ah . ,1986 )). The emerging system- 
atics on the very small suppression factors / supports 
the statement that isospin is quite pure at high exci- 
tation energies (several 10 MeV) of the CN. Other ex- 
periments on isospin mixing in the Ericson regime that 
belong to the first category are measurements of the 



isospin-forbidden neutron decay of the giant dipole res- 
onance in medium-heavy (^°Ni) to heavy nuclei (^®Sr, 
89 ^^Zr). One finds suppression factors ranging from 
/ = 0.48 to 0.84 which indicate a fairly sizable isospin 
mixing at nuclear excitation energies of E^ ~ 20 MeV. 
Information on isospin mixing in highly excited com- 
pound nuclei is also deduced from evaporation spectra 
in (a, a'), (p,p'), {p,a') and (a,p') reactions. In short, if 
the cross-section ratio R = Oaa' cfpp' / <^ap'Cfpa' is approx- 
imately equal to unity then the isospin selection rule does 
not play any role. This follows from the Bohr hypoth- 
esis (independence of formation and decay of the CN). 
More generally, an expression for R can be obtained and 
compared with the experimental results for any degree 
of isospin mixin g from the general ized Hauser-Feshbach 
expression (|M1) (jHarnev etai\ . \\m^ ). 

Another highly sensitive test of isospin violation is 
provided by a comparison of cross-section fluctuations 
in the Ericson regime for pairs of isobaric mirror chan- 
nels (channels that are linked by the neutron o pro- 
ton transform a tion) . Such reactions were studied by 
ISimpson et al\ (jl978[ ) and form the second class of nu- 
clear reactions that test isospin violation. An example 
is the reaction ^''N -I- -'^^C leading to the highly excited 
CN ^^Al* which subsequently decays into ^"^Mg + "^H or 
into the mirror channel ^^Na -I- '^He. The most sensitive 
test for isospin violation is provided by the value of the 
cross-correlation function Ctt'(e) at e = for the two 
mirror channels t and t' . If isospin is a good quantum 
number (so that the mixing parameter z = 0), the cross 
sections in the mirror channels are strongly correlated. If 
isospin mixing is complete (so that Z/M2 — >■ 00), the two 
cross sections should fluctuate in an uncorrelated way, 
as does any pair of cross sections pertaining to different 
final states (jEricson and Maver-Kuckiifl Il966l : iRichterl . 
I1974I) . In several pairs of mirror channels the measured 
cross-correlations are significantly larger than for arbi- 
trary pairs of cross sections. 

In the third class of CN reactions sensitive to isospin 
mixing, the shape of the cross-section autocorrelation 
function in the Ericson regime is used as a test. The 
relevant observable is the correlation width F. The 
Lorentzian form of the autocorrelation function is ob- 
tained only for very strong isospin mixing. In the case 
of strict isospin conservation, a superposition of two 
Lorentzians is expected. The theoretical expressions for 
the case of pa rtial isospin symmet ry breaking obtained 
from Eq. (^5)) ( Harney et a^.l [l986l) have been compared 
with very precise cross-section fluctuation data mainly 
from the 32Si(d,a)3°P reaction (Spijkervet, 1978) lead- 
ing to final states with isospins T = and T = 1. 

Charge-exchange reactions like {p,n) or {n,p) that 
populate the isobaric analogue state of the target (see 
Section IVII.A.2P form the fourth class of CN reactions 
sensitive to isospin effects. If isospin were totally con- 
served, such reactions could be viewed as elastic scat- 
tering processes in isospin space. In the Ericson regime, 
the elastic enhancement factor has the value W = 2, see 
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Eq. (|67)) . That same value for W is expected for a charge- 
exchange reaction if isospin is conserved while W = 1 for 
strong isospin mixing. It seems that this interesting test 
has not been used so far. 

We have briefly summarized what is known 
about isospin-symmetry breaking in CN reac- 
tions ( Harney et all Il986l) . The data show that in the 
Ericson regime, isospin-symmetry breaking is neither so 
weak as to be altogether negligible nor so strong as to 
reduce CN scattering to a Hauser-Feshbach situation 
without any reference to the isospin quantum number. 
Isospin mixing is expected to become weaker with in- 



creas ing excitation energy, see (iSokolov and Zelevinskv 



1997 ) and references therein. As shown in (IHarnev et < 
19861 ). data in the Ericson regime can be used to de- 
termine the average Coulomb matrix elements (y^)^/^ 
and spreading widths r\ — 2TTz^d2 or = 27rz^(ii 
for isospin violation in nuclei. While the values of 
the average Coulomb matrix elements vary over many 
orders of magnitude, the values of the spreading widths 
are nearly constant ver sus excitation energ y and mass 
number, see Fig. 15 in ( Harnev et all . Il986l ). This fact 
provides a meaningful consistency check on both theory 
and data analysis. Similarly to the experiments reviewed 
in Section IVI.AI and to the study of sy mmetry breaking 
in the regime of isolated resonances ( Alt et all . 1 19981 : 
iDietz etall . |2006() , further subtle aspects of the theory 
in the regime T ^ d might be tested with the help of 
experiments on two coupled microwave billiards. 



2. Fine Structure of Isobaric Analogue Resonances 

If isospin T were a good quantum number, isospin mul- 
tiplets consisting of degenerate states with fixed T but 
with different z-quantum number ( "isobaric analogue 
states") would exist in nuclei with the same mass num- 
ber A but different neutron and proton numbers. The 
degeneracy is lifted by the isospin-breaking interaction 
(mainly the Coulomb interaction between protons), and 
the energies of the members of a multiplct increase with 
increasing proton number. For proton numbers Z ss 20, 
the energy difference between neighboring members of a 
multiplet is of the order of 10 MeV. In a medium-weight 
nucleus with ground-state isospin Ti, the lowest state 
with next-higher isospin T2 = Ti + 1 (an isobaric ana- 
logue state of the "parent state" , here: the ground state 
of a nucleus with the same mass number but one pro- 
ton replaced by a neutron) typically has an excitation 
energy of several McV. Higher-lying states with isospin 
T2 follow with typical spacings of several 100 keV. These 
states may be unstable against proton decay. The result- 
ing resonances ("isobaric analogue resonances", lARs) 
are then observed in elastic proton scattering. The pro- 
ton channel does not have good isospin and couples to 
both the lARs and the numerous background states with 
isospin Ti. The situation is schematically illustrated in 
Fig. [TSl The parent state is the ground state of the nu- 



cleus z{-^+^)n+i with isospin T2. In its ground state 
the CN nucleus z+i{^+^)n has isospin Ti. The isobaric 
analogue state in that nucleus occurs at an excitation 
energy of several MeV. 
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FIG. 18 Level scheme showing the parent nucleus, the iso- 
baric analogue state, and a state of lower isospin (schematic). 

Low-lying lARs correspond to simple states of the shell 
model. Therefore, their elastic widths (typically sev- 
eral keV) are much larger than those of the complicated 
CN background states. On the other hand, the elastic 
width of the lAR is typically small compared to the av- 
erage spacing between lARs which will thus be consid- 
ered as isolated. The isospin-breaking interaction mixes 
each lAR with CN background resonances. The mixing 
strongly enhances the elastic widths of all CN resonances 
that occur in the vicinity of the lAR. The mechanism is 
similar to that of a doorway state, see Section I.II.G, 
except that now all states are resonances, and that the 
mixing mechanism is rather special. The resulting "fine 
structure" of an lAR is a topic of special interest. While 
it is sometimes possible to apply nuclear-structure the- 
ory to individual lARs, the background states are too nu- 
merous and too highly excited to allow for anything but a 
random-matrix approach. Thus, the theoretical descrip- 
tion uses for the Hamiltonian the matrix (|95p with the 
proviso that the submatrix 7?^^^ is taken from the GOE 
while the submatrix H^"^^ has dimension one. 

lARs as resonances in the CN were discovered by the 
Florida State group whi ch partially resolved a proton s- 
wave lAR in 92Mo(p,p) (|Richard et a/.l . ll964l ). If the CN 
background resonances overlap only weakly, the excita- 
tion curves fluctuate strongly, and the fine structure of 
an lAR can be completely resolved experimentally. That 
is possible mainly in the nuclear \f-2p shell, see Sec- 
tion I. IV. A. The fine structure is investigated in elas- 
tic and inelastic proton scattering and sometimes in the 
(p, n) reaction. The latter is isospin forbidden and gives 
direct evidence for symmetry breaking. Depending on 
the mean level spacing d of the background states, we 
distinguish three cases: almost all of the original proton 
strength is retained by the analogue state (weak mixing). 
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the strength is spread among many states with no state 
being dominant (strong mixing) , and cases between these 
extremes (intermediate mixing). The background states 
have their own (small) proton widths. Thus, the ob- 
served fine structure is the combination of two amplitudes 
which may display interference effects. Special interest 
was shown in the resulting asymmetry ( Robsonl . Il965[) 
that is found in many fine-structure distributions. An- 
other phenomenon in nuclear physics with s imilar mixing 
patterns is that of fission doorways ( Lvnnl . [T969[ ) . 

The observation of fine structure of lARs requires ex- 
cellent energy resolution for the incident proton beam. 
Since most of the fine-structure data were obtained by 
the Triangle Universities' Nuclear Laboratory (TUNL) 
group, we limit discussion of the experimental techniques 
to a brief description of their method. The requirements 
of very good beam-energy resolution (needed to resolve 
the fine structure) and of high beam intensity (needed to 
have good statistics) seem contradictory, particularly be- 
cause of the time-dependent fluctuations in beam energy. 
The method adopted by the TUNL group to resolve the 
resolution-intensity problem uses two beams. One high- 
intensity (H"*") beam is used to perform the experiment, 
while the other (HH"*") beam is used to generate a feed- 
back signal that follows the beam-energy fluctuations. 
This signal generates a voltage difference which is applied 
to the target, thus canceling the time-dependent energy 
fluctuatio ns. The experimenta l details are covered in the 
review bv lBilpuch ^td\ (|1976[ ). 

After early work by the Florida State group on an 
s-wave lAR in 92Mo(p,p) (jRichard et aLl . Il964[ ). later 
meas urements on t his a nalogue with better resolu- 
tion ( Bilpuch et all Il974[ ) provided the fine-structure 
pattern with the largest number of individual states ever, 
seeFig.[THl Prior to these data t he existence of fine struc - 
ture was most c l early shown by iKevworth et. al. I (|1966D : 
IKevworth et al. I (|1968D . see Fig. ^ who used a window- 
less gas target. Their data definitively established the 
essential correctness of the view that in the A w 40 mass 
region, the analogue (doorway) state is mixed into the 
CN background states with a spreading width (see Sec- 
tion LILG) of the order of 10 keV. The bulk of the fine- 
structure data consists of 15 elastic proton scattering ex- 
citation functions on thin solid targets (of order 1 ^g) in 
the mass region 40 < A < 64. There are also some data 
on other channels, including the (p,p'), (p,7), and (p,n) 
reactions. Almost all of t hese r esults are included in the 
review by iBilpuch et al. 1 (I1976D . 

Many data on stren gth func t ions o f lAR display the 
"Robson asymmetry" ( Robsonl . Il965l ) , characterized by 
a dip in the strength function located above the energy 
of the resonance maximum. For simplicity we illustrate 
the origin of that asymmetry for the case of a single open 
channel (the proton channel). That case suffices to dis- 
play the very peculiar nature of the doorway mechanism 
in lARs. The lAR is mixed with the background states 
by two mechanisms: The matrix elements of the Coulomb 
interaction V appearing in the non-diagonal blocks of ex- 
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FIG. 19 Partial widtlis of individual resonances (lower panel) 
and their integral (upper panel) versus proton bombard- 
ing e nergy E„ in the ela stic proton scattering on ^^Mo. 
From (|Bilpuch et all llQT^ ) . 



pression (^5]) . and the elements F^i/ of the matrix pi)) 
connecting lAR and background states. The latter do 
not vanish automatically because isospin is not a good 
quantum number in the proton channel. Moreover, in 
contrast to most other applications of the statistical the- 
ory reviewed in this paper, the elements of the shift ma- 
trix are not small, and the characteristic features of the 
lARs are in fact due to the real part of i^^^. This is be- 
cause lARs occur below or at the Coulomb barrier where 
the energy dependence of the matrix elements Wfj,c{E) 
cannot be neglected. To display the effect most clearly 
we follow Robson, neglect V (that approximation is often 
referred to as "no internal mixing"), and consider only 
mixing due to F^^ ("purely external mixing"). In other 
words, isospin mixing is solely due to the proton channel 
to which both the lAR and the background states are 
coupled. 

We omit the channel index and express the scattering 
function S{E) in terms of the iiT-function, see Eqs. 
and We replace -v/^ttW^ by 7^ and have 



K{E) 



7. 



E-E., 



(99) 



The parameters 7^ and iJ^ of the if-function are exper- 
imentally obtained by a multi-level i?-matrix fit to fine- 
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FIG. 2 Same as in FIG. [191 but for the target nucleus ""Ar. 
From (|Kevworth et aZ.LIl968l '). 



structure data, see Section UlI.CI The object of interest 
is the strength function (7^)/d obtained as an average 
over a number of neighboring resonances, see Figs. [19] 
and [201 We replace in Eq. ([M]) E hy E + il and obtain 
-{1/'k)'^K{E + il) « {il)/d, see Section HvZl The av- 
eraging interval / should contain many CN resonances 
(to reduce the statistical error), but be small compared 
to the spreading width of the lAR defined below (to dis- 
play the resonance enhancement and asymmetry) . In the 
analysis of actual data it may be hard to meet both re- 
quirements. To calculate K{E) we drop V in Eq. ([95]) 
and use Eqs. (|¥T|) and (j^^ keeping the shift function F. 
After a little algebra we obtain 



i^(ii;)=7r^I^,(S-^),,iy, 



(100) 



The matrix B has the same form as expression . The 
indices (i, j) take the values 1 to (for the background 
states) and zero (for the analogue state) while and v 
run from 1 to A'^ as before. Explicitly we have 



„ _ f (E ~ £^,)5f^^ 



{E — Eq + Fqo) 
(1) 



(101) 



We have assumed that the matrix H^nJ + SRF^i/ has been 
diagonalizcd. The resulting eigenvalues are denoted by 



but the notation on the transformed matrix elements 
and 3ffF^o has not been changed. The are random 
Gaussian variables which also appear as arguments of the 
integrals defining the matrix elements 3?F^o- Thus, W^j, 
and SftFi^o are correlated ioT n = ly. The energy of the 
unperturbed lAR is denoted by Eq. The matrix element 
Wo is not random. 

For a common doorway state, the spreading width F^ 
is defined as F^'' — 2nv'^/d, see Eq. (1.26). Here d is the 
mean spacing of the background states and is the aver- 
age squared coupling matrix clement. Using the analogy 
between Eq. poip and the matrix description Eq. (1.24) 
for a doorway state, we define the spreading width of an 
lAR as 



F^ = 27r((5RFo^)2)/d . 



(102) 



This equation shows once again that isospin mixing is 
due to the proton channel. 

The strength function is obtained from Eq. (|100p by 
replacing E hy E + il and assuming d I <^ T^. The 
explicit calculation uses the statistical assumptions men- 
tioned ab ove and may, for instance, be fou nd in Chap- 
ter 13 of ( Mahaux and WeidenmiillMl . ll969| ). The result 



IS 



(E - Eof 



iE-Eo + ^Foi 



(l/4)(Fi)2 



(103) 



Here s^^ is the strength function of the background states 
in the absence of the lAR. The lAR enhances the strength 
function in the vicinity of the analogue state. Since F"^ ^ 
d the resonance is completely mixed with the background 
states and is seen only through the enhancement factor 
in Eq. (|103p . That factor has the shape of an asymmetric 
Lorentzian and approaches the value unity far from the 
resonance. The width of the Forcntzian is given by the 
spreading width. The strength function vanishes at E = 
Eq where there is no isospin mixing. The zero occurs 
above the resonance energy Eq — ^Fqq: Because of the 
Coulomb barrier, the main contribution to the integral 
defining SRi^oo stems from states with energies larger than 
Eq, and 3fiFoo is, therefore, positive. 

In general (several open channels, both external and 
internal mixing) the asymmetry of the strength func- 
tion in channel c in the vicinity of an analogue state 
is reduced. For purp oses o f fittiii g data the expression 
given, for instance, b y Fand ( 1969f ) can be written in the 



form (|Bilpuch et all Il976( ) 



.b9 



2g^9Ac(£o - E) 
= ' (£o-i?)2 + (l/4)(Fi)2 

7c'rV(27r) 
(£o-i?)^ + (l/4)(F^)2 ■ 



(104) 



Here sj!^ is the background strength function in channel 
c, Eq is the resonance energy, Ac is the asymmetry pa- 
rameter in channel c, and 7^ is the total reduced width 
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of the analogue state in channel c, i.e., the sum over the 
fine-structure contributions. 

Eq. (|104|) was originally derived under the assumption 
of strong mixing, 3> d. However, the experimen- 
tal data indicate intermediate or weak mixing in essen- 
tial ly all cases . That case was theoretically considered 
m (|Lane et all \l974 . It was shown that the form of 
Eq. ()104p remains unchanged and that (7^^) /d as given 
by Eq. (|104p must be considered an ensemble average of 
the strength function, all members of the ensemble hav- 
ing the same physical parameters. To analyze the data 
one avoids the need of using energy-averaging intervals 
and fits the accumulated strength to its ensemble average 

Due to the limited number of fine-structure states and 
to the width fluctuations (we recall that the latter follow 
the Porter-Thomas distribution, see Section I.II.D.l), 
the extracted parameters have rather large uncertain- 
ties. The general status of these parameters is as fol- 
lows: The best-fit values for the background strength 
function s[l^ (E) are generally in agreement with those ex- 
pected from systematics. The values of the energy Eq of 
the analogue resonance agree with systematics on shifts 
betw een analogue state and parent state, see ( Janeckel 
I1969D . The values of the reduced widths 7^ can be used 
to calculate the proton spectroscopic factors Sp. The 
latter essentially measure the probability to find a pro- 
ton of fixed angular momentum in the projection of the 
resonance wave function onto the target nucleus. When 
compared with the neutron spectroscopic factors S„ of 
the parent states, the analogue spectroscopic factors are 
significantly lower (30-50 %) than expected, even after 
Coulomb corrections. To the best of our knowledge, that 
discrepancy has never been satisfactorily resolved. The 
statistically significant best-fit values for the asymmetry 
parameter Ac are all negative and usually agree (within a 
factor of 2) with the values predicted bv iRobsonI ( 1965() . 
More detailed analysis indicates that the effects of the 
inelastic channels do not dominate the elastic channel, 
but are not negligible either (as assumed by the Robson 
model). Instead these effects are comparable to that of 
the elastic channel for most analogues. The best-fit val- 
ues of the spreading width agree rather well with the 
Robson model (only external mixing that occurs only in 
the elastic channel). There is evidence that the inelastic 
channels contribute significantly to in some cases. 

The analogue state can decay via several inelastic 
and/or capture-gamma channels and is then a doorway 
common to these channels. Such decay processes pro- 
vide unique information on the analogue state. F or an 
isolated doorway common to two channels c and c' iLand 
( 197l[ ) predicted that the reduced width amplitudes 7c^ 
and 7c'^ of the fine-structure resonances labeled fi should 
be maximally correlated. More precisely, the normal- 
ized linear correlation coefficient /9(7c,7c') defined in 
Eq. (1.29) should be equal to unity and the product 
7cp7c'/i shoul d have the sa.me si gn fo r all fine-structure 
resonances fi. iGraw et al. I (|1974[ ) and lDavis efaj] (Il975l ) 



indirectly showed the expected constancy of the rela- 
tiv e phase. The two p redictions were proved directly 
bv iMitchell efaj] (I1985D . We consider the fragmented 
3/2~ analogue state at proton bombarding energy Ep = 
2.62 MeV in ^^Sc as an example. 

That state has a strong decay to the 2+ first excited 
state in ^'^Ca. In the channel-spin representation, the 
two p-wave proton inelastic decay channels have spins 
s = 3/2 and s = 5/2. These two channels also dis- 
play a well-developed fine structure pattern (see Fig. [2T|) ; 
both show very clearl y the Robson a s ymme try. With the 
method described bv IMitchell et al\ (|l985[) . one can de- 
termine the correlation between the decay amplitudes. 
The measured value p{'ys=3/2,ls=5/2) = 0-93 is in agree- 
ment with the predicted value of unity. The relative 
sign of the decay amplitudes is the same for the 15 
conse cutive resonances that make up th e analogue. To 
quote iMahaux and Weidenmiilleil ( 19791 ) , "isobaric ana- 
logue resonances provide the best-understood example 
of isolated doorway states. " 
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FIG. 21 Fine structure of the analogue resonance with spin 
3/2 in the CN *^Sc. The partial width amplitudes 7 are in- 
dexed by the channel angular momentum and spin. The Rob- 
son a symmetry is clearly displayed. From (jMitchell et al\ . 
Il985l ). 
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B. Parity 

In this Section wc discuss data on parity violation ob- 
tained by low-energy neutron scattering, and their sta- 
tistical analysis. We keep the discussion short since a 
general survey of experiments on parity violation with 
neutro ns and their analysis was given by iMitchell" etali 
( I999I ). A comprehensive review of the experiments and 
analysis of the TRIPLE (Time-Reversal Invariance and 
Parity violation at Low Energi es) collaboration was pre- 
sented bv Mitche ll et all (|200ll) . 

The study of parity violation in nuclei has a long his- 
tory. Following the initial discovery of parity violation 
in beta decay, efforts towards detecting and understand- 
ing the induced parity-violating nucleon-nucleon inter- 
action focused on precise measurements of observables 
indicating parity violation. Most of these involved elec- 
tromagnetic transitions between nuclear states at exci- 
tation energies of a few MeV. In spite of many efforts 
and elegant experimental results, the work has been only 
partially successful. The difficulty is in the theoreti- 
cal analysis - calculating the induced effective parity- 
violating interaction required a more precise knowl- 
edge of the wave functions of the nuclear states in- 
volved than can be attained w it h pre sent-day nuclear 
theory (lAdelberger and Haxtonl. Il985l) . The situation 
changed when Sushkov and FlambaumI ( 1980[ ) predicted 
two enhancement factors which together would lead to 
large parity-violating effects in the CN scattering of low- 
energy neutrons, and when statistical concepts were used 
to analyze the data. 

The two enhancement factors are usually referred to 
as dynamical enhancement and ki nematical enhance- 
ment. Dynamical enhancem ent ( Blin-Stovld . 119601 : 
ISushkov and FlambaumI . Il980l ) arises because in heavy 
nuclei the average spacing of neutron resonances of op- 
posite parity is small and typically 10 eV or so. With 
V the parity-violating part of the nucleon-nucleon in- 
teraction and El and E2 the energies of two levels |1) 
and 1 2) of opposite parity, the mixing of the two states 
is given by {l\V\2)/{Ei — i?2). The mixing obviously 
increases with decreasing spacing {Ei — £'2). The in- 
crease is not inversely proportional to the spacing, how- 
ever, because the complexity of the wave functions of 
the states |1) and |2) also increases with increasing level 
density, reducing the overlap in (1|T^|2). The result 
is an increase of the mixing that is inv ersely propor- 
tional to the square root of the s pacin g (jFrench et all 
Il988al lbl: ISushkov and Flambaun] . Il982l ). With typical 
spacings of states of opposite parity in the ground-state 
domain around 100 keV, the resulting enhancement fac- 
tor is « 10^ and does not change rapidly with excita- 
tion energy. Kinematical enhancement arises because of 
the unequal resonance strength of the states mixed by 
the parity-violating interaction. At low neutron ener- 
gies only resonances with orbital angular momentum 
or 1 are populated. Because of the angular momentum 
barrier, the s-wave resonances have much larger widths 



than the p-wave resonances. Thus, s-wave decay of a p- 
wave resonance with an s-wave resonance admixture is 
enhanced over the regular p-wave decay of that resonance 
by a factor given by the ratio of the two barrier penetra- 
bilities - approximately {kR)~^. Here R is the nuclear 
radius and k the wave number. The factor (kR)~^ is 
strongly energy dependent, significant enhancement oc- 
curs only near neutron threshold. The product of the 
two enhancement factors is about 10^. Since the weak 
interaction is approximately a factor of lO'^ smaller than 
the strong interaction, the total enhancement leads to 
expected parity-vio lating effects of the order of per cent. 
That prediction of ISushkov and FlambaumI (|l980[ ) was 
confirmed shortly afterwards at Dubna (Alfimenkov et 
ai, 1982, 1983). The results at Dubna were extremely 
interesting. However, there were both experimental and 
theoretical limitations. 

The experimental limitations at the Dubna facility 
were severe: the neutron flux dropped dramatically above 
a few tens of eV, limiting the experiments on parity viola- 
tion to neutron energies below 20 or 30 eV. Subsequently 
the TRIPLE collaboration was formed to extend the ex- 
periments on parity violation to higher energies. These 
experiments were performed at the Los Alamos Neutron 
Science Center (LANSCE). A moderated and collimated 
beam of neutrons is produced by spallation. A beam of 
longitudinally polarized protons is produced by scatter- 
ing from a longitudinally polarized proton target. The 
neutron spin direction is reversed by a system of mag- 
netic fields. The neutrons pass through the target and 
are detected in a highly segmented system of liquid de- 
tectors located at a distance of approximately 57 meters. 
CN resonance energies are determined by the time-of- 
flight method. A sample result is shown in Fig. [22| The 
de tails of the experimen tal system are given in the review 
bv lMitchelf et al. I (I2001D . 




FIG. 22 Neutron transmission spectra for two helicity states 
near the 63 eV resonance in ^^^U. The resonance appears as a 
dip in the transmission curve. The transmission at the reso- 
nance differs significantly for the two helicity states and par- 
ity v iolation is apparent by inspection. From (|Mitchell et all 
I1999I ). 
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Although the observation of large parity violation in 
neutron resonances was certainly very impressive, the 
Dubna results were first considered of only anecdotal in- 
terest, since the resonance wave functions were too com- 
plicated to be theoretically accessible. This problem was 
overcome by the statistical approach. Since s~wave neu- 
tron resonances obey GOE statistics, it is safe to assume 
that p-wave resonances do, too. Then the matrix ele- 
ments of the parity-violating interaction connecting s- 
wave and p-wave resonances have a Gaussian distribu- 
tion with mean value zero. The variance of the distribu- 
tion (or the mean-squared matrix element for parity 
violation) can be determined from a set of p-wave res- 
onances with parity violation in a given nuclide. This 
is described below. However, is not a good mea- 
sure for the strength of the parity-violating interaction, 
since it decreases rapidly with increasing complexity of 
the wave functions and, thus, with increasing mean level 
density. A measure that is roughly independent of exci- 
tation energy and mass number is the spreading width 
r-l- = 27r^/d (see Sections I.II.G and IVII.Ap . The con- 
vention is to adopt for d the mean spacing of the s-wave 
resonances. The expected size of the weak spreading 
width can be estimated from the spreading width for 
the strong interaction (which is experimentally of order 
MeV) and by adopting for the ratio of the square of the 
strength of the weak interaction to the strength of the 
strong interaction the value 10"^^. Thus one expects 
r-^ « 10"^ eV. 

The observable for parity violation is the longitudinal 
asymmetry (often simply referred to as "asymmetry" ) 



respectively. The matrix ([95)1 takes the form 



P 



(105) 



where cr^ is the total p-wave cross section for neutrons 
with helicities ±. Clearly we have P = if parity is con- 
served. Asymmetries for a set of resonances were deter- 
mined separately for each run of approximately 30 min- 
utes and a histogram created for each nuclide and each 
resonance measured. The mean of this histogram was the 

value adopted for P. 

We illustrate the analysis that determines and the 
weak spreading width by considering a target nucleus 
with spin / = and positive parity. This case illustrates 
most of the principles involved. The s~wave resonances 
have spin and parity 1/2"*" and the p-wave resonances 
l/2~ or 3/2~. Only the l/2~ resonances are consid- 
ered as only these are mixed with the 1/2+ resonances 
by the parity-violating interaction. We use the formal- 
ism of Section IIII.BI We neglect direct reactions, use 
the scattering matrix as given by Eqs. dH]) and (|^ . 
and replace the Hamiltonian i^'^^E -^^ gyj by an 
expression of the form (^5]) . The upper indices 1 and 2 
now stand for states with positive and negative parity, 
respectively, and V denotes the induced parity-violating 
nucleon-nucleon interaction. We use the diagonal repre- 
sentation of _ff and Jf(^) and denote the eigenvalues by 



i?^^^ and Elf' and the eigenfunctions by \lfi) and \2i/), 



p(2) 



(i/i|y|2p) 



{2<j\VM Ef^5p„ 



(106) 



In the diagonal representation of iJ'^' and i/'^', the par- 
tial width amplitudes of the states with positive and neg- 
ative parity are denoted by g'lP and g^ff' , respectively. We 
assume that the p-wave resonances are isolated and fo- 
cus attention on a single one. We take the bombarding 
energy E in the center of that resonance, E = Ep . We 
neglect the total width of that resonance and of the ad- 
mixed s-wave resonances since in all cases investigated so 
far, the spacings \E^^'^ — E^ p^ \ are larg e compared to the 
total widths. This yields (iBunakov and Gudkovi . il98H 
ISushkov and Flambaun] . Il980[ ) 



{MV\2p) 



E\ 



(2) 



(1) (2) 
-p(n) 



(107) 



where T 



in) 



{g^p'Y is the the partial width for neu- 



tron decay of the p-wave resonance with label p. The 

ratio 5^i^''<7p'^Vr|)"'' = dji^ /g^^^ contains the kinematical 
enhancement factor, and the first term on the right-hand 
side, the dynamical enhancement factor. 

For spin-zero target nuclei, the resonance parameters 
are usually known. For the s-wave resonances, the infor- 
mation is available from previous work on s-wave neu- 
tron scattering. For the p-wave resonances, most of 
the information was obtained in the framework of the 
TRIPLE experiments. In practice one may assign the 
spin value 1/2 to a p-wave resonance by the presence 
of parity violation. Unfortunately one cannot determine 
the individual matrix elements (l/i|y|2p) since there are 
too few equations and too many unknowns. But us- 
ing the fact that the matrix elements (l^|y|2p) have a 
Gaussian distribution with mean value zero and a sec- 
ond moment given by u^, we write Eq. ()107|) in the 
form Pp = ^^^(l/^l^l^p); with coefficients App = 



{2/{eI^'> - Ej,^'')){gl!>/g}f>). Then, Pp is a linear combi- 
nation of equally distributed Gaussian random variables 
and, therefore, is itself a Gaussian random variable with 
mean value zero. The variance of Pp with respect to both 
fi and the ensemble generated by a sequence of runs is 
given by v^, where = {l/N)J2pp^l,p and where 
N is the number of p-wave resonances. It follows that 



,(l)/.(2) 



— var(P^, 



(108) 



Eq. (|108p is the central result of the statistical approach. 
It yields v-^ from the data in spite of the fact that the 
signs of the partial width amplitudes in Eq. (|107|) are 
usually not known. The analysis is more difficult for 
target nuclei with non-zero spin values. Moreover, usu- 
ally some but not all spectroscopic information is avail- 
able. Suitable methods of analysis were developed for all 
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such cases ( Mitchell et all l2Q0l[) . but the spirit of the 
approach is the same. It yields the values of and of 
although the wave functions of the individual nuclear 
states are not known. 

For the actual determination of i;^, a maximum like- 
lihood approach was adopted. The value of Pp is a re- 
alization of a random variable. For a number of inde- 
pendent resonances the likelihood function is the prod- 
uct of the functions for the individual resonances. One 
inserts the values of the experimental asymmetries and 
their errors, determines the spectroscopic term A from 
the known resonance parameters, and calculates the like- 
lihood function. The location of the maximum gives the 
most likely value wl of the root mean square matrix el- 
ement V = \f%P- . A maximum-likelihood plot for ^ssy 
is shown in Fig. [231 The likelihood function is found to 
have a well-defined maximum so that v is rather well de- 
termined. When that analysis was applied to ^'^^Th, the 
mean value of (l/i|F|2p) was found to differ from zero, 
and a Gaussian with non-zero mean value was required 
to fit the data - all parity violations had the same sign. 
This result raised serious questions about the statistical 
approach and led to much theoretical activity. However, 
the anomaly only occurred in Thorium, and subsequent 
studies at higher energie s in ^^^Th revealed pa rity vio- 
lations of opposite sign (jSharapov et all |2000[ ). It was 
concluded that this "sign effect " was due to a local door- 
way state (see Section I.II.G). 




V (meV) 

FIG. 23 Plot of the maximum-likelihood function L(v) versus 
the root-mean-square matrix element v for '^^^ U. The spins of 
all th e p-wave resonances are known. From ([Mitchell et all 
119991 ). 



Parity violation was studied in 20 nuclides - especially 
in those regions of mass number A where the p-wave 
strength function has maxima. (The strength function is 
defined and discussed in Sections [IV. Dl and IVI. A[) . Max- 
ima occur near A = 238, i.e., for 232rp]^ 238^ (-j^^x- 
imum of the 4p strength function) and near A = 100 
where most of the other nuclides were studied (maxi- 
mum of the 3p strength function). Parity violation was 



observed in all but one of these targets (^^Nb). For 15 
nuclei sufficiently many resonances with parity violation 
were observed to determine the weak spreading width; for 
the other nuclides only very approximate values or limits 
could be determined. The values of lie around 10~^ 
eV as expected and are approximately independent of A, 
with some indications of local fluctuations. Such fluc- 
tuations have been observed in th e spreading width f or 
isospm mixing (Section[VlLXland ([Harnev et all\\%%<^ \. 

In summary: Except for the Thorium anomaly, the 
data are consistent with the statistical model. With 
the help of a statistical analysis, it is possible to deter- 
mine the root-mean-square matrix element v for par- 
ity violation and the weak spreading width F-'' without 
knowledge of the wave functions of individual nuclear 
states. The values found for are consistent with ex- 
pect ations based on t he str ength of the weak interaction, 
see ( Tomsovic et all [2000l ) for an analysis. For lack of 
space we have n ot discussed experi ments on parity vio- 
lation in fission ( Kotzle et al ]. [2000[) . 



C. Time Reversal 

Time-reversal (T) invariance implies symmetry of the 
^-matrix, Sah = Sha and, hence, detailed balance, 
\Sa\i^ = \Sha? ■ Tests of T invariance compare resonance 
scattering cross sections for the two-fragment reactions 
a —i' b and 6 — !• a at the same center-of-mass energy and 
aim at establishing an upper bound on the strength of 
the T-invariance violating interaction in nuclei. Such ex- 
periments have been performed both for isolated and for 
strongly overlapping resonances more than twenty years 
ago and are well documented in the literature. Thus, we 
can be brief. 

To test detailed balance, one compares the reaction 
rates Ai + A2 ^ Bi + B2. The fragments ^1,^2, -81,-82 
are in their ground states and unpolarized. A difference 
in the rates a — >■ 6 and 6 — >■ a indicates a v iolation of 7~ 
invariance. Detailed balance was tested bv [Driller et all 
([1979[) in the reactions ^^Al + p ^ ^^Mg + a populating 
an isolated J'^ — 2+ resonance at an excitation energy 
= 12.901 MeV in the CN ^^Si. For a reaction through 
an isolated resonan ce, detailed balance would no rmally 
hold automatically ([Henlev and Jacobsohnl [l959l ) . The 
present case is different because for the fragmenta- 
tion ^^Al -|- p, the partial waves with angular momenta 
I = and I = 2 and channel spin s = 2 inter- 
fere ( Pearson and Richteil . [19751 ) . The cross sections for 
the reactions a ^ b and b ^ a were found to agree within 
S = 0.0025±0.0192%. This resuh is consistent with 5 = 
and, thus, with T invariance. 

[Bunakov and Weidenmiilleil ( 1989[ ) pointed out that in 
detailed-balance tests that use two close-lying CN res- 
onances in the regime T ^ d, large enhancement fac- 
tors amounting to several ord ers of magnitude for T- 
invariance violation may arise. [Mitchell et all ([1993[) in- 
vestigated specific experimental possibilities and showed 
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that the difference S of the two reaction cross sections 
depends sensitively on energy, angle, and on the param- 
eters of both resonances and may vary by many orders 
of magnitude. These theoretical predictio ns have been 
parti ally tested experimentally in billiards ( Dietz et all . 
120071 ). 

In the most precise test of detailed balance in CN 
reactions so far, the reactions 2''Al(p, Q;o)^^Mg and 
2^Mg(Q , iPn)^^Al were comp ared in the Ericson regime 
r ^ d (iBlanke et at . 1983). As in a predecessor of this 
experiment ( von Witsch et fflEl . ll968() . both reaction rates 
were measured at a scattering angle of « 180° and were 
normalized at a suitable cross-section maximum. The 
results were then compared at a cross-section minimum. 
This was done in order to maximize the sensitivity for 7~- 
invariancc violation. The result is shown in Fig. 1241 The 
measured relative differential cross sections agree within 
the experimental uncertainty d = ±0.51 % and are, thus, 
consistent with 7~-invariance. From this result an upper 
bound ^ < 5 • 10^"* (80 % confidence) for a possible T- 
noninvariant amplitude in the CN reaction was derived. 
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FIG. 24 Normalized cross sections for the reaction 
^^Al(p, ao)^''Mg (solid line) and for the inverse reaction 
^*Mg(Q,po)^^Al (open circles) near and at a deep minimum. 
The t wo cross sections are equal, and detailed balance holds. 
From (|Blanke et alUWSA ). 

At the time of the experiment, there was no adequate 
theoretical framework to interpret that upper bound. 
The development of the st atistical theory mad e such 
an interpretation possible ( Boose et all Il986allb[ ) and 
yielded upper bounds for both, the strength of and the 
spreading width F-'' for the T-invariance violating part 
of the nuclear Hamiltonian. Both bounds compare favor- 
ably with those derived from an anal ysis of spectr a l fluc- 
tuatio ns (sec Section I. III. D. 4 and (jFrench et all Il985l 
ll988bD '). As in the cases of violation of isospin and of 
parity in Sections IVII. AI and IVII.Bl the fundamental pa- 
rameter is the spreading width, and the upper bound on 
that quantity is F-^ < 9 • 10"^ eV. The upper bouiid fron i 
the detailed balance experiment of ( Blanke et all Il983l ) 
has been improved by two orders of magnitude in an 
experiment on polarized neutro n transmission throug h 
nuclear-spin-aligned Holmium ( Huffmann et all . Il997 ). 



The measurements test reciprocity and could possibly im- 
proved by another order of magnitude with more intense 
neutron beams at new spallation sou rces and with other 
targets ( Barabanov and Bedal l2005l) . The expected up- 
per bounds for the strength of possible parity-conserving, 
time-reversal-violating interactions will, however, still 
be several orders of magnitude larger than the ones pro- 
vided by the upper limit on the electric dipole moment 
of the neut ron for parity and time-reversal violating in- 
teractions (|Baker et all 120061: iHarrisj, |2007[) . 

To include violation of T-invariance in the statistical 
theory, the T-invariant Hamiltonian ensemble TJ^ioe 
the right-hand side of Eq. (|47)) is generalized, see Sec- 
tion I.III.D.4, 



^GOE ^ ^ ^ 



(gGOE + ,"^). (109) 

v/1 + (l/iV)a2 ^ VN ' ^ ' 



The independent elements of the real, antisymmetric, 
A'^-dimensional matrix A are uncorrclated Gaussian- 
distributed random variables with zero mean values and 
second moments given by A^^Apcr ~ [X^ /N){5^p5i,a — 
5p„5i,p). The parameter a measures the strength of T- 
invariance violation. As explained in more detail in Sec- 
tion I.III.D.4, the normalization factor N~^/^ is chosen 
so that significant invariance violation on the scale of 
the mean level spacing occurs for a ^ 1. With the re- 
placement (|109p in Eq. (|T7)) . the calculation of measures 
for T-invariancc violation in the statistical theory is a 
formidable task. In the Ericson regime, iBoose et al\ 
( 1986al lb[) solved the problem by combining a perturba- 
tive treatment of A with an asymptotic expansion in pow- 
ers of d/ r. A treatment valid for all values of F/d was 
given in (jGerland and Weidenmiille"il . ll996l : |Pluhar et all 
Il995l ) and further developed for the analysis of scatter- 
ing data on microwa ve resonator s with i nduced violation 
of T-invariance in (jPietz et all l2009a[) . These papers 
were based on progress in underst anding the GOE — » 
GUE cross over transition in spe ctra ( Altland et a/.l . [r992l 
Il993l ). In (jPietz et all l2009al) . the parameter a IS re- 
placed by TT^. The resulting expressions for measures of 
T-invariance violation are complex and not reproduced 
here. We confine ourselves to a discussion of the results. 

In Sections IVI.AI and IVI.BI it was shown that many 
properties of chaotic scattering can be studied with the 
help of microwave resonators. That statement applies 
also to the violation of T-invariance. In a flat chaotic mi- 
crowave resonator — a quantum billiard — shown schemat- 
ically in Fig. 1141 T-invariance violation can be induced 
by placing a ferrite (magnetized by an external fleld) into 
the resonator. The spins within the ferrite process with 
their Larmor frequency about the magnetic field. This 
induces a chirality into the system. The magnetic-field 
component of the radio frequency (rf) field in the res- 
onator can be split into two circularly polarized fields 
rotating in opposite directions with regard to that static 
magnetic field. The component that has the same rota- 
tional direction and frequency as the rotating spins, is 
attenuated by the ferrite. This causes T-invariance vi- 
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olation. The strongest effect is expected to occur when 
Larmor frequency and rf frequency coincide. Connecting 
the resonator to two antennas, one defines a scattering 
system. The violation of T-invariance then causes the 
scattering amphtudes S12 and 6*21 to differ. 

Experiments on induced 7~-reversal invariance viola- 
tion in microwave billiards offer two advantages. First, 
by a measurement of phases and amplitudes of the re- 
flected and transmitted rf signals it is possible to test the 
reciprocity relation Sat = Sta while experiments with 
nuclear reactions typically test only the weaker detailed- 
balance relation |S'a6p = jS'haP- Second, such experi- 
ments offer the unique chance of a stringent and detailed 
test of the statistical theory that is otherwise often taken 
for granted and used to analyse data on CN reactions. 

Induced violation of T-invariance in microwave bil- 
liards has so far been stu died in two s catter ing exper- 
iments. In the first one (jPietz et all |2007[) . a vector 
network analyzer was used to measure magnitudes and 
phases o f 5'-matrix elements in a fully chaotic "annular" 
billiard (jPembowski et all I2OOOI: iHofferbert et al. I. I2005D 
in the regime T <^ d. By interchanging input and output 
antennas, both S12 and 521 were measured. For all eight 
isolated resonances (singlets) that were investigated, the 
complex element 5*12 agr ees with ,S'2i , i.e., reciprocit y 
holds, in agreement with (jHenlev and Jacobsohnlll959t) . 
For the three pairs of partially overlapping resonances 
(doublets) that were studied reciprocity was found to be 
violated. The dependence of the T- violating matrix el- 
ements of the effective Hamiltonian for the microwave 
billiard on the magnetization of the ferrite could be de- 
termined with the help of E qs. (|109p and (|41I47)) . 

In the second experiment ( Dietz et aLl . l2009al) . induced 
violation of T-invariance was investigated in the regime 
of weakly overlapping resonances. A small cylindrical 
ferrite was placed within the fully chaotic tilted-stadium 
billiard shown on the inset of Fig. [TOl The ferrite was 
magnetized by an external magnetic field. Again the ele- 
ments 5*12 and S'21 of the complex-valued S'-matrix were 
measured versus resonance frequency. Figure [55] shows 
that magnitude and phase of the S'-matrix elements fluc- 
tuate strongly, and that reciprocity is violated. The value 
of the normalized cross-correlation coefficient 



Ccross(O) — 5R- 



\Sl2\^ |S2l|2]l/2 



(110) 



serves as a measure of the strength of T-invariance vi- 
olation. T-invariance holds (is completely violated) for 
Ccross(O) = 1 (Ccross(O) = 0, respectively). The upper 
panel of Fig. [26] shows that Ccross(O) depends strongly on 
frequency. Complete violation of T invariance is never 
attained. The lower panel shows the value of the param- 
eter ^ = a/TT for T-invariance violation deduced from 
the data. 

The data were used for a tho rough test of the under- 
lying theory ( Dietz et'ol} . l2009al) . The parameters of the 
theory (the transmission coefficients Ti and T2 in the two 
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FIG. 25 Magnitude and phase of (solid lines) and of 
S'21 (dashed lines) measured in the frequency interval from 
16 to 16.5 GHz for a fi xed magnetization of the ferrite. 
From (jPietz et o/.l . l2009al ). 
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FIG. 26 Upper panel: Values of Ccross(O) calculated from the 
fluctuating matrix elements 5i2 and S'21 of Fig. 1251 Lower 
panel: The par ameter ^ v ersus frequency as deduced from 
the data. From ([Dietz et a/.. . .2009a ). 



antenna channels, a parameter describing Ohmic absorp- 
tion by the walls of resonator and ferrite, and ^) were fit- 
ted to the Fourier-transformed S'-matrix elements. As in 
Section [VI.Al a goodness-of-fit test was used to establish 
the quality of the fit, with excellent results. Moreover, 
values of the elastic enhancement factor versus frequency 
were predicted correctly by the theory without further fit 
parameters. We recall that for T-invariant systems, that 
factor takes values between 2 and 3, see the remarks be- 
low Eq. ([65]). Experimentally, values well below 2 were 
found in some frequency intervals. These are possible 
only if T invariance is violated. 



VIII. SUMMARY AND CONCLUSIONS 

The Bohr assumption on independene of formation and 
decay of the CN was the first of several insightful conjec- 
tures concerning compound-nucleus scattering. It laid 
the ground for all later developments and was followed 
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by the Hauser-Fcshbach formula, the Weisskopf estimate, 
Ericson's idea of random fluctuations of nuclear cross sec- 
tions, and the generalization of the Hauser-Feshbach for- 
mula for the case of direct reactions. By way of justifica- 
tion, all of these developments referred to the fact that 
the resonances that dominate compound-nucleus reac- 
tions, have statistical properties. The work anticipated 
general features of quantum chaotic scattering discovered 
only many years later. 

The challenge to actually derive these conjectures from 
the statistical properties of resonances was taken up 
early. But it took several decades and the efforts of 
many people until a comprehensive theory of compound- 
nucleus scattering was established. The theory is based 
on a description of the statistics of resonances in terms 
of Wigner's random matrices. That description applies 
generically to resonances in chaotic quantum systems. As 
a result, the statistical theory of nuclear reactions is, at 
the same time, a generic random-matrix theory of quan- 
tum chaotic scattering. For systems with few degrees 
of freedom, it competes with the semiclassical approach 
to chaotic scattering. The latter incorporates system- 
specific features in the form of short periodic orbits. In 
the statistical theory, such features enter as input param- 
eters. Unfortunately, the semiclassical approach has so 
far not been extended to many-body systems. 

The development of the theory posed essentially two 
major challenges: (i) How to incorporate random-matrix 
theory into scattering theory and (ii) How to calculate 
moments and correlation functions of S'-matrix elements 
from that input. The first problem required a formula- 
tion of scattering theory in terms of the Hamiltonian gov- 
erning the quasibound states that turn into resonances 
as the coupling to the channels is taken into account. 
The second problem was solved with the help of meth- 
ods adopted from quantum field theory. Such methods 
have found wide application in condensed-matter the- 
ory. The connection to this area of physics is not one 
of technicalities only. It actually connects the statistical 
theory of nuclei with the statistical mechanics of many- 
body systems. Conceptually, it shows how a separation of 
scales is achieved in the limit of large matrix dimension 
TV: Universal features govern the system on an energy 
scale given by the mean level spacing d. The unphysical 
properties of random-matrix theory (i.e., the shape of 
the average spectrum) matter on the scale Nd. 

The ensuing statistical theory is a complete theory that 
uses the minimum number of input parameters and has 
predictive power. As is typical for many applications 
of random-matrix theory, the statistical theory predicts 
fluctuations in terms of mean values. The latter com- 
prise the values of the energy-averaged elements of the 
scattering matrix and, in the case of correlation func- 
tions, the mean level spacing d. It is here that system- 
specific features enter. Examples are the optical model 
and the strength function which reflect properties of the 
nuclear shell model. The theory predicts the values of 
moments and correlation functions of S'-matrix elements. 



These determine mean values and fluctuation properties 
of cross sections and other observables. The resulting 
expressions vindicate the early conjectures, define the 
limits of their applicability, and yield expressions that 
hold under more general circumstances. The S'-matrix 
autocorrelation function is a case in point. In nuclei, 
the range of energies where the statistical theory applies 
is limited by the underlying assumption that random- 
matrix theory correctly describes the statistical proper- 
ties of resonances. That is true only when the nuclear 
equilibration time is shorter than the average life time of 
the compound nucleus and holds for bombarding energies 
up to 20 MeV or so. Beyond that range, precompound 
processes modify the reaction dynamics. 

The theory has been the object of stringent tests, 
largely performed with the help of microwave billiards. 
These have been extremely successful. There is no rea- 
son to doubt that the theory adequately accounts for all 
aspects of chaotic scattering even though one aspect of 
the theory has received little attention so far and has 
not been seriously tested: The unitary transformation 
in channel space that takes account of direct reactions. 
The reason is that direct reactions and CN processes are 
almost mutually exclusive. When one is important, the 
other one typically is not, and vice versa. 

The theory has found numerous applications both 
within and outside the field of nuclear physics. These 
have only partly been reviewed. We have paid spe- 
cial attention to violations of parity, of isospin, and of 
time-reversal invariance. In the first two cases, the the- 
ory allows for the determination of the strength of the 
symmetry-violating interaction. In the case of time re- 
versal, it yields an upper bound on that strength. 

The calculation of moments and correlation functions 
of the S-matrix from the statistical theory is not as com- 
plete as one may wish. It would perhaps be unrealistic to 
expect complete knowledge of all moments and all corre- 
lation functions. The theoretical effort grows immensely 
with increasing order of such expressions. Moreover, mo- 
ments and correlation functions of low order only can 
reliably be determined from the finite range of data typi- 
cally available in experiments. But it would be desirable 
to have theoretical expressions for cross-section correla- 
tion functions even though approximate expressions are 
available. That poses a continuing challenge for theorists. 
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